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Deconvolution  of  a  single  convolution  equation  is 
usually  an  ill-posed  problem.  This  has  been  sufficiently 
illustrated  in  the  literature.  The  shortcomings  of  linear 
and  of  non-linear  deconvolution  methods  can  be  found,  for 
instance,  in  the  very  clear  review  paper  .  Advances  in 
the  theory  of  holomorphic  functions  of  several  complex 
variables  led  Berenstein,  Taylor  and  Yger  to  realize  that 
systems  of  convolution  equations  could  be  deconvolved 
exactly,  thus  avoiding  the  above  ill-posedness.  Their 
preliminary  papers^X^r^T^eventually  led  to  this  project.  The 
practical  interest  of  this  observation  is  that  whenever  such 
a  set  of  convolution  equations  represents  a  set  of  physically 
realizable  devices  (e.g.  transducers,  sensors)  then  one  has, 
by  use  of  a  digitally  implemented  inverse,  essentially  an 

arbitrary  bandwidth  device .  . Let  us  make  these  ideas  exact  by 

. 

recalling  the  following  definitions  and  deconvolu  r-  a  scheme. 

Let  be  a  collection  of  distributions  of 

compact  support  in  Rn  ("convolvers’').  Their  Fourier 

A 

transforms  y ^  <  z )  are  holomorphic  functions  of  n  complex 
variable  in  the  Paley-Wiener  class.  We  say  that 

i  n 

are  strongly  coprime  if  there  are  constants  e>o,  C>0  such 
that  for  every  7  e  c"  : 


1 


*  ^j(z^  -E  exp(~C  (log(l  +  |*|)  +  |  Im  z|  }  ) . 


In  that  case,  one  can  prove  the  existence  of 
"deconvolvers"  *  That  is,  distributions  of 

compact  support  such  that 

u  *  v,  + _ +  y  *  v  =6. 

11  mm 


Once  this  v  have  been  explicitly  found,  and  herein  lies  the 
difficulty  we  set  out  to  solve,  one  has  the  following  scheme 
to  determine  an  unknown  signal  <{> 


♦ 


sampling 

— © — 

&  noise 

— 0 - 

— © - 


V 


v2*  — © 

V  ir  - 
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A  very  thorough  discussion  of  the  meaning  of  this 
scheme,  and  its  implementation  in  the  1-d  case,  can  be  found 
in  the  final  report  to  a  previous  ARO  supported  project  [4], 

In  this  project  we  set  ourselves  the  following  tasks: 

(i)  find  explicit  and,  relatively  easy,  formulas  for 

the  deconvolvers  in  the  1-d  and  2-d 

situations. 

(ii) show  that  in  the  above  deconvolution  scheme  the 
overall  system,  in  the  presence  of  noise,  is 
entirely  dependent  on  the  sampling  rate  and  noise 
characteristics,  i.e.  no  "inherent"  bandwidth 


2 


limitations.  This  means  to  show  that  the  above 


block  diagram  behaves  like 


4 — ►  overage  over  smell  delector  =  4>, 


(iii)  construct  1-d  simulations  of  the  deconvolution 
scheme  to  show  its  feasibility  beyond  purely 
theoretical  considerations. 

These  three  objectives  have  all  been  accomplished  and 
documented  as  it  will  be  explained  now.  Furthermore,  the 
payoff  of  the  powerful  techniques  developed  for  these 
questions  has  also  taken  place  in  two  other  areas:  explicit 
solution  of  algebraic  equations  and  their  complexity,  and  the 
local  Pompeiu  problem,  these  will  be  discussed  below. 

First  of  all,  we  have  two  schemes  to  find 
deconvolvers.  An  analytic  one  is  presented  in  Chapter  1 
where  everything  depends  on  the  following  identity 
illustrated  for  the  2-d  case  with  three  convolvers: 

viz)  -  Y  - 

rj  _  'Vv w  _j 
1  V'l  •  '  v<2 

*  * 

J(0  *  Jacobian  determinant  of  C 


ij(Oi2(C)  i3«) 
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Here  u  can  be  taken  as  a  mollifier  of  small  support 
and  the  estimate  $  of  the  signal  $  will  then  be 

<J>  =  Vj*  iJjj  +  ...  +  v3*  ^3, 

where  j  =  (  p  *  +  noise)  *u  =  "smoothened"  data.  The 

resolution  level  wilx  just  correspond  to  the  sampling  rate. 
This  chapter  has  appeared  in  print  in  [5]. 

In  other  words,  we  can  say  we  construct  analytically 
deconvolvers  ,  once  for  all,  and  then  we  implement  a 
numerical  integration  scheme  (choice  of  u  and  sampling) . 

There  is  a  second  method  to  find  numerical 
deconvolvers  directly,  explained  by  Dr.  Taylor  in  the  recent 
West  Point  ARO  Conference  on  Computational  and  Applied 
Mathematics.  In  this  second  method,  for  a  fixed  grid 
(resolution)  we  compute  the  deconvolvers  numerically  and  the 
existence  of  the  analytic  deconvolvers  guarantees  that  when 
the  mesh  of  the  grid  decreases  and  the  noise  decrease  the 
corresponding  estimates  $  converge  to  4>  We  do  not  yet 
have  sufficient  numerical  data  to  compare  these  two  methods. 

In  Chapter  2,  we  present  an  evaluation  of  the  effect 
of  noise  on  deconvolution,  background  limited  detectors,  and 
amplifier  limited  detectors  are  considered,  a  criterion  is 
given  to  evaluate  and  design  different  systems  of  convolvers. 
We  prove  that  deconvolution  of  strongly  coprime  systems 
behave  well  with  respect  to  noise.  This  is  crucial  for 
practicality  and  for  the  actual  engineering  design  of  such 
systems.  This  chapter  constitutes  the  paper  [6]  which  has 
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been  submitted  by  invitation  and  accepted  for  publication  in 
the  forthcoming  IEEE  Proceedings  in  Multidimensional  Signal 
Processing.  It  is  an  expansion  of  part  of  the  Ph.D.  thesis 
of  E.  Vincent  Patrick  written  under  the  direction  of  Dr. 
Berenstein. 

In  some  applications,  we  find  ourselves  with  occluded 
near  regions  where  data  could  be  collected  by  sensors  but  no 
deconvolution  in  the  above  sense  performed.  This  is  the  so 
called  Local  Pompeiu  problem.  Let  us  say  a  signal  f  is 
collected  by  averaging  its  values  over  any  square  of  side  a, 
but  that  this  can  be  done  only  inside  some  disk  of  radius  R. 
What  is  the  precise  relation  between  R  and  a  (if  there  is 
any)  such  that  the  data  collected  determines  f  (this  would 
allow  to  decide  whether  f  is  what  we  are  looking  for  or  not) . 
Furthermore,  can  we  reconstruct  f  from  this  data? 
Surprisingly,  the  answer  to  both  questions  is  yes,  at  least 

when  R  >  a~  for  the  first  problem  and  R  >  2  ^  a  for  the 

2 

second.  This  is  the  content  of  Chapters  3  and  4, 
respectively.  It  is  not  necessary  to  emphasize  the 
importance  of  this  work  in  the  problem  of  Automated  Target 
Recognition.  Chapter  3  has  been  expanded  and  appeared  in 
print  as  [8]  and  Chapter  4,  [9],  has  been  accepted  for 
publication. 

Let  us  conclude  this  introduction  with  the  observation 
that  the  formulas  from  Chapter  1  have  had  an  unexpected 
impact  in  problems  about  systems  of  algebraic  equations, 
their  complexity  and  properties  of  the  solutions  of  the 
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algebraic  Bezout  equation.  In  this  sense,  this  paper  and  its 
follow  ups  [10,11]  have  been  cited  in  numerous  papers  on 
complexity  theory  and  related  algebraic  problems.  See  [12] 
for  an  (already  slightly  outdated)  introduction  to  this 
subject. 
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CHAPTER  1 


ANALYTIC  BEZOUT  IDENTITIES 


by 

* 

Carlos  A.  Berenstein  and  Alain  Yger 


Introduction 


1 . 

In  a  number  of  papers  12,3,4,6,7,12)  the  problem  of  finding 

explicit  solutions  for  the  Bezout  equation: 

f,h,*»-***+f  h  *  1  has  been  considered.  If  f,,***,f  are 
11  mm  l  m 

complex  polynomials  in  n  variables  and  they  have  no  common  zeros 

in  Cn,  the  existence  of  explicit  analytic  expressions  for  the 

corresponding  polynomials  hj,*‘*,h  has  a  number  of  applications 

to  control  theory  and  commutative  algebra.  The  most  notable 

application  being  that  of  Brownawell  [11]  where  the  essentially 

best  possible  estimate  for  the  degrees  of  the  h,  ,***,h  was 

l  m 

obtained  using  such  explicit  analytic  expressions.  Up  to  date  no 
purely  algebraic  proof  of  these  bounds  has  been  found. 

Similarly,  in  the  deconvolution  problem,  this  time  the 

functions  f,,**»,f  being  the  Fourier  transforms  of  a  strongly 
i  in 

coprime  family  •  •  ,fJm  of  distributions  of  compact  support, 

one  searches  for  a  procedure  to  compute  explicitly  distributions 

of  compact  support  such  that  /z,  *v .  +•  •  •  +/z  *v  =6. 

i  in  linn 

(Here  v  play  the  role  of  h^  In  the  equivalent  formulation 

v  el).  This  question  arises  in  problems  of  robust 

filtering,  image  processing,  etc.  [10).  In  [7]  we  wrote  down 

formulas  for  a  solution  v,,  ••*,*'  of  the  deconvolution  problem 

i  m 

in  terms  of  interpolation  series.  The  problem  we  have  faced 
recently  is  that,  while  for  the  1-d  case  these  formulas  can  be 
easily  Implemented,  in  the  higher  dimensional  case  they  are  far 

too  cumbersome,  some  of  them  seem  to  be  beyond  the  range  of 
symbolic  languages  like  MACSYMA  on  which  we  had,  perhaps  too 
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optimistically,  relied  upon  to  carry  through  the  computations 
Involved.  For  that  reason  we  present  here  a  new  version  of  our 
original  deconvolution  formulas  which  assumes  extra  conditions  on 
the  family  •••,/*  but  has  as  a  payoff  a  very  simple  formula 

for  the  deconvolutor6  •••,»>.  We  0*ve  herein  simple  examples 
where  these  extra  conditions  are  satisfied. 

The  problem  of  finding  an  efficient  algorithm  to  compute  the 
above  mentioned  solutions  to  the  algebraic  Bezout  equation  being 
still  open,  we  also  analyze  here  the  particular  case  in  which 
those  polynomials  can  be  computed  in  terms  of  interpolation 
formulas.  Finding  an  algorithm  with  a  low  complexity  for  this 
problem  will  have  many  important  applications  in  the  theory  of 
distributed  parameter  systems  and  in  robotics. 

We  have  also  found  that  a  language  barrier  prevented  our 
work  [7]  to  be  more  easily  available  to  some  engineers,  and  hope 
that  the  present  paper  will  overcome  those  shortcomings. 

The  first  author  would  like  to  express  his  gratitude  to  the 
University  de  Bordeaux  I  for  its  hospitality  while  this  work  was 
carried  on.  His  work  on  the  algebraic  aspects  of  this  paper  are 
inspired  on  the  questions  raised  by  the  AFOSR-URI  project  at  the 
University  of  Maryland.  We  will  also  like  to  thank  our  friends 
Dal€  Brownawell  and  B.  Alan  Taylor  for  many  helpful  remarks. 
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2 .  Analytic  case 

We  will  consider  only  entire  functions  f  of  n  complex 
variables  satisfying  inequalities  of  the  form 

m  H<Im 

(1)  If (z) I  *  A ( 1+ ] z J )  e  .  z  c  Cn  , 

Im  z  «  (1m  z,,***,Im  z  )  c  Rn  where  H  is  a  convex  continuous 
in 

function  in  Rn,  homogeneous  of  degree  1  (i.e.  H(Xx)  =  XH(x)  when 
X  >  0).  We  call  such  a  function  H  a  supporting  function.  By 
the  Paley-Wiener  theorem  [18]  there  is  a  distribution  p  of 
compact  support  in  Rn  such  that  f  *=  p,  the  Fourier  transform 
of  p.  Furthermore,  the  supporting  function  Hq  of  cv  supp  p 
will  satisfy  i  H  (Here  cv  denotes  the  convex  hull). 

Conversely,  if  f  *  p  we  can  take  H  *  Hq  in  (1)  and  m  is 
related  to  the  order  of  p  in  a  simple  manner.  Hereafter  we  will 
just  write  f  c  S'  (=  e'(Rn))  if  f  satisfies  (1). 


For  simplicity  denote  p(z) :=log(2+jz| )  +  | Im  zj.  A  family 

f,,’**,f  of  functions  in  S'  is  said  to  be  strongly  coprime  if 
l  m 

there  is  a  constant  c  such  that 

m 


(2) 


]T  1  f  J  (z)  I  2  *  e“cp(z)  ,  z  €  Cn. 
j  =  l 


It  is  well  known  [14]  that  (2)  is  a  necessary  and  sufficient 

condition  for  the  existence  of  functions  ' *  • ,h  e  t'  -such  that 

m 


(3) 


\ 

2L 


fjhj 


l . 


3-1 

In  other  words,  a  strongly  coprime  family  is  precisely  a  family  for 
which  the  analytic  Bezout  equation  (3)  has  a  solution.  If  we 
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consider  (3)  In  terms  of  the  distributions 

such  that  ■  f  ^  and  then  we  have  the  identity 

(4)  »l  •  +•••+  nm  •  «  6. 

i.e.  v  solve  the  deconvolution  problem  stated  in  17]. 

1  in 

We  will  8ay  sometimes  that  the  family  of  distributions  p  , 
is  strongly  coprime. 

It  might  be  useful  to  explain  why  is  (4)  called  a  deconvolu¬ 
tion  problem.  If  we  have  an  unknown  signal  (function  or  even 
distribution  or  random  process)  p  then  the  usual  data  one 
measures  would  be  given  by 

j.  m 

(3)  Pl  •  ^ffi  *  *>  • 


The  way  to  recover  *>  is  by  deconvolution  (which  is  still  given 
here  by  convolution  with  distributions  of  compact  support) . 


(6) 


P  m  V .  *  ?>, 

11  mm 


As  we  have  mentioned  in  the  Introduction  our  problem  is  to 
find  easily  computable  functions  g^  and  corresponding  distribu¬ 
tions  Vj  solving  (3)  and  (4)  respectively.  We  note  that  under 
the  strongly  coprime  condition  (2),  or  even  under  the  weaker 
assumption  that  (2)  is  only  satisfied  for  real  values  of  z  (zcRn), 
there  are  readily  available  tempered  distributions  solving 

the  deconvolution  problem,  namely  let 

f  Az) 


(7) 


Vz>  ■  nr 


n 


,  Z  e  R  . 
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The  problem  Is  that  the  o  do  not  have  compact  support  and 

furthermore,  the  themselves  sre  not  so  readily  computable 

(except  by  inverting  the  Pourier  transform) .  Nevertheless  there 
are  many  situations  where  these  are  still  very  useful,  among 

other  reasons  because  they  minimize  the  noise  amplification  of  the 
deconvolution  process  (6)  (see  (17)  for  an  example  of  implementa¬ 
tion  in  2-d) .  On  the  other  hand,  in  many  applications  it  is  often 
not  necessary  to  obtain  an  exact  solution  to  (4)  but  one  is 
allowed  to  replace  the  Dirac  6  in  (4)  by  a  (sufficiently)  smooth 
function  u  with  small  support,  i.e.  an  approximation  to  6 .  It 
is  this  approximate  deconvolution  problem  that  is  more  readily 
solvable,  even  with  very  good  knowledge  on  the  support  of  the 
distributions  Vy  which  will  turn  out  to  be  (reasonably)  smooth 
functions . 

Since  the  method  we  use  relies  on  Koppelman  type  formulas, 
like  those  developed  in  [1,9],  we  need  the  following  explicit 
relation  whose  proof  is  an  immediate  verification. 


Lemma  1 .  Let  p  be  a  distribution  of  compact  support  in  Rn, 

1  s  k  s  n  and  <;  *  (<  »•••,£ n)  e  Cn.  The  holomorphic  function  of 
2n  complex  variables  gk(z, K)  defined  by 


(8)  gk(z.t) 


*(z. 


k+1 


'cn>  - 


u  (z. 


'k-l 


•£n> 


is  the  Pourier  transform  (for  K  fixed)  of  the  distribution 

denoted  I  «  I(//,<,k),  which  evaluated  at  pc  C®(Rn)  has  the 

o 

value 
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(9)  <!,«>>: 


-  1  |  Q  •  .0)  e 


ic^u-v 


du 


“1(tk+lCk+1+ ’* ’ +tn*n)  . 

e  dM  t )  . 

(By  abuse  of  language  we  have  written  Jy(t)dp(t)  to  denote  <u,f>) 

We  note  that  for  the  distribution  I  we  have  cv  supp  I 
£  cv  supp  //.  Furthermore  the  collection  of  functions  g  ,***,g 
satisfies 


(10) 


g1(z,c) (Zj-C jH* * *+gn(z,o (zn-*:n>  ■  /*(z)  -MO 


Associated  to  these  holomorphic  functions  we  have  a 

differential  form  g  in  the  variable  <  given  by 

n 

(ID  g  •  g(z,  0/0:  =  ygv(z.(  )dtv. 

k=l 

Given  a  family  of  a  entire  holomorphic  functions 
zero  set  Z  is  defined  as 


(1,0) 


i  1 8 


(12)  Z:«  (z  c  Cn:  f  1  (z)- •  • -fm(z)  *  0  ). 

In  our  applications  we  will  only  consider  the  case  where  the  set 
Z  is  discrete.  We  say  that  2  is  almost  real  if  there  is 
constant  A  >  0  such  that 

(13)  Z  £  (z  e  cn:  Jim  z|  s  A  log(2+)zj)). 

It  is  well-known  that  an  almost  real  zero  set  Z  Is  discrete  (8], 
(15).  For  a  discrete  set  Z,  r  >  0,  we  can  define  a  counting 
function  n(Z  ,r)  :*#(ZnBr) ,  Br  ■  (zeCn:  |‘z(  <r }  *  Euclidean  ball  of 
center  0  and  radius  r.  The  distance  function  is  d (z,Z):« 


■ln(  1  ,ain{  I  z-t  I  :f  e7 } ) . 

«  •  I  * 

Given  a  family  of  n  distributions  of  coapact  support  in 

Pn,  let  us  denote  by  Hj  the  supporting  function  of 

n 

cv  U  supp  v  . .  that  is 
1  J 

(13)  H.(0):c  max  max(x0  :  xc  suppji.)  ,  (©ePn)  , 

is  j<n  J 

x-e  *  x.e , +•  •  •  +x  b  . 
ll  n  n 


Definition  1.  A  family  of  n  distributions 

compact  support  in  Pn  is  well-behaved  if  there  exists  positive 

constants  A,B,N,x  and  a  supporting  function  Hq  such  that 

0  i  H  s  H. .  such  that  the  zero  set  Z  of  the  functions 
o  l 

f,  f  *  .  is  almost  real, 

ll  n  n 

(14)  n(Z,r)  *  0(rA), 
and,  denoting 

^  1/2 

(15)  |  f  (z )  |  :  -  (2jVZ,j2) 

1 


the  following  inequality  holds*. 


(16) 


|f(z)  J 


i  8  d(z.Z)* 

(l+|z| ) 


Ho(Im  z) 

e 

N 


Definition  2.  A  well-behaved  family  is  very 

well-behaved  if  there  are  constants  Cj.M,  Cj>0,  such  that  for 
every  f  e  Z  we  have 

|J(0|:-  |d«t[^(C)]  |  *  Cjdtjt;  |  )~M  . 

1  1 ,  J 


<17) 


This  condition  implies  that  the  common  zeros  of  *i'**"'*n 
are  simple,  that  we  can  take  m  *  l  in  (16)  and  that  if  C,C'  c  2, 

.U* 

C  0  C'  then  *  c2(14|i;J)  for  some  positive  constants 

c2.M'. 

We  will  say  also  that  functions  f^,***,^  are  (very)  well- 
behaved  if  the  above  properties  hold. 

Given  a  family  f,,***,f  in  6'(Fn),  m  >  n,  with  no  common 

x  m 

zeros  we  introduce  the  following  functions  and  differential  forms. 
First,  let  g-*  *  g ^(z,t,Jfj),  fj  =  p  j ,  be  the  (1,0)  differential 

forms  in  C  given  by  (11),  we  write  g^  *  ^9kd^k’  Reca11  the 

k=l 

coefficients  g^  are  holomorphic  in  both  z  and  C .  Let  F  be 
the  vector  valued  holomorphic  function  F:  »  ( f j , • * • , fft) ,  we 
write  B 

(rtoi :-(  Y  ifjioi2]172  . 

3-1 

CD 

which  is  a  nowhere  vanishing  C  function  of  K •  Let 

m 

(18)  *  «  Viz.Z)  s  »  ]>  fj(C  )  f  j  (z)  / 1 F  (C  )  |2  , 

m 

(19)  Q  *  Q{z,U  :  «  ]>  f  j(t)g*(z,t)/|FU  )  |2  • 

Therefore  •  is  a  C*  function  of  *  1  and,  as  a 

function  of  z,  p  is  a  linear  combination  of  the  fj.  Q  is  a 
(1,0)  differential  form  in  K,  Its  coefficients  are  C®  in  (z,S) 

<SD 

and  holomorphic  in  z.  Finally,  the  n+1  functions  Aj,  C  in 
(z,<)  and  holomorphic  in  z  are  defined  by  the  identities 


1  S  j  s  n 


(20)  Q1A*  •  •  *AgnAQ  «  A  jdC  jA*  •  •  Adt;n  , 

(21)  o1*--V  .  An+J«,A---Adt:n  . 

It  is  clear  that  the  A^  are  6imply  n  x  n  determinants  whose 
entries  are  obtained  from  the  coefficients  of  g*,»**,gn,Q. 
Therefore,  as  functions  of  z,  they  are  finite  linear  combina¬ 
tions  of  products  of  n  among  the  functions  g^,  l<ksm.  Note 
that  these  products  are  just  Fourier  transforms  of  convolutions  of 
n  distributions  of  the  form  I(*/j,<,k).  (see  (9)). 

In  order  to  obtain  simple  and  easily  computable  deconvolution 

formulas  we  need  to  asstime  that  a  strongly  coprime  family  of 

distributions  //„,*••,//  contains  a  (very)  well-behaved  sub- 

l  m 

family  z/ .,•••,//  .  Furthermore,  we  need  some  control  on  the 
l  m 

relation  between  the  support  of  all  the  versus  the  supports  of 
the  first  n.  Let 

(22)  H 0(0)  *  max  max{x-©  :  x  e  supp  p.)  ,  (0eRn)  . 

l<jsn  ? 

One  6uch  relation  between  the  supporting  functions  is 

given  by 

(23)  H2  s  2Hr  and 

(24)  2 (n-1 )Hj (© )  +  H2(9)  +  H2(0)  <  2n  HQ(e)  if  0  x  0. 

The  last  condition  is  equivalent  to 

(25)  3rc  >  0  such  that  rQ|©|  s  2nHo(0)  -  2  (n-1  )H1  (© ) -H2  (©  )  . 

With  all  this  notation  in  place  we  are  now  ready  to  state  the 
first  deconvolution  formula. 
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Theorem  1  Let  //j. be  a  strongly  coprime  family  of  distri¬ 
butions  such  that  is  a  very  well-behaved  subfamily. 

Assume  further  that  (23)  and  (25)  hold.  For  any  u  c  C^R11)  with 

*  o 

supp  uc  (xc  Rn  ;  JxJ  s  ro)  one  can  write 

V'  *  An+l {z,<*  * 

(26)  u(z)  w(C)  — -  *>(z,0 

CeZ 


Proof  of  Theorem  j .  It  follows  the  lines  of  Theorem  3  from  our 
paper  (7).  It  uses  the  Koppelman  type  generalization  of  the 
Cauchy  Integral  representation  formula,  specially  in  the  version 
due  to  Berndtsson-Ander8son  (1).  One  Introduces  first  a 
parameter  c  >  0,  a  function  •  and  two  (1,0)  differential  forms 
In  K  as  follows: 
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n 


]Tf  jK)f  j(z)+c 


(28) 

♦>  U.O  :  «  - $ — 

c  { f (t )  |  2+r 

n 

(29) 

s(z,c)  :  *  y  K  i~zi)dC  i 

3  =  1 

n 

j(C  )gj(2,t;  ) 

(30)  Qr(z,C)  :  -  - 5 -  • 

£  |f(Ur+£ 

n 

where  as  before  f  *  (fj,***/^)*  Jf(C)j2  =  ^|fjK)|  •  The 

3=1 


procedure  from  ([9],  p.402  and  p.409)  gives  two  kernels  K£,P£ 

(i.e.  differential  forms  in  the  variable  K  of  type  (n,n-l)  and 
(n,n)  respectively)  such  that  if  v  is  a  holomorphic  function  in 
a  neighborhood  of  B^,  z  e  BR  then 


(31) 


v(z) 


1 

(2ni)n 


v(i:)K£(z,C) 

R 


v«)P£(z.U 


B. 


« 

-  • 

- 


These  two  kernels  are  defined  as  follows.  Let  G^(t) 
6j(t)  *  t,  we  denote  for  any  a  e 


tn  and 


(32) 


(33) 


o<a)  -  c'a|(z.n  -  gs  «,Jt^Ut0 


where  •  is  given  by  (18)  and  f£  by  (28).  With  Q£  defined  by 
(30)  and  Q  by  (19),  we  define 
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(34)  KC(Z,<) 


(a.)  (a.)  a  a.  a, 

3  G1  G2  °a(§Q£)  2a(5q)  2 


ao+<1l+<12*n"1 


ax  !a2 ! 


2(a0«.j 


(35) 


-I  V^TT  1 


(a,)  (a  )  _  a  a 

G1  G2  (dQc]  A  (*Q) 


Ctl+a2Bn 


where  a0'ai*a2  e  Everywhere  the  variable  z  is  considered  as  a 
parameter  and  the  d  derivative  is  taken  with  respect  to  t;  .  Due 
to  our  choice  of  function  G^,  the  index  a 2  can  only  take  the 
values  0  and  1.  For  this  reason  the  expression  for  P£  becomes 
particularly  simple 


(36) 


P£  *  *>(3Q£)n  +  n9>£{dQc)n~1  a  dQ. 


The  terms  (dQ£)n  1  and  (dQ£  )n  must  be  computed,  for  instance 


(37)  (dQ£)  *  3Q£A**-AdQ£  (n  times)  «  (2i)n  n!An+1 


JK  )c 


2 . . . n+1 


(|f(C)|  +c) 


where  dX  *  dX  (K )  -  Lebesgue  measure  in  C  .  (Me  have  eliminated 
the  variables  (z,<)  where  they  were  evident,  we  will  use  this 
convention  freely  in  the  rest  of  the  paper.) 


It  i6  clear  that  f£  and  Q£  are  singular  when  c  =  0 
precisely  at  the  points  K,  c  z  .  The  expression  (37)  shows  what  is 
the  strength  of  this  singularity  in  one  of  the  terms  of  P£  .  The 
strategy  of  the  proof  is  to  try  to  get  very  singular  terms  so  that 
when  c  — >  0  the  volume  integrals  in  (31)  become  sums,  while  the 
boundary  integrals  tend  to  zero  when  we  set  v  «  u  and  let  R  — »  ® 
over  a  conveniently  chosen  sequence.  The  reason  this  idea  works 
is  the  following  lemma  ([7],  Corollary  4.1.1): 
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Lemma  2  Let  o  be  the  measure  which  is  the  sum  of  Dirac  masses 
-  o 

od  n  \  ’ 

at  the  points  of  2,  i.e.  for  f  c  CQ(C  )  we  have  J*doo  *  >  *K) 

t;ez 


Then,  the  family  of  measures  o given  by 

£ 


(38) 


do€<C) 


(|m)|2+on41 


dl  u  ) 


converges,  when  c 
(39) 


0,  to  the  measure 


n 

nT 


n  do 


2  ' 


|j«>i 

where,  as  always,  J  denotes  the  determinant  Jacobian  of 


rl*  ' rn* 


From  (37)  we  see  that  the  first  term  in  (36)  is  amenable  to 
Lemma  2.  The  second  term  is  not  singular  enough,  therefore  it 
will  be  transformed  using  Stokes’  formula  in  the  corresponding 
integral  of  (31).  Namely,  due  to  type  considerations ,  one  obtains 
the  first  part  of  the  following  identity 

dt.{w[t.)<t>e(z,K){dQ£{z,K))ri'1  a  Q(z,U) 

=  d[v<pc(d Qc)n_1  a  Q)  *  v*c(aQc)n-1  a  a'Q  +  v  d<Pc A(aQ£)n_1  A  Q. 

The  last  identity  follows  from  the  fact  that  v  is  a  holomorphic 
function  in  K  and  the  (2n-2)  form  (aQ£)n  1  is  a-closed. 

Using  this  identity  the  representation  formula  (31)  becomes 

(40)  v(z)  -  — — -  f  v(r )  (K  ♦  n®  (5q  )n-1  a  Q) 

<*■*»"  JaBR 

♦  — - — -  f  v«)(*(dQ  )n  -  n  a®  a  (6>Q  )n_1  A  Q), 
(2ni)n  J  c  c  c 

®R 
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where  the  Integration  is  In  the  variable  <  and  we  have 
suppressed  the  dependency  on  [z.K)  of  the  kernels. 


Lemma  3  The  following  identity  holds 


(n+1  )n 


(41) 


d?cA(dQc )n_1  -  (n-1 )!(-!)  2 


n 


( Jf {<  ) |2+c )n+1 


^](-l)J(f i(z)-fi(K  ))  A  9k  a  A 
j  =  l  Vj*  1  k=l 


where  the  wedge  products  in  (41)  are  to  be  taken  in  their  natural 

order,  e.g.  /\  g  =  g  A«**Ag  . 
k>  l 

Proof  of  Lemma  3  We  start  by  rewriting 
n  n 

^fjUJfjU)  ♦  c  ^fjKMf^z)  -  f  (O)  +  |  f  ^  (H  )  |  2  ♦  c 
_ .  _ 


<t> 


£  1  f  («; )  1 


jfiur  +£ 


V'  ^j(^) 

-  1  ♦  >  ltAz)-tAO)  * - s —  . 

pi  1  3  |f«)|3« 

Denote  f^  *  fj(Z)  and  -  f^K)  :=  (  |  f  (<  )  |  2+c  )  f  ^  .  Then  we 

have 


n 


=  1  + 


Z«V»  '  W  V  Zr3®; 

j'l  J*1 


Therefore, 

a  {5Q£  )n-1  -  [X(fj(z)  "  fj)5vj)  A  A  flj) 


n-1 
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*  Z,f3u>  -  vK  a  ( lA a  5k]n'1] 

j*l  k*  J 


since  dfj  a  «  0.  Using  that  the  2-forms  3^  a  g  commute 

and  that  the  product  of  two  of  them  with  the  same  index  vanishes, 


we  have 


(Z  5rkA91 

k-J 


(n-1 ) !  /\  (dr  AgK) 
k*j  * 


(n-1 )!(-!) 


(n-l  )n 

2  (A  gk)A(A  ^v>k) . 

kr  j  k*  j 


Hence 


'j  a  f  Za%AB‘ 

k«J 


n(n+l )  n 

(n-l)'(-l)  2  (A^k>A<A^v)- 

k*j  k=  1 


Now.  we  have  =  (|f|  +C  )  dt ^  -  (|f|  +c  )  f  ^  d  |  f  j  .  Therefore 

—  2  —  2 

we  can  use  that  d J f |  a3 J f j  *  0  and  obtain 


n 

n  - 

&  A '  Tjti^o" 


If  we  now  expand  3  If  I  * 


j<ks 


(|f j2+£)n+1 


I'k 


we  see  that  only  the  term 


fjSTj  remains  in  the  triple  product  above.  Hence 

n 


n  , 

*k  ( )n 


j  =  l  *  _  r 

‘  '  7m7 


This  concludes  the  proof  of  Lemma  3. 


D 


Lemma  3  tells  us  that  Bfc a (§Qf )n_lAQ  Is  the  product  of  a 
measure  with  a  smooth  density,  independent  of  c  ,  and  the  function 
c  ( J f (< ) | 2+c )  n  1.  Lemma  2  can  now  be  invoked  to  see  that  the 
volume  Integral  in  (40)  reduces  to  a  sum  when  c  — ►  o.  In  fact, 
let  us  choose  R  so  that  |f(C)|2  *  | f (< ) | 2+* • • ♦ | f (< ) | 2  m  o  when 
|  *  R.  This  choice  is  always  possible  since  2  is  a  discrete 
set  by  assumption.  In  this  case  none  of  *>c  ,  ,  Q£  have 

singularities  when  c  «  0  and  C  e  3br.  We  set  ,K  Q  to  be 


o'  o,^o 


the  correspondient  quantities.  Therefore 


v(z>  ‘  TTTTn  1  v(,;nKo  +  '*0(<sVn‘1  A 

(2"1)  SB. 


i—  11b  f  v«)L(SQ,  Jn  -  n  S.  a(SQ,  l"'1  a  q! 
ni)n  £-*0+  „  L  C  C  C  J 


(2*1) 


Recall  that  fj(<)  ■  0  if  K  e  Z  and  1  s  j  <  n.  Using 
Lemma  2  and  the  definition  (20)  and  (21)  of  the  A  j  we  can 
compute  explicitly  the  limit  and  obtain 


(42) 


v(z) 


— -  f 

<2nl)n  SB, 


v(()(Ko  ♦  tvpo((aQo)n_1  A  Q) 


X  V(C,^1 


+y.-i.n+i-^(z)  y  v(t )* 3 ( z_— . 

J-l  (c/Tb„  J(I:  1 


Up  to  this  moment  we  have  only  used  that  2  n  3B_  ■  $  and 

K 

that  v  is  holomorphic  in  a  neighborhood  of  B_.  To  let  R  — ►  ® 
we  have  to  choose  a  sequence  R  — ♦  ®  judiciously.  Recall  that 
u(2,r)  «  #(2nBr)  s  CrA  for  some  positive  constants  A,C  and  all 
r  *  1.  Let  M  be  the  smallest  Integer  t  C(R+1)A  +  1,  divide  the 
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shell  ®R4j\br  into  H  concentric  subshells  by  choosing  the 
boundaries  to  be  9B(R+j/M)'  0  *  3  s  M-  There  is  at  least  one 
such  subshell  that  is  free  from  points  of  Z,  choose  R'  to  be 
the  mid-radius  of  this  subshell,  then  d [K  ,Z )  *  (2M)1  if 
|C  |  *=  R'.  Starting  from  the  sequence  R  =  q  *  1,2,*  ••  we 
construct  a  sequence  R^,  q  <  Rq  <  q  +  1 .  such  that  for  some 

positive  constants  A  , 


-N  H  (Im  K  ) 

(43)  |fK>|  *  A1  q  e  °  if  \<\  *=  Rq. 

This  follows  from  (16)  and  the  choice  of  Rq . 

We  are  now  ready  to  estimate  the  terms  in  the  boundary 

integral  of  (42)  for  R-Rq*  We  will  assume  |z|  s  Cq  <  od  and 

consider  those  q  such  that  R  *  C  +1. 

q  o 

First,  let  us  observe  that  the  functions  g£,  Is  j  s  m, 

1  s  k  s  n,  satisfy  an  estimate  of  the  form 

,  M.  M  H,(Ii  z)  +  H  (Im  r ) 

(44)  |g£(z,*:)|  *  Cjd+Jzj )  1(1+|<|)  e 


for  some  constants  Mi*c1  >  0.  If  1  s  j  s  n  we  can  replace  H2 
by  Hj .  We  can  now  estimate  the  coefficients  of  differential  from 
Q.  Denote  || Q ( z , t; )  |j  the  largest  absolute  value  of  the  co¬ 
efficients  of  d<k  at  the  point  (z,{).  We  proceed  as  follows. 
First, 

-M.  H  (Im  <  ) 

(45)  |F«)|  *  |f(t)|  *  AAq  e  °  if  |C  |  -  Rq . 


Therefore 


which  leads  to 


|| q ( z , < ) n  s  pT*yj-[  £||sJ(z,<)H2] 

3-1 


1/2 


Cjd+jzl)  2e 


M  H  (Im  z)  M,  H„(Itt  K.  ) 

“  (1+Rq)  e 


A1  q 


-n1  hou*  t; ) 


(46) 

||Q(z.<  )  H  s  C3q  *  e  *  ° 

The  constant  C3  depends  in  fact  on  z  ,  but  jzj  s 

Co 

and 

N2  - 

+  Nj  (In  fact  Cg  can  be  estimated  in  terms 

of 

e! Im  2I 

and  polynomials  in  | z ( . ) .  Similarly,  with  possibly 

different 

values 

for  the  constants  C3>K2  *PPear*ng  below  we  have 

(47) 

N,  H_ ( Im  K)  -  2  H  ( Im  K  ) 
JJaQ( z )  Jj  s  C3  q  Ze  2 

(48) 

N,  2H1(Im  <)  -  2  H  (Im  K  ) 
||aQo(z,«:  )  ||  S  C3  q  Ze  1 

(49) 

N.  -Hflm  K) 

|4>0(z,U|  *  C3  q  e 

(50) 

H  -H  (Im  C) 

*  {fItH  s  c3q  * 

To 

estimate  K  we  recall  that  a.  can  only  take 

O  2 

the 

values 

0  and  1  in  (34).  In  case  a  =  0,  we  have 

to 

estimate 

n-a  _  ax 

terms  of  the  form  f  )  ,  with  0  s  a,  s  n  -  1 .  There 

o  o  l 

are  powers  of  q  that  we  will  disregard,  the  estimate  is  then  the 
following  functions  of  Im  K 
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n-a  0  -(n+n  +  1)H  *2o  H 

»*„  1*<a<30>  'll  «  «  1  °  ‘  • 

Since  H.  *  H  #  the  worst  case  estimate  occurs  when  a.  «  n  -  1. 
1  o  i 

Hence  the  terms  corresponding  to  a2  ■  0.  aQ  +  cij  ■  n-1,  in  the 

definition  of  K  can  all  be  estimated  by 

o 


N,  -2nH  (In  o  ♦  2(n-l)H.(Itt  t) 

(51)  C3  q  3e  °  .  M  -  R,  • 

The  terms  with  a  ■  lf  ao  +  ai  “  n  “  2  correspond  to  the 

n-a .  _  a. 

estimate  of  |jf>Q  1  (dQ  )  AdQjj.  The  worst  case  occurs  this  time 
when  cij  *=  n  -  2  and  we  obtain  an  estimate  of  the  form 


(52) 


c3  q  e 


N3  -2nHo(Im  <)+2(n-2)H1(Im  C)+2H2(Im  O 


■1*1  -  V 


In  (42)  we  have  one  more  term  to  estimate  for  |<  |  *=  Rq, 


(53)  |k0(aQ0)n_1  a  Q ||  s  C3  q 

The  conditions  (23)  and  (25)  imply  that  the  largest  exponential 
factor  in  (51),  (52)  and  (53)  i9  the  one  in  (53)  and  it  satisfies 


N3  -2nHo(Im  C)+2(n-l)H1(Im  t)+H2(Im  K) 


(54) 


-2nH0(Im  K  )  +  2(n-l)H1(Im  K  )  +  H2(Im  Us-  rjlm  * 


<JD  .  - 


since  we  have  assumed  that  u  e  C  (B_  )  we  have 


o '  r 


-N  -2n  r  | Im  C | 

(55)  |u(U|  s  C4(l+|i:  (  )  J  e  °  ,  t  c  Cn 

which  allows  us  to  conclude  that,  if  v  ■  u  , 

(56)  11m  f  u(u[k  (z,U  +  (z,U(3Qo(2,U)n'1  a  Q(z,C)1 

q-tQD  1  L  -* 


0. 


db. 
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To  conclude  the  proof  of  Theorem  1  we  need  to  show  that  the 
series  appearing  In  the  representation  formula  (26)  converge 
absolutely  and  uniformly  In  compact  subsets  of  Cn  to  functions 
In  €'.  Since  we  have  assumed  that  Z  is  almost  real,  the 
estimates  of  all  the  terms  Aj.e.J-1  are  in  terms  of  powers  of 
\K  | .  Using  that  n{Z,r)  «  0(rA)  and  |u(<)|  decreases  as  fast 
"N4 

as  J <  J  for  t,  e  Z  we  have  the  desired  convergence  once 

is  chosen  sufficiently  large.  The  support  of  all  the  distribu¬ 
tions  thus  obtained  is  contained  in  the  convex  set  K  whose 
supporting  function  is  (n+l)H2-  o 

Remarks  1.  One  can  see  that  the  condition  211^  *  H2  cannot  be 
relaxed  if  the  other  conditions  of  theorem  remain  in  the  same, 
otherwise  the  exponent  in  (52)  would  become  positive  and  we  would 
not  be  able  to  prove  (56). 

2.  A  way  to  weaken  the  conditions  in  Theorem  1  is  to 
impose  6ome  better  lower  bound  on  { F |  than  (45)  that  depends  only 
on  the  first  n  functions.  We  will  do  so  in  Theorem  2  below 

N  - 

3.  It  is  clear  that  one  needs  only  u  e  C  (B  )  for  N 

°  ro 

sufficiently  large  to  obtain  (26). 

The  following  example  shows  how  Theorem  1  simplifies 
enormously  the  computation  of  the  deconvolution  formula  proposed 
In  (7J. 

Let  «  ,n  be  the  characteristic  functions  of  the 
r  1  r  2  r  3 

squares  centered  at  0,  of  sides  parallel  to  the  axes  and  of  length 
2/3  ,  2/2  ,  2  respectively.  One  can  easily  show  (7],  (17)  that 
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Is  •  very  well-behaved  family  with 

H0(lm  C  }  -  /2(|lm  Cj  ♦  j  Im  K  2\ ) 


Here  H^Im  O  *  HjIIb  t)  *  /3{|lm  1. 1 1  ♦  |  Im  C2|).  In  this  case  the 
main  hypotheses  (24)  reduces  to  verify  that  4/2  -  3/3  >  0  !. 

Since  Jxjj  +  Jx2|  *  /Xj+Xj  for  x  e  R  ,  we  have 

rQ  =  4/2  -  3/3  *  0.2.  The  variety  Z  is  this  case  is  given  by 


z 


j.k  e  Z 


>  {£•§)' 


j. 


k  e  Z 


(There  were  about  forty  different  types  of  terms  in  (7) 
compute . ) . 


Before  we  proceed  to  state  Theorem  2  we  need  to  point  out 

that  the  representation  formula  (31)  does  not  depend  on  the 

particular  choice  of  the  differential  forms  g^  we  have  chosen, 

rather  on  the  fact  that  (10)  is  satisfied.  That  is, 

n 

]Tgi[(z,C){zk-i:k)  -  fj(z)  -  fjK). 

k=i 

Now,  let  f  *  and  h(z)  =  sinBZj  for  some  B  >  0,  and  denote  by 
g  the  differential  form  associated  to  f  by  (8)  and  (11).  Let 
us  define  a  differential  form  r  by 

sin  Bz.  -sin  Bt; 

(57)  z(z,t):«=  f(t)  - - -  dtj  +  h(z)g(z,t), 

writing  f  «  we  1,8 ve 

r1lz,K)  (Zj-C  j)+*  •  •+zn(z,C  )  (zn-Cn)  -  f(z)h(z)  -  f(()h(t) 

therefore  we  can  associate  r  to  the  product  f*h.  It  is  also 
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clear  that  as  a  function  of  z  the  are  Fourier  transforms  of 

distribution  of  compact  support  easily  computable  terms  of  v  and 

B.  Obviously  we  can  replace  slnBZj  by  sinBz^  without  any 

problems,  hence,  given  a  family  f,,***,f  we  can  construct  an 

augmented  family  t ^ , • • • , fB, fB+1 :•* j *  sin  BZj  ,  •  •  •  , 

f-  :=f  •  sin  Bz,  ,  •  •  •  ,  f  .  .  :=f  ■  sin  Bz„  .  The  corresponding  g^ 

m  i  i  n+ 1  j  m  m  n 

for  j  2  m+1  are  computed  following  the  procedure  (57).  It  is 

clear  that  if  f,,***,f  was  strongly  coprime,  the  augmented 

l  m 

family  remains  strongly  coprime.  If  form  a  very 

well-behaved  family  we  will  keep  the  notation  Ho,HlfH2  to 
indicate  the  support  functions  corresponding  to  the  m  original 
members  of  the  augmented  family  fJ , • • • , f (n+1 )a  • 


Theorem  2 

Let  fx. 

•  •  •  f 

'  rm 

be  a  strongly  coprime 

family 

such  that 

the  subfamily  f^,* 

•  •  f 

*  rn 

is  very  well-behaved. 

There 

are 

constants 

Bo  *  0, 

r  >  0 
o 

such  that  for  any  B 

1  Bo' 

and  any 

u  e  CQ(Br  )  the  representation  formula  (26)  is  valid  for  the 
o 

augmented  family  f ^ , • • • , f^, • • • , f (n+1 defined  above  if  either  of 
the  following  two  conditions  holds: 


(58) 

H2  s  2Hx 

and 

2{n- 

1)H1  <  (2n-l)HQ 

(59) 

2HX  <  H2 

and 

2  ( n- 

2  JHj  +  H2  <  ( 2n 

Proof  The  proof  is  exactly  the  same  as  that  of  Theorem  1  except 
for  improvements  on  the  estimates  (46), (47)  and  (50)  for  the  new 
Q  and  p.  Recall  that  it  is  there  where  all  the  functions 
fl'*"'f<n+l)m  *PP®ar-  Let  pi5*(fi'*** 'f(n+l)m)  and  keeP  the 


notation  F  ■  (f 2  .• • •  .f,)  •*  before.  We  have 

(60)1?^)!,  «  |(1+  |®in  4-|sinBTfcl2)1/2  ■  l  ^ 

It  is  clear  that  for  some  positive  constants  cn*cn  we  have 


d(C)  *  (l+JsinBCjl  +••  •♦|sin  BCn|  )' 

'  2B |  In  t,  j  2B|Im  (  |' 

*  c'  e  ♦•••+© 

n 

* 

2B | Im  C| 1  2B j Im  < { 


2.1/2 


*  c  e 
n 


*  c  e 
n 


n  n 

where  |  Im  t,  |  ^  ^  1 1®  K  ^  |  *  I Im  ^  I  =  I Im  *•  j  I  ] 


1/2 


We  estimate  first  g^,  j  *  m+1,  for  ) z-K J  *  1  since  that  is  the 
only  case  that  appears  in  the  proof  of  (56).  As  it  follows  from 
(57)  we  have,  for  some  i,k  (lsisn,  lsksm)  ,  the  estimate 


B I  Im  z .  j  r  9i/2  lc  1 

||  g 3  (z,K  )||  s  Ce  1  (j  f^tt  )  |  ( 1  +  |  sin  BC  A  |  )  '  +  |fg  (z,C)||J 


N.  N.  B | Im  Z.|  ' 

s  C  ( 1+  J  z  |  )  1(1+|C  |  )  1  e  |f(U|^(C)  ♦ 

•» 


H2(lm  z)  H2(Im  K) 
e  e 


It  follows  that 


HQU.Oll  s 


(n+1  )m  1/2 

X  llg3<*.Oll2 

J=i  J 


*  ( i+ 1 z j )  1  < i+K  | ) 


Hx  B  j Im  z|+H2(Im  z] 


H2(Im  K) 


FK  )  |»K  ) 

I'jii)! 
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Therefore,  for  a  positive  constant  C2  depending  on  z,  we  obtain 

W  r  H(I®  ^)-H  (In  5  ) 

(61)  ||Q(z,C ) ||  s  C2(l+jC|)  *  e  2  °  h  M. 


Similary , 
(62) 

Pinally , 


B, 


|F  (z)|  N  H  (I®  K)  -S|I»  K  J. 

i^u.ui  *  yrjen  5  1  }  e 


|j3f,U>||  fl  (2,C  )|J 

{|3q(z,C)|(  s  C  max  - * - - — - 

lO-b(n+l)»  |Pj  (Ur 


Now,  for  j  *  m+1  we  have  for  some  is  1  s  n,  l  s  *  s  ® 


llBf  )H  *  IjsinlKj-af^O  +  f,(t;  )B-cos  Bi^dCjj 


N 


?  H5(lm<) 


*  Cjll  +  jtl)  (C  )  e  4  J  PU)}»(0>. 


It  follows  that 


**2  r 


iisQ(*,t)ii*  c  i^.im 

iFjKir 


S  C2u+I<|) 


N. 


H 


,  „  .H.dmU  2H  ( Im  T  )] 

|F(C)|  •«)  ♦2|PU)|d(U S  *  +«U)e 


B 


l  +  2e 


2  ( Im  C  )  -Hq  ( Im  C  )  2H2<  In  0“2H0(I®  ?  )-j|  Im  C  ^ 

^  0 


Choose  B  *  0  such  that 
o 


(63)  2H2(Im  K  )  -  2H0(Im  K)  -  ^  | Im  C  \1  s  0. 

When  B  *  Bq  we  will  have  j|Q{z,C)||  *  C2^1+KI)  1  and 

N  H  {Im  r ) ~H  (Im  *  ) 

(64)  !!*Q(z,r)jj  s  c2(i  +  jrj)  e  ° 

where  C2  still  denotes  a  constant  depending  on  z  of  the  form 
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H  H  CI»  *)*B|Xm  *| 

Cj  *  const. (14 | x | )  *  •  *  . 

Me  csn  now  return  to  the  proof  of  Theorem  1  st  the  point 
where  we  obtained  the  estimates  (51) , (52) , (53) .  Ignoring  powers 
of  q  the  exponential  factors  are: 

(SI* )  exp(-2nHo(Im  C)  ♦  2  (n-1  )R1  ( Im  K  ) -jj  |  Im  K  |  j  ) 

(52)  exp(-(2n-l)H0(Im  O  ♦  2  (n-2  )Hj  ( Im  t,  )+h2  ( Im  K) ) 

(53)  exp(-(2n-l)H0(Im  U  +  2(n-2)H1(Im  K) ) . 

Under  the  hypothesis  (58)  the  largest  of  these  three  is  (53') 
and  Its  exponent  satisfies 

(53")  -(2n-l)Ho(Im  <)  +  2(n-l)H1(Im  t)  s  -  rjlm  K\, 

for  some  tq  >  0.  If  the  hypothesis  (59)  holds,  then  the  largest 
exponent  Is  (52 ')  and  we  define  rQ  >  0  by 

(52")  -  ( 2n-l )  Hq(  Im  <)  +  2  (n-2 )  Hj  ( Im  <  )  +  H^Xa  <)  s  -  r0|lm  <  |  . 

In  either  case  the  rest  of  the  proof  is  the  same  as  that  of 
Theorem  1 .  o 

Example  As  shown  in  [7]  the  family  Pj,p2,p3  obtained  by  taking 
p  *=  characteristic  function  of  the  unit  square  ,p 

a  rotation  of  Pj  by  36* ,  and  p^  a  rotation  of  Pj  by  45*  , 

satisfies  the  first  conditions  of  Theorem  1  and  Theorem  2  with 
Hq(6)  ■  |$|  since  the  squares  contain  the  unit  disk.  One  can 
easily  convince  oneself  that  the  hypothesis  (24)  does  not  hold 
(e.g.  take  f«(t,0),t>0.)  On  the  other  hand  one  can  take 
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H  (9 )  «  h  (0 )  »  /2|0|  since  mil  the  squares  sre  contained  In  the 
disk  of  radius  /?  .  He  are  In  the  situation  of  hypothesis  (58) 
and  Its  verification  reduces  to  the  fact  that 

r  «  3-2/2  >  0 
o 

Furthermore  B  -  4(/2-l)  works  In  this  case, 
o 
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3.  Polynomial  case 

The  conditions  on  Theorem  1  and  2  imply  that  the  convex  set 
defined  by  Hq  contains  a  ball.  If  we  want  to  prove  an  algebraic 
version  of  (26),  the  fact  that  this  condition  i6  not  satisfied 
plays  a  role.  Such  a  representation  was  stated  in  [3,4]  without 
proof.  We  analyze  here  the  conditions  under  which  it  Is  valid. 


Theorem  3 .  Let  p  ,***,pB  be  a  family  of  polynomials  in  Cn 
without  common  zeros,  suppose  further  that: 


a)  D:  =  max  deg  p.  *  deg  p.  for  1  s  i  s  n. 

b)  Z  *  (z  e  cn:  p^z)  ®  0;  1  s  i  s  n)  is  discrete. 

c)  J ( z ) :  =  Jacobian  determinant  of  p1#**»,p  at  z  is  n  0 
for  all  z  c  Z 

d)  P1.*,*.Pn  have  no  zeros  at  infinity,  i.e.  #Z  *  Dn. 


Then 


n 


r — >  d  .  -  (  Z  #  C  )  r — »  n+1-H  V — '  4  .(Z,C) 

<”>  1  -  2.  -TTTi — *u.t)  +  V*’2-  -JTT!— 

<cz  3. a  (ez 


where  are  defined  as  in  (18)  -  (21)  with  respect  to  the 

polynomials  p,,***,p  . 

l  m 


Remarks  1.  The  functions  defined  by  (8)  are  obviously 

polynomials  of  degree  D-l.  It  follows  that  (66)  has  the. form 

(66)  p  ( z )  A  (z)+-“*+P  ( z )  A  (z)  ■  1 

11  mm 

for  some  polynomials  Aj  s  C[z1#***,zn]  of  degree  at  most  n(D-l). 
This  follows  from  the  fact  they  are  given  as  n  *  n  determinants 
involving  the  g^  . 
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2.  As  before,  the  esse  *  »  n  ♦  1  leads  to  a  particularly 


pleasing  fora  of  (66) 


(67) 


i «  y  _i _ i _ 

A  ju>  P^U 


B*U.U 


Pj  (2). 


n+1 .  „  % 

»n  «*>'» 

pn*l(*> 


3.  The  two  statements  in  condition  (d)  above  are  really  a 
form  of  Bezout's  theorem  (19).  The  meaning  of  the  expression 
,,P1,***,Pn  have  no  common  zeros  at  infinity"  is  that  if  we 

D  f2l 

introduce  homogeneous  polynomials  H.(z  .  ***,z  ):  =  z  p.  — 

j  o  n  o  j  z  7 

o 

then  the  subset  of  Cn+1  defined  by  {z  =0,  H=* • • =H  =0)  is 

o  l  n 

(0).  This  is  equivalent  to  the  statement: 


(68)  |  P  ( 2 )  |  -  [X,PJ(Z)|2) 

J«1 


1/2  *  C|z|D  if  | z |  *  R  si. 


If  we  call  Pj  the  leading  homogeneous  polynomial  of  Pj 
then  the  statement  (68)  is  also  equivalent  to: 


(69) 


(z  c  Cn:  p°(z) 


=Pn(z)  *  0)  «  (0) 


We  also  note  that  (69)  implies  (b)  above.  That  is,  condition  (d) 
above  implies  condition  (b). 


Proof  of  Theorem  3.  The  proof  of  the  same  as  that  of  Theorem  1. 
This  time  we  take  v*l  and  R  arbitrary  *  Rq  (cf.(68)).  One  can 
estimate  ||  Q  <  z ,  <  )  ||  *  C/\K\.  ||5Q(z,t)||  s  C/|C|2.  |f>0(2,U|  *  C/|<|D, 
||dQ  (z,C  )||  *  C/ 1<  | 2  if  |<  |  -  R  and  |z|  s  K  s  R-l,  with 
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C  *  C(k)  >  0.  These  utlaitn  Imply  that  the  boundary  Integral  in 
(40)  tends  to  zero  when  F  — *  ® .  o 

o 

If  Pi****»Pn  mre  *^ch  that  their  leading  terms  satisfy 

(69)  but  their  degrees  deg  p^  *  Dj  are  not  all  equal  or 
max{D^:l  *  J  i  n)  is  smaller  than  D-max  (D^lsjsm),  then  we  can 
still  prove  a  version  of  Theorem  3.  That  corresponds  to  the 
analytic  counterpart  of  Theorem  1,  that  is  to  Theorem  2.  For  the 
moment  we  continue  to  assume  that  J(z)  *»  0  VzeZ*  (p1«*  •  •  cPn*0) . 

Let  Ljfz)  «=  ui2i  +  ‘ *  * +unzn  a  **near  homogeneous 

polynomial  with  generic  coefficients.  The  condition  that 

(z  e  Cn :  L1=p°*=***P2  *  0)  *  (0)  is  an  algebraic  condition  on  the 

coefficients  of  .  Therefore  we  can  choose  L ^  such  that  for 

d 

any  integer  *  0  we  have  (p^L^  *  P2*:***'=Pn  =  °)  *  (°)  • 
Continuing  in  this  fashion  we  can  choose  Lj,«**,Ln,  dir*"'dn 
such  that  for  any  choice  of  constants  c^•'”•cn•  we  define 

Pj  **  Pj  ^ )  .  then  P°*"**-P°  «  0)  *  (0)  and  deg  p  ^  =  D  for 

1  i  j  s  n.  It  is  clear  now  that  for  most  choices  of  c ^  we  still 

have  that  all  common  zeros  of  p°,#,,,p°  are  6lmPle  and 

p?, • • • ,p°,p^. ,  ,  •  •  •  ,p_  have  no  common  zeros.  Theorem  3  can  now  be 
i  n  n+ i  & 

applied  to  this  new  family,  one  obtains  polynomials  Aj 

m 

AjPj  *  1,  deg  Aj  s  n(D-l) 

J-l 

and  such  that  they  have  a  representation  of  the  type  (65). 
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Me  remark  that  a  representation  such  as  (65)  cannot  be  valid 
If  the  p  ,‘««#pn  have  coaaon  zeros  at  Infinity.  For  Instance. 
In  the  example  of  Masser-Phlllppon  [11]: 


P1  “  zl'p2 

one  knows  that 
the  degrees  of 

polynomials 


P«- 


n-1 


P 


1-«B  izn'2 


=  Dn-Dn  1  Is  the  best  estimate  possible  for 
solving  the  polynomial  Bezout  equation.  The 


.n-2 


A  *  -  z 
n-1  n 


K-2  'Cl  1 


n-1 


Zn-2~Zn-l 


V 


_n- 1  , 
,  D  6 

1-z 

n-1  n 


have  exactly  this 


degree.  On  the  other  hand  if  we  had  a  representation  like  (65)  we 
could  conclude  that  there  are  solutions  A  of  the  polynomial 
Bezout  equation  with  deg  A.  s  n(D-l)  like  in  (70). 


We  would  like  now  to  show  that  the  condition  (c)  of  the 
simplicity  of  the  zeros  in  Theorem  3  is  not  necessary.  Regret¬ 
fully,  we  only  know  how  to  do  this  in  the  case  where  m  «=  n  +  1. 


Theorem  4 .  Let  p  .•••,pn+1  be  a  family  of  polynomials  in  C 

without  any  common  zeros,  D  *=  deg  p^'^deg  pn  i  deg  Pn+1  and 

{p°=...=p°  s  o)  =  (0).  Then  we  can  find  polynomials  A ^  of 

•n+1 

degree  s  n(D-l)  satisfying  the  identity  ^  AjpJ  =  1 '  T^e 

1 

coefficients  of  the  A ^  can  be  written  in  terms  of  the  values  of 

p  and  values  of  derivatives  of  p  and  the  coefficients  of 

n+1  n+1 

the  g£+1(z,t;)  (when  considered  as  polynomials  in  z),  all  of 
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these  evaluated  at  the  points  of  ?  ■  (f  :p  **mP  (^ )  “  0)* 

1  n 


Proof  Por  0  <  «  1  the  function  l/pn+J  1*  holomorphlc  In  a 

neighborhood  of  D  ■  (IPjl  *  cj»***»|Pn|  1  #et  *•  • 

compact  polynomial ly  convex  set.  (The  compaclty  follows  from  the 
condition  (66)).  By  Sard's  theorem  one  can  choose  the  c ^  so 
that  the  sets  (|Pj|s£j)  ®re  real  analytic  submanifolds  of  Cn. 
(In  fact,  we  only  need  It  in  a  nelgborhood  of  D.)  For  any 
v  c  X(ff)  we  have  that  Integral 


(71)  Res?  (vdC  jA*  •  •  AdCn) 


( 2n  1 ) 


« \n  J 


v(C ) 


dt  2a< 


Adt; 


n 


Pj(<  )*  *  *Pn(^  ) 


ipii“CTi 

see 

Ip  I  =o 

1  *n'  n 

is  independent  of  choice  of  ai'***'CTn  as  l°nS  as  °  <  aj  s  cj 

and  the  (|PjJ  *  <7^)  are  smooth.  Furthermore,  if  v  is  in  the 

ideal  generated  by  p1,*«*,pn  *n  *(D)  then  this  residue  is 
zero.  Therefore,  it  depends  only  the  values  of  v  at  Z  and  a 

certain  number  of  derivatives  of  v  at  2  (as  it  follows  from 

the  Nullstellensatz  as  presented  e.g.  in  [13J,  (16J).  In  other 
words,  the  integral  (71)  can  be  considered  as  an  operator  defined 
by  a  certain  linear  combination  of  the  Dirac  masses  b^  and  their 


derivatives 


la! 


3C 


b  ,  <  e  Z,  applied  to  the  holomorphic  function 


v.  This  operator  is  very  hard  to  compute  explicitly  except  in 
very  simple  cases  but  it  is  perfectly  defined  as  the  common  value 
of  all  the  Integrals  (71).  It  is  called  the  residue  curren/  of 
Z  . 
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Let  us  consider  now  the  polynomial  B(z)  defined  by: 


(72)  B(z)  -  Res 


gj(z,t )  •••  gJ4l(z,C  ) 


fl£(z,0  •••  g"+1(z.C) 


PjU) 


•••  *W2) 


n+1 


dC  ,a*  •  •  AdC  _ 

1  n 


This  polynomial  16  in  fact  of  the  fora  ^A^(z)Pj(z),  with 

3=1 

polynomials  of  degree  s  n(D-l)  whose  coefficients  are  computed 
in  terms  of  the  values  of  derivatives  of  Pn+1  and  the 
coefficients  of  g^  (as  polynomials  in  z)  evaluated  over  Z. 

The  only  problem  is  to  show  B(z)  si.  We  fix  values  0  y0<°  j<c  y 
Consider  complex  numbers  ai'**’'an  sufficiently  small  and  so 
chosen  that: 

1)  All  the  common  zeros  of  p,-o, ,»»-#p  -a  are  simple  and 

i  l  n  n 

lie  in  <|Pl|  <  §<V’*,'|Pnl  <  K1’ 
ii)  Pj-Oj  ,  •  •  *  'Pjj-^n'  ^n+1  have  no  common  zeros. 

Note  that  Pj-ki^  *  ’  *  *  * Pn^n  still  do  not  have  any  common 

zeros  at  od.  Let  us  denote  Z  =  (z  c  Cn:  p„-a  =•  •  •  «=p  -a  *0). 

o  li  n  n 


Then 


(73)  Res_  ( vdC -a*  •  •  AdC  )  =  - -  [  v(C  ) 


dC  y  •  •  AdC  n 


(2"ir|p*  |-o 
•  1  •  1 

Ip  I "cr 

1  *n'  n 


where  J(0 


^(Pl"ai,# *  * ^(P,»***.P«) 


n  n 


n 


8  (t  j  ,  •  •  •  it  n) 


and  the  last 


Identity  follows  from  Stokes*  theorem.  (Replace  the  contour  by 
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little  sphere*  about  the  distinct  points  of  T  ).  This  lest 

o 

Identity  Is  the  effective  computation  of  the  residue  current  of 
2a,  namely 

Res?  -  V 

The  first  Identity  In  (73)  shows  that  Res?  — *  Res2  as  o  — »  0, 

a 

i.e.  the  residue  currents  at  0*0  are  continuous  when  acting  on 
holomorphic  (n,0)  forms.  On  the  other  hand,  by  Theorem  3  we  have 

gJ(z#C)  •••  g"+1(z,<) 

•  • 

•  • 

•  • 

siu.o  •••  o"+1u.<> 
pl(z,-‘V  pn-H(z) 

(Note  that  the  g^  corresponding  to  p^  and  to  p^  - 
coincide).  By  continuity  we  obtain  B  *  1.  This  concludes  the 
proof  of  Theorem  4.  □ 

Remarks  l.  The  can  obviously  obtain  the  saune  result  without 
assuming  the  degrees  of  Pj,***,p  coincide  or  that  they  are 
larger  or  equal  than  that  of  Pft+1 • 

2.  The  reasoning  of  Theorem  4  extends  to  a  strongly 
coprime  family  of  n+1  elements  whose  first  n  members  from  a 
well-behaved  family.  Under  the  other  conditions  of  Theorem  1  or 
Theorem  2,  we  obtain  a  series  representation  of  the  solutions  of 

the  Bezout  equation  which  we  computed  in  terms  of  the  residue 
current  associated  to  Z.  This  time  the  series  converges  after 
grouping  of  terms. 
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3.  The  Interest  of  the  theorens  In  this  section  lies  In 

the  sesrch  for  explicit  algorithms  to  obtain  eolutlons  A. 

a  3 

for  the  algebraic  Bezout  equation  «  1  which  satisfy 

1 

u 

Brownawell's  estimate,  deg  A^  <  3  //nD  ,  //  «  mln(n,a). 

4 .  Conclusion 

We  have  shown  how  explicit  solutions  to  the  analytic  and 

algebraic  Bezout  equations  can  be  obtained  under  natural 

restrictions  on  the  original  functions  f,,***,f  .  This  work  has 

l  m 

applications  to  the  implementation  of  deconvolution  for 
multidetector  systems. 
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CHAPTER  2 


EXACT  DECONVOLUTION  FOR  MULTIPLE  CONVOLUTION  OPERATORS  - 
AN  OVERVIEW  PLUS  PERFORMANCE  CHARACTERIZATIONS  FOR  IMAGING  SENSORS* 

by 

t  tt 

C.A.  Berenstein  and  E.V.  Patrick 


To  Appear  In  IEEE  Procurement  in  Multidimensional  Signal  Processing 


1.  INTRODUCTION 


Throughout  the  last  several  years  Mathematical  results  have  been  present¬ 
ed  which  form  the  foundations  for  the  use  of  multiple  (parallel)  linear  opera¬ 
tors,  each  given  by  convolution  with  a  distinct  kernel  (or  impulse  response), 
in  place  of  the  use  of  a  single  such  linear  operator  or,  equivalently,  in 
place  of  the  use  of  multiple  (parallel)  operators  each  with  the  identical 
kernel  (1]  -  [13].  See  Fig.  1.  The  mathematical  results  cited  above  describe 

the  conditions  under  which  compactly  supported  distributions  p, ,  p_ .  p 

12  m 

have  associated  to  them  compactly  supported  distributions  v, ,  i>_ .  v 

12  m 

such  that 


m 


Iv*i 

i  =  l 


(1) 


where  6  is  the  Dirac  distribution  on  Rn  and  where  •  denotes  convolution. 
We  often  refer  to  the  p^  as  convolvers  and  to  the  as  deconvolvers. 

To  introduce  some  explicit  multiple  operators  and  their  role  in  deconvo¬ 
lution,  let  us  look  at  an  example  and  the  results  of  a  computer  simulation 
[14].  Our  example  is  for  a  case  in  which  equation  (1)  holds  for  m  =  2,  and 
it  is  outlined  in  Fig.  2. 

In  the  center  of  Fig.  2  the  block  diagram  of  Fig.  1(a)  is  reproduced  for 
the  case  m  =  2.  (The  operator  $  in  Fig.  1(a)  will  be  discussed  later.)  We 
shall  be  considering  in  this  example  the  one-dimensional  case,  that  is,  func¬ 
tions  of  one  variable.  To  the  left  of  the  block  diagram  is  the  graph  of  the 
input  signal,  that  is,  the  function  f,  which  here  consists  of  the  sum  of 
translated  Gaussian  pulses.  (This  function  is  one  that  is  frequently  used  for 
the  evaluation  of  deconvolution  algorithms;  see  for  example  [15].) 

The  input  signal  f  is  acted  on  by  two  convolution  operators,  one  with 
kernel  p}  and  one  with  kernel  p^.  The  resulting  output  functions  are  g^ 
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and  (£,  respectively: 

gj(x)  *  (f-PjHx)  *  f(y)Mj(x-y)dy,  1  *  1,2,  x  c  R. 

**R 

In  this  example,  the  kernels  Pj  and  p2  are  essentially  the  simplest  pos¬ 
sible  kernels  with  compact  support:  the  kernel  functions  are  each  constant 
over  some  bounded  Intervals  and  zero  elsewhere.  (The  support  of  the  kernel 
function  refers  to  (the  closure  of)  the  set  on  which  the  function  is  non-zero; 
the  support  Is  called  compact  if  this  set  is  bounded.)  The  kernel  functions 
are  shown  in  Fig.  2:  is  supported  by  the  interval  1-1,1),  and  p2  by 

the  interval  l~V2,y/2]. 

The  output  functions  and  g2  are  also  shown  in  Fig.  2.  We  often 

refer  to  these  as  the  data  functions,  for  in  applications  gj  and  g^  are 
typically  the  results  from  two  distinct  measurement  operations  which  Eire 
modeled  by  the  "convolvers"  p1  and  p2-  An  application  that  is  a  physical 
realization  of  the  situation  presented  in  Fig.  2  is  spectroscopy.  In  that 
case,  the  function  f  is  an  unknown  density  function  with  the  wavelength  of 
the  optical  radiation  as  the  variable:  f(x)  is  the  energy  per  unit  wave¬ 
length  at  wavelength  x.  To  make  a  measurement  of  nonzero  energy  from  this 

density  function  it  is  necessary  to  integrate  f  over  some  range  of  wave- 

x+L 

lengths,  for  example,  g  (x)  =  f(y)dy.  We  could  also  make  the  second 

1  'x-A 

x+v^2A 

measurement  g  (x)  =  f(y)dy.  A  linear  scaling  of  the  wavelength  vari- 

2  J  X-V2L 

ables  converts  these  integrals  into  exactly  the  example  in  Fig.  2.  In  this 
interpretation,  the  peaks  in  the  unknown  spectrum  f  have  a  separation 
smaller  than  the  lengths  of  wavelength  intervals  over  which  we  integrate. 
Consequently,  the  data  functions  and  g2  retain  no  obvious  Indication  of 

two  peaks. 
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Of  course,  the  classical  deconvolution  problea  Is  to  recover  the  unknown 
Input  function  f  from  only  one  data  function,  say  gg  »  pQ*f.  The  term 
multiple  operator  deconvolution  is  used  to  emphasize  that  in  recovering  f  we 
may  use  more  than  one  data  function.  In  fact,  part  of  the  problem  is  deciding 
how  many  are  sufficient.  More  precisely,  the  issue  is  how  many  of  which  con¬ 
volvers  p^  are  sufficient  for  the  recovery  of  f.  As  we  shall  see  later, 
the  two  we  have  chosen  for  our  example  in  Fig.  2  are  in  fact  sufficient.  They 
are  in  fact  sufficient  in  a  very  strong  manner.  Not  only  is  f  uniquely 
determined  by  and  g^,  but  in  fact  there  exist  linear  operators  given  by 

convolution  with  distributions  and  v both  with  support  no  greater 

than  the  larger  of  the  supports  of  and  p^,  such  that 

f  *  *  Vg2 

We  will  give  the  explicit  formulas  for  v ^  and  i>2  later.  Here,  to  finish 
our  example  and  discussion  of  Fig.  2,  we  have  shown  to  the  right  in  Fig.  2  the 
result  of  recovering  f  from  gj  and  for  a  digital  simulation  [14): 

discrete  versions  of  "deconvolvers"  i/j  and  were  constructed  and,  using 

discrete  samples  from  the  data  functions  gj  and  g  a  discrete  approxima¬ 
tion  of  f  was  constructed  by  the  sum  of  discrete  convolutions.  A  quite  good 
reconstruction  was  obtained.  (We  attribute  the  asymmetrical  ringing  in  the 
result  to  a  combination  of  truncation  error  and  the  use  of  a  loop  in  the  simu¬ 
lation  that  recalculated  a  quantity  that  should  have  been  calculated  only 
once. ) 

We  would  like  to  emphasize  one  final  point:  the  example  is  based  on 
mathematical  results  regarding  relationships  between  convolution  operators  and 
smooth  functions  on  Rn.  Whenever  we  mention  the  use  of  discretization,  it  is 
always  in  the  sense  of  a  discrete  approximation  to  a  smooth  function.  When  we 
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refer  to  discrete  date,  we  aesn  staples  of  the  underlying  data  functions.  The 
reconstruction  of  the  unknown  input  depends  directly  upon  how  well  the  dis¬ 
crete  data  can  be  used  to  approxlaate  the  data  functions  for  a  given  choice  of 
interpolating  functions.  In  Fig.  2  It  should  be  clear  that  if  gj  and 
were  sampled  at  a  rate  of  one  sample  per  unit  interval  12~  tick  marks),  one 
could  hardly  expect  an  accurate  estimate  of  f.  The  sampling  rate  for  the 
results  on  Fig.  2  was  50  samples  per  unit  interval  (20  samples  per  tick  mark 
interval ) . 

Thus  we  always  view  the  theory,  the  applications,  and  the  algorithms  in 
the  following  order.  First  the  problems  and  results  are  stated  in  terms  of 
continuous  domains  that  model  the  applications.  The  accuracy  of  the  algorith¬ 
mic  implementation  of  the  results  then  is  understood  to  depend  on  the  sampling 
rate,  with  convergence  as  sampling  rates  increase.  The  specification  of  a 
possible  minimum  sampling  rate  for  a  given  reconstruction  problem  is  typically 
a  problem-specific  task  in  interpolation  error  estimates.  In  fact,  the 
Nyquist  sampling  rate  is  sufficient  for  band  limited  functions  f. 

The  example  in  Fig.  2  illustrates  why  equation  (1)  is  of  interest  for 
applications  in  which  the  convolver  p.  must  correspond  to  a  physical,  analog 
device  wherein  the  impulse  response  is  dictated  by  a  solid  state  or  biological 
process.  It  is  frequently  possible  to  select  such  analog  convolvers  which 
satisfy  approximately  the  multiple  operator  criteria  such  that  equation  (1) 
will  hold.  Then  each  associated  deconvolver  can  be  digitally  implemented. 

The  fact  that  the  deconvolvers  act  linearly  and  have  compact  support  means 
that  their  implementation  is  straightforward;  their  action  as  continuous 
linear  operators  implies  stability.  Most  importantly,  the  evident  high  band¬ 
width  of  the  overall  operator  Is  accomplished  without  any  essential  change  in 
the  response  functions  of  the  analog  devices.  The  term  overall  operator 
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refers  lo  the  operator  given  by  the  kernel  distribution  £  p  »v  -  6.  Of 

1-1  1  1 

ccise,  because  of  practical  constraints  such  as  analog  and  digital  approxima¬ 
tions  and  computation  time,  the  design  objective  for  the  overall  operator 
would  not  be  the  identity  operator  with  impulse  response  6  but  rather  a  high 
bandwidth  approximation  of  the  identity  operator  given  by  an  impulse  response 
In  terms  of  the  distributions  in  equation  (1),  since  convolutions  commute, 

m  m 

^Tcp^J-Vj  =  =  <f>.  (2) 

i=l  i=l 

In  a  sense  <p  can  be  considered  to  be  made  up  of  ''parts,"  each  of  which 
arises  from  one  of  the  practical  constraints  just  listed,  along  with  a  special 
part  that  is  deliberately  added  to  control  the  noise  power  spectrum  of  the 
output  of  the  overall  operator. 

The  publications  on  this  subject  have  appeared  primarily  in  the  mathema- 
cal  literature.  The  following  issues  regarding  (1)  have  been  addressed: 
sufficient  conditions  for  the  existence  of  solutions  [1)  -  [4],  [16];  examples 
of  sets  of  distributions  that  satisfy  the  sufficient  conditions  [5]  -  [7]; 
construction  of  explicit  solutions,  that  is,  explicit  formulas  for  the  decon¬ 
volvers  [7]  -  [9];  and  construction  and  evaluation  of  approximate  solutions 
[9],  [11].  [14]. 

Only  recently  have  specific  applications  of  (1)  been  mentioned.  The  work 
of  Berenstein,  Kr ishnaprasad,  and  Taylor  [14]  was  one  of  the  first  times  that 
(1)  and  contemporary  mathematical  methods  for  understanding  the  equation  were 
applied  to  physical  problems.  This  work  also  discussed  the  question  of 
additive  noise  and  the  question  of  the  continuity  of  the  overall  operator  with 
respect  to  the  distributions  p,.  p0,...,  p  .  The  noise  question  is  in  regard 
to  noise  added  following  the  action  of  the  operators  defined  by  the  p  , 
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while  the  continuity  question  is  in  regard  to  the  dependence  of  the  overall 
performance  on  either  the  actual  analog  approximations  of  the  or  the 

digital  approximations  of  the  i>j. 

The  approximation  methods  of  (111  were  motivated  by  this  work  of 
Berensteln  et  al.  These  methods  exploit  the  approximation  in  (2).  In  con¬ 
junction  with  the  analysis,  a  computer  simulation  for  F2  was  performed. 

This  simulation  dramatically  illustrated  (2)  for  imaging  devices  in  which  the 
analog  convolvers  were  solid  state  photodetectors.  With  these  results  there 
was  an  increased  interest  in  imaging  applications.  This  led  to  the  considera¬ 
tion  of  not  just  detectors  but  of  linear  systems  consisting  of  sequences  of 
operators  with  each  operator  of  the  multiple  operator  type.  These  activities 
led  to  the  need  to  answer  basic  systems  analysis  questions. 

In  what  follows  three  topics  in  multiple  operator  deconvolution  are  dis¬ 
cussed.  The  first  is  that,  in  a  sense  that  is  relevant  for  applications,  mul¬ 
tiple  operators  are  necessary  for  the  deconvolution  problem  to  be  well-posed. 
We  present  examples  to  illustrate  the  i 1 1-posedness  of  single  operator  decon¬ 
volution,  and  a  theorem  which  shows  that,  except  for  an  uninteresting  case,  m 
must  be  greater  than  one  in  our  problem  statement. 

The  second  topic  is  that  of  identifying  convolvers  p, . p  for  which 

1  m 

equation  (1)  holds  and  the  construction  of  the  deconvolvers. 

The  final  topic  is  the  major  one  here:  measuring  the  utility  of  a  mul¬ 
tiple  operator  design.  We  describe  the  result  of  our  application  of  standard 
methods  of  linear  systems  and  random  signals  to  the  multiple  operator  type  of 
system  of  equations  (1)  and  (2).  This  analysis  was  necessary  if  one  was  to 
seriously  consider  multiple  operator  designs.  While  the  extended  bandwidth 
was  well  understood,  analyzed,  and  even  illustrated  in  simulations,  the 
consequence  of  the  Introduction  of  noise  and  of  design  errors  was  not  fully 
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understood.  It  was  clear  that  since  the  operator  was  1  Inear  and  continuous 
that  there  would  be  no  lnstabllty  due  to  noise  (at  least  for  smooth  approxima¬ 
tions),  which  Is  already  an  Improvement  over  the  case  of  single  operator 
reconstruction  methods  (5),  1 14],  115].  However,  the  performance  needed  to 
be  expliclty  described  so  that  standard  tools  such  as  resolution,  equivalent 
bandwidth,  and  signal  to  noise  ratio  wou)  .  be  available  for  systems 
engineering  design  studies. 

This  Investigation  was  motivated  in  large  part  by  the  potential  appli¬ 
cation  of  these  multiple  operator  methods  to  electo-optics,  especially  to 
imaging  devices.  We  have  in  mind  imaging  devices  that  are  for  the  detection, 
transformation,  and  display  of  electromagnetic  radiation  for  a  human  observer 
as  well  as  such  devices  for  artificially  intelligent  "observers."  The  prob¬ 
lems  and  the  desired  solutions  have  the  flavor  of  this  application.  While  the 
analysis  and  the  results  are  in  a  sense  general,  much  is  framed  and  guided  by 
the  motivating  problems. 

This  interest  in  electro-opt ics  is  made  explicit  in  the  last  sections  of 
♦his  paper.  There  we  examine  convolvers  and  additive  noises  of  the  type 
encountered  in  imaging  detector  arrays.  The  method  developed  here  for  general 
performance  comparisons  is  used  to  compare  conventional  detector  arrays  with 
arrays  configured  such  that  equation  (15  holds  (up  to  the  resolution  limita¬ 
tion  imposed  by  a  choice  of  optical  sam.pl  ing  rate).  In  terms  of  the  familiar 
parameters  of  resolution,  modulation  transfer  function,  and  signal  to  noise, 
we  show  that  the  performance  of  the  strongly  coprime  design  exceeds  that  of 
the  conventional  design.  Moreover,  the  comparison  has  a  quantified,  function¬ 
al  form  suitable  for  system  design  trade-off  studies. 
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2.  THE  NECESSITY  OF  MULTIPLE  OPERATORS 


Recall  our  mathematical  problem  statement  from  the  Introduction:  under 

what  conditions  do  the  compactly  supported  distributions  p,,p0,...,p  on  Rn 

l  c  m 

have  associated  to  them  compactly  supported  distributions  u..u_ . v  such 

l  c.  m 

ii 

that  J)  p  *i/.  =  A? 

1  =  1 

Here  we  offer  some  observations  and  a  theorem  which  motivate  this  state¬ 
ment.  A  brief  demystification  of  "compactly  supported  distributions"  Is  given 
In  the  Appendix. 

In  Fig.  1(a)  and  in  the  problem  statement,  the  integer  m,  the  number 
of  distributions  p.  ,  is  not  specified  beforehand.  The  number  of  distribu¬ 
tions  els  well  as  the  distributions  i>  ,  ...,y  are  parts  of  the  problem.  In 

1  m 

the  example  above  for  F  it  was  stated  that  m  =  2  sufficed.  When  the  u 

l 

are  characteristic  functions  of  cubes  on  Fn.  then  there  are  choices  for  the 
Pj  such  that  m  =  n  +  1  is  sufficient. 

Fig.  1(b)  represents  a  different,  more  restricted  problem:  it  repre¬ 
sents  a  restriction  of  the  problem  statement  above  to  m  =  1.  We  show  below 
that  there  is  only  one  trivial  case  in  which  a  single  compactly  supported  dis¬ 
tribution  pQ  has  an  associated  compactly  supported  distribution  i>Q  such 
that  Mq'I'q  =  5.  Before  presenting  this  general  result  we  briefly  review 
three  well  known  observations  which  help  to  motivate  the  problem  statement. 

For  these  remarks  let  pQ  be  the  characteristic  function  of  the  interval 
(-1,1)  In  !R  (Pq  takes  the  value  one  or.  the  interval  and  zero  elsewhere). 

Observation  1.  There  can  be  no  such  vQ  if  we  permit  the  input  signals  f 
to  belong  to  any  space  that  contains  C(F),  the  continuous  functions  on  F. 

If  such  a  existed  then,  for  any  continuous  f.  i>0»(pQ*f)  =  f.  But 

for  f  continuous,  pQ«f  is  (represented  by)  a  continuous  function  and. 
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tt  Is  ssslly  checked  that  P0*fn 


hence,  for  f  (x)  ■  sln(nxx),  n  »  1,2,..., 
n 

Is  zero  everywhere. 

The  problematic  functions  f^  in  Observation  1  can  be  eliminated  by  con- 

2 

sidering  only  functions  in  L  (R).  For  in  this  case  we  can  consider  Fourier 

transforms  of  f,  of  (pQ  is  in  L^R)),  and  of  n0»f  (in  L2(R)). 

Let  un  denote  the  Fourier  transform  of  u_.  Since  u_(u)  =  2  -*n  u,  u  €  R, 

0  0  0  u 

then  ( Mq* f )  =  Mpf  determines  f  almost  everywhere.  Since  the  Fourier 

2  2  ' 

transform  is  an  isometry  from  L  (R)  to  L  (R),  f  is  uniquely  determined  by 
Pg»f.  But  still: 

Observation  2.  There  can  be  no  such  v ^  even  if  ve  restrict  the  input 

2 

signals  f  to  belong  in  L  (R). 

If  such  a  existed,  then  let  be  the  support  of  u  .  Let 

Zu  ~  be  the  characteristic  function  (-(M+2) , M+2)  and  let  f  (x)  =  sin  nnx, 
n+c  n 

n  =  1,2,...  .  Then  M0»(fn*M+2)  is  2ero  °n  ( "(MM ) ,  M+l ) ,  hence 
W(fn*M*2)  is  zero  on 

The  difficulty  in  Observation  2  is  due  to  the  fact  that  we  seek  i>Q  with 
compact  support.  What  if  we  dropped  this  requirement?  This  would  be  proble¬ 
matic  for  applications,  for  then  to  get  better  estimates  of  f  from  M0«f  we 
would  have  to  process  larger  and  larger  subsets  of  the  domain  of  p0»f.  But 
ever,  if  we  could  tolerate  noncompactness,  there  is  still  a  problem  of  bounded¬ 
ness  or  continuity. 

Observat ion  3.  There  can  be  no  ( non-compact  1 y  supported)  continuous  linear 

2  2 

operator  L  from  L  (R)  to  L  (R)  such  that  L(MQ«f)  =  f  for  all  f  in 
L2(R). 

p 

The  difficulty  is  that  convolution  with  does  not  carry  L  (R)  Into 

2 

L  (R).  For  example,  let  ^  be  the  approximate  identity  defined  by  ^(x)  = 

x  €  *•  Then  the  l1  norB  of  *c  ls  *  l’  whlle  *^cl2  =  y 
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2  2 
for  the  L  nor*,  and  #c*f*0  converges  to  pQ  in  L  .  If  L  existed  es 

described,  then  by  continuity 

L(p  )  «  L(  1 1®  *pQ)  *  He  L(*  *p  )  «  lie 

c-*0  c-K)  u  c-*0  c 

These  observations  illustrate  ways  in  which  the  problem  of  inverting  the 
convolution  or  solving  the  convolution  equation  pQ»f  =  g  f°r  a  specific  pQ 
fails  to  be  well-posed  (that  is,  that  the  solution  f  to  pQ*f  c  g  exists, 
is  unique,  and  depends  continuously  on  g)  even  with  a  variety  of  restric¬ 
tions  on  f.  Our  problem  statement  in  part  seeks  to  find  m  and  p  ,  .  ..,p 

l  m 

such  that  the  inversion  of  the  simultaneous  distribution  equations  p^*f  * 

g, . p  • f  =  g  is  well-posed,  along  with  an  extra  condition  on  the  form  of 

I  m  m 

the  inverse.  But  possibly  our  observations  were  due  entirely  to  an  unfortu¬ 
nate  choice  for  the  distribution  p^.  Possibly  there  is  a  compactly  supported 
Pj  for  which  there  exists  a  compactly  supported  distribution  v ^  such  that 
^l"Ul  =  lssue  Is  settled  by  the  following  theorem. 

Theorem.  Let  p^  and  be  compactly  supported  distribut  ions  on  Fn  such 

that  p^*^  =  5,  the  Dirac  delta  distribution  or  unit  impulse  (at  the  origin 

in  Then  there  exist  a  e  Fn  and  C  e  C,  C  *  0,  such  that  p.  =  C5  , 

1  a 

the  Dirac  delta  distribution  translated  to  a  e  Fn. 

Proof.  See  Appendix. 

These  observations  and  this  theorem  are  primary  examples  of  the  difficul¬ 
ties  that  are  avoided  whenever  we  can  use  multiple  operators  p  ,...,p  for 

1  m 

which  equation  (1)  holds.  Of  course,  there  is  a  vast  literature  and  many 
approaches  to  address  the  difficulties  of  ill-posed  problems  such  as  the 
inversion  of  a  single  convolution  equation  1151,  117)  and  (18].  Our  interest 
is  in  exploiting  those  cases  in  which  using  multiple  operators  we  have  a  well- 
posed  inverse  problem.  In  the  next  section  we  identify  some  of  those  cases. 
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3.  EXISTENCE  AND  CONSTRUCTION  OF  DECONVOLVERS 


Our  problem  statement  Is  typically  viewed  In  two  parts.  The  first  part 

Is  existence,  the  Identification  of  sets  of  convolvers  pj . p^  such  that 

deconvolvers  exist.  The  second  part  Is  the  explicit  representa¬ 

tion  of  the  deconvolvers. 


These  problems  are  addressed  using  methods  from  the  study  of  analytic 
functions  of  several  complex  variables.  For  if  p^  Is  the  distribution  on 
Rn  with  compact  support,  then  its  Fouriei — Laplace  transform  p^  is  an  analy¬ 
tic  function  on  C°.  Moreover,  Ip^Cz)!,  z  e  Cn  increases  for  large  |z|  = 

1/2 

£  Iz. I  j  in  a  manner  which  completely  characterizes 

these  analytic  functions  which  sire  transforms  of  distribution  of  compact 
support  (the  Paley-Wiener-Schwartz  theorem,  see  Appendix).  Thus  for  a  given 


set  of  convolvers  p, . p  ,  by  taking  transforms  the  existence  of  the 

l  m 

deconvolvers  is  equivalent  to  the  existence  of  solutions  i>,,...,v  of  the 

1  m 

analytic  Bezout  equation 


m 

^p.fzli^Cz)  =  1,  z  €  Cn,  (3) 

i  =  l 


with  y . u  each  in  the  Paley-Wiener  class  of  functions. 

l  m 

In  this  form  one  quickly  sees  that  a  necessary  condition  on  the  convol¬ 
vers  p^  is  that  they  have  no  common  zeros.  Moreover,  because  the  are 

to  satisfy  certain  growth  conditions  for  large  \z\,  one  has  a  stronger 
condition  which  turns  out  to  be  both  necessary  and  sufficient. 

Theorem  [1),  { 16).  For  the  compactly  supported  distributions  p  ,...,p  on 

1  m 

Rn  there  exists  compactly  supported  distributions  v , , . . . ,  v  such  that 

1  in 

in 

E  p,*v,  *  S  if  and  only  if  there  exist  positive  constants  c.  and  c_  and 
1*1  1  1  12 
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a  positive  Integer  H  such  thot 


(z)|2  *  Cle’C2,I“  z,(l*|zlfN.  z  c  Cr 


A  set  of  convolvers  that  satisfies  the  inequality  in  the 

theorem  is  often  referred  to  as  strongly  coprime. 

The  following  sets  of  convolvers  are  known  to  be  strongly  coprime: 

2 

i)  on  R  ,  the  characteristic  functions  of  two  disks  wherein  the  ratio 
of  the  radii  of  the  disks  has  any  integer  value  from  2  to  200 


ii)  on  R  ,  the  characteristic  functions  of  n+1  cubes  with  side 

2  2 

lengths  s. . s  ,  such  that  s,,...,s  .  are  integers  that  are 

1  n+ 1  1  n+ 1 

pairwise  relatively  prime  and  at  most  one  of  these  integers  is  a 
perfect  square  of  integers  (7),  Ill]. 

The  example  of  the  introduction  falls  within  case  ii.  Another  case  which 
is  known  is  the  characteristic  functions  of  certain  rotations  of  three  squares 
in  the  plane  [7 ] . 

Explicit  formulas  for  the  deconvolvers  have  been  found  for  the  disks  on 

2 

F  (case  i)  and  for  n  =  1,  and  n  =  2  in  case  U  [5]  -  [9]  and  [14].  For 
example,  the  explicit  solutions  for  the  deconvolvers  for  n  =  1  in 

case  ii  is  the  following.  For  i  =  1,2,  a^  >  0,  let  be  the  character¬ 

istic  function  of  the  interval  (-a^.a^),  with  a^  =  V^a^.  Let  #  be  any 
smooth  function.  Then  and  v 2  such  that  =  5  are  given 


2  3 

V*  =  8a~ar*2“*  *  24a~  a2  ~2  (*2**5  +  73  ( ( *2K1  * *** 
12  1  dx  dx 


where 
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and  v  and  are  given  by  the  permutation  of  indices.  Note  that  the 

support  of  i/j  is  contained  in  that  of  z2»  8111(1  tlie  support  of  v 2  is 

contained  In  that  of  Zj. 

Throughout  the  foregoing  we  have  emphasized  that  we  seek  deconvolvers  v 
which  have  compact  support,  and  we  have  described  necessary  and  sufficient 
conditions  on  the  convolvers  for  these  to  exist.  The  explicit  determina¬ 

tion  of  whether  the  p^  are  strongly  coprime  and  the  explicit  construction  of 
the  deconvolvers  is  typically  too  restrictive  an  approach  for  general  engin¬ 
eering  applications,  and  new  numerical  optimization  approaches  are  being 
developed  [10]. 

There  is,  however,  a  middle  ground  that  we  use  when  we  have  a  set  of 
strongly  coprime  convolvers  but  do  not  have  explicit  representations  for  com¬ 
pactly  supported  deconvolvers.  These  alternatives  are  due  to  the  fact  that 
the  deconvolvers  are  not  unique.  To  see  this  nonuniqueness,  consider  the  one¬ 
dimensional  example  of  the  Introduction.  Let  p^  and  be  compact 

supported  and  strongly  coprime,  and  let  and  u be  compactly  supported 

deconvolvers.  Let  X  be  any  compactly  supported  distribution.  Then  it  is 

readily  seen  that  both  v '  =  u.  +  X»p^  and  v '  =  i>  -  X»p^  are  compactly 
2 

supported  and  £  p .•!>'  =  6.  Going  further,  it  was  not  necessary  that  X  be 
i  =  l  1 

a  compactly  supported  distribution.  That  is,  we  can  construct  deconvolvers 
that  are  not  compactly  supported. 

We  frequently  use  approximate  deconvolvers  even  when  they  are  not  com¬ 
pactly  supported.  For  example,  let  f>  be  any  integrable,  smooth  function 
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lx 

with  integral  over  It  equal  to  1.  Then,  for  •  (x)  »  -*(-),  ».*•  is  an 

T  T  T  IT 

approximate  deconvolver  and  converges  to  as  x— >0.  We  can 

1 

choose,  also,  that  the  v^mfr  be  L  functions.  Let  p  have  compact 
support  and  be  sufficiently  differentiable.  (The  function  p  is  usually 
called  a  mollifier. } 

Finally,  we  can  construct  approximate,  non-compact ly  supported  deconvol¬ 
vers  on  Rn  solely  from  the  knowledge  that  the  convolvers  are  strongly 

coprime.  Let  u~  )  =  Re  z  *  (Re  z,,Re  z_, ...,Re  z  )  and  for  what 

in  l  c.  n 

follows  we  will  need  to  consider  the  Fourier  transform  only  on  the  real  sub¬ 
space  Rn  of  Cn.  Since  e  C^tR11),  and  with  p  €  C^R11),  with  p  com¬ 
pactly  supported,  let 

a 

D,  (w)  =  - ,  h.  =  D.(u)p(u),  i  =  1 . m.  (5) 

l  m  *  ?  i  i 

l  iM,(w)r 

In  n 

This  defines  h.  c  L  (R  }  such  that  J]  h.  »p.  *  p  (z  denotes  the  complex 

i  =  l  1 

conjugate  of  z).  This  is  the  class  of  deconvolvers  that  are  discussed  in  the 
last  section,  because  from  among  these  one  can  find  optimal  deconvolvers  when 
there  is  additive  noise. 


4.  PERFORMANCE  OF  OPTIMAL  DECONVOLVERS. 

While  (5)  is  exhibited  essentially  by  inspection,  the  result  can  be 
obtained  in  a  more  systematic  fashion  as  well  as  in  a  more  general  form.  We 
first  recall  some  standard  tools,  apply  these  tools  to  a  simple  case,  and  then 
proceed  to  the  more  general  form.  The  diagram  in  Fig.  3  represents  an 
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operator  L  acting  on  a  function  f.  Let  (temporarily)  f  be  bounded  and  in 
C**(Rn).  Let  Hy  p2 . be  an  arbitrary  set  of  a  distributions  with  com¬ 
pact  support.  For  each  linear  operator  defined  by  p^  let  be  a  sample 

function  of  a  zero  mean,  wide-sense  stationary  random  process  that  is  added  to 
the  output  of  p^,  let  ijj  €  L*(Rn),  and  let  (Nj  2  0)  be  the  noise 

power  spectral  density  of  the  process  (see,  for  example,  119)  Ch.4,  Ch.6). 

For  each  distinct  i  and  J  let  tj  be  independent  of  t>  and  let  each  i) 

*  J  J 

be  independent  of  f.  Let  u.  be  defined  by  =  D^,  where  0,  e 

p  n  * 

C  (F  ),  ^  has  compact  support,  and  r  is  sufficiently  large  so  that 

^D.*j  e  L!(Rn).  Let  g  €  L"(Rn)  be  defined  by 

m 

g  =  Lf  =  ^(p^f  +  (6) 


In  the  usual  manner,  with  E  denoting  expectation. 


E(g}  =  £Vf.CV„. 


Let  denote  translation  by  y,  r(x)  =  x  ♦  y,  let  '  denote  inverse 


Fourier  transform,  and  let 


denote  the  LH  norm.  Directly  from  the 


definition  of  wide-sens>t  staL i un^ry  and  n^ise  power  spectral  density  it 


follows  that 


•  m 

E|(g-E{g>)^(g-E{g)).rJJ  «  XN?IDi|2'^2  (y)* 

4  =  1 


and,  for  y  »  0,  that 


K'8'^8"2}1  ^Th|Iniidi|2|J|2 
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The  simplest  configuration  for  L  is  all  distributions  equal,  all  decon¬ 
volvers  trivial,  and  all  randoa  processes  Identically  distributed: 

Mj  =  Wj  *  6.  *  N^.  for  1  *  1,2 . a.  (10) 

Then 


E<g>  =  m»0***f.  E|(g-E(g>)2J  =  _J?_||N0lil||2.  (11) 

The  utility  of  (7)  -  (9)  or  of  (11)  is  that  if  L  is  followed  by  a 

linear  operator  II  with  kernel  u  (which  could  model  a  specific  "end-lise") 

2 

then  one  cam  compare  the  function  (li(E(g>))  with  the  constant  function 
2 

E(  (U(g-E(g> ) )  ).  In  the  case  of  the  simplest  configuration,  (10)  and  (11), 
there  are  the  following  formulas  and  bounds. 


(U(E<g>))2  =  EOig}2  =  (u.(mp0.^.f))2 


(2  n)n 


A  A  A 

llup0*f|l 


f-1"-  -J2  Ilf  11?  when  f  e  L2(Rn); 

l  (2n)nJ  0  2  2 


(12) 


and 


eT ( U ( g  -  E{g> )  )21  =  — !1_|!uN  ill2 

J  (2n)n  U  * 


(13) 


p 

The  function  E(Ug)  is  referred  to  as  the  signal,  its  square  E(lig>  is 

2 

referred  to  as  the  signal  power  or  energy,  and  E{ (U(g-E(g) ) )  )  is  referred 
to  as  the  noise  power.  Typically  the  ratio  of  E{Ug}2  to  E(  (lKg-E(g) )  )2) 
is  considered  or.  alternatively,  the  positive  square  root  of  the  ratio.  Here 
we  shall  consistently  use  the  latter.  If  this  ratio  is  evaluated  at  some  dis¬ 
tinguished  point,  the  value  defines  a  "signal  to  noise  ratio."  We  denote  by 
P  the  evaluation  of  the  absolute  value  of  the  function  at  so.ne  distinguished 
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point.  CIven  L  and  for  a  given  choice  of  #,  f,  1,  and  $  define  the 
signal  to  noise  ratio 

W*(UL)  *  - P-U(E(-g}J - -rjj.  (14) 

^(u(g-E{gn)2JJ 

For  a  fixed  choice  of  0,  f,  11,  and  J>.  two  operators  L  and  L'  can  be 
compared  and  ordered  by  (14). 

On  the  other  hand,  for  a  choice  of  f,  U,  and  J),  (14)  is  determined 

for  the  case  of  the  trivial  operator  in  (10)  by  the  pair  of  functions 

mpg  and  i/in  Nq.  (15) 

In  general,  let  operators  L  and  L'  (for  example,  as  in  Fig.  3)  have  trans¬ 
fer  functions  and  noise  power  spectral  densities  p,  N2  and  p' ,  N'2, 
respectively.  For  a  choice  of  U  we  shall  say  that  UL|UL'  (i.e.,  "UL 
divides  UL' " )  if  there  exists  a  "quotient"  function  q  €  L°°(IRn)  such  that 
aqp  =  ip'.  If  UL|UL'  and  S  |u|2N'2,  we  say  that  UL  >  UL' . 

This  definition  is  motivated  by  the  following.  As  usual,  let  4>  be  such 
that  a  linear  operator  Q  with  kernel  q  can  be  associated  with  q  by 

considering  q <f>.  Let  8  be  any  continuous,  translation  invariant,  linear 

!fA7(  ( 1L5QL ) 

operator.  For  fixed  11  if  UL  >  UL' ,  then  TlIsL *  )  “  *  '  Consequently, 

sup  yA'fldiSL)  >  sup  yAX(USL'). 

8  8 

Next  consider  the  operator  L  diagrammed  in  Fig.  3  for  the  case  in  which 

p,,  P-.....P  are  distinct  and  strongly  coprime  (i.e.,  satisfy  (4)).  An 
1  c  m 

m  .  2 

obvious  consequence  is  £  Ip  (w)|  >  0  and,  equivalently, 

i  =  l  1 

0  *  (p  (w) , p_(w) , . . . , p  («))  €  Cn.  u  €  Fn.  (16) 

1  c  m 

Consequently  we  can  visualize  (16)  as  is  shown  in  Fig.  4a.  A  similar 
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Illustration  can  be  used  to  visualize  flvHfijCtO.Pglw) . .  {«))  » 

*  *  a  *  Sk  tk 

(f(w)p1(w),f(w)p2(u) . Tlu)fijiu)),  except  the  "curve"  passes  through  the 

origin  if  and  only  if  f(w)  *  0.  The  power  spectral  densities  are  real  and 

2 

nonnegative  (thus  we  write  and  choose  1  0).  Assume 

Nj (w)  >0,  w  c  Fn.  1  -  1.2 . a  .  (17) 


We  can  visualize  (17)  as  is  shown  in  Fig.  4b.  The  case  of  strongly  coprime 

multiple  operators  has  the  useful  feature  that  the  consideration  of  (16)  and 

(17)  pointwise  in  conjuction  with  (7)  -  (9)  uniquely  determines  an  alternative 

choice  for  the  Di  of  (5).  This  choice  will  be  optimal  in  the  sense  it  has 

2 

the  smallest  E((g-E(g})  >  among  all  sets  of  deconvolvers. 


Proposition.  Let  N.  e  L°°(Rn) ,  N  (w)  >  0  for  u  €  Rn,  1  =  1,2 . m.  Then 

there  is  a  function  D  :  Rn— »Cm  uniquely  determined  (almost  everywhere)  by 

the  condition  that  for  each  u  e  Rn,  D(«)  =  (D,l«),D_(u) . D  (u))  solves: 

l  d  m 


Minimize 


|2N?(w) 


on  the  set 


I  Vi-  = 


In  fact 


D.  (u) 

l 


Pj  (w) 

N2(w) 

m  |p  (w)|2 

i—i — 

i=l  N  (u>) 


(19) 


Proof.  Any  z  that  satisfies  (18)  is  clearly  contained  in  the  linear 
subspace  of  Cm  determined  by  the  span  of 

<(*!.(«).  0 . 0),  (O.M-(w).0 . 0) .  (0,0.  ...0,p  (w))>.  (20) 

i  c.  m 

That  is,  2^  *  0  if  Pj(w)  «  0.  Equivalently,  there  exists  A  « 
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Then  (18)  reduces  to 


Let  £' 
1 


denote  £ 

.1  =  1 

,  |p  (u>)  |2 

Minimize  £'  |Xj  |2 1^ (w)  |2  on  — jSnJT 


(22) 


From  this  it  follows  that  the  X^  are  all  real,  so  that  (22)  has  the  form 


.  ?  f  *  iMi  (WM  1 

Minimize  £'  (Xj  |p.  ( u) \y  on  |E' I (“)  I  —  *l|- 


(23; 


With  the  new  variable  X^|p.(w)|,  it  is  elementary  to  see  that  (23)  has  the 
unique  solution 


(for  MjJu)  *  0).  (24) 


Consequently,  from  (21),  the  unique  z  corresponding  to  the  minimum  is  D(d) 
as  in  (19).  ■ 

In  addition  to  hL  >  0,  i  =  l,2,...,m,  we  shall  assume  NQ  >  0. 

Further,  we  shall  assume  that  the  are  sufficiently  differentiable  and 

that  ~  =  0( |u | p)  for  some  integer  p,  i  =  0,1,2,.. . , m  (that  is,  1/N 

1  1 

_  A  _  / 

does  not  grow  faster  than  |u|p).  With  this  we  can  find  0  =  0(|u|  p  )  so 
that  (D^0)  €  L2(Fn)  and  for  $  sufficiently  smooth  and  with  compact 

support  then  (D^)’  e  L^F0). 

Corol larv.  For  the  choice  of 


Ai |Mi (w) |  = 


|M.(u) j 

Ni  ( w) 


m  |fi  («))' 


j=l  N^(w) 


from  the  Proposition, 


(25) 


Let  Ln  identify  the  trivial  configuration  of  L  in  (10)  and  let  L 

identify  the  strongly  coprime  configuration.  Unless  explicitly  indicated  to 

the  contrary,  Lg  indicates  that  the  deconvolvers  of  (19)  are  used.  The 

first  of  the  functions  in  (25)  is  the  transfer  function  for  L  and  the 

s 

second  is  the  square  root  of  the  noise  power  spectral  density.  The  corres¬ 
ponding  functions  for  LQ  are  (15).  The  mollifier  <f>  is  suppressed  but 
understood.  From  (25)  obviously  l_s  1 L.  for  any  operator  L.  (With  an  abuse 
of  notation  we  use  L,  LQ  and  Lg  to  denote  configurations  consisting  of 

linear  operators  auid  additive  noises.)  From  (15)  and  (25)  the  quotient  q 
*  2 

for  L  =  Lq  is  mpQ.  Let  Ng  denote  the  noise  power  spectral  density  of 
Lg.  In  the  sense  discussed  earlier  let  denote  the  linear  operator  asso- 

A  A 

dated  with  .  That  Ls  |L^  with  quotient  means  LQ  =  aQLs-  Then 

QgLs  has  functions  corresponding  to  (25)  (transfer  function,  square  root  of 
noise  power  spectral  density)  given  by 


m 


woZ"iDi 

i  =  l 


r  m  "i 

1/2 

r  m  - 

4 

=  mp0.  |mp0|Ns  =  m|M0| 

Zim2"? 

4  =  1  J 

m 

1  z- 

j=l 

h'2  . 

1/2 

vGiN  .  (26) 


By  definition  llLg  2  UL^  if  ImpglhMu)  <  v'mNg(u)  on  the  support  of  u,  and 

UL_  2  11L  if  ULn|UL  and  |mp  |N  (u)  2:  vfaL(u)  on  the  support  of  u. 
u  s  u  s  u  s  u 

Thus,  whether  UL^  2  ULg  or  ULg  2  UL^  holds  depends,  in  part,  on  whether 
one  of  the  following  inequalities  holds  on  the  support  of  u  :  from  (15)  and 
(26) 
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mp0|Ng(w)  ^  }/mN0(v) 


(y  "g<“> 


N?(u)  . 


1/2 


V^b|m0(u)|.  (27) 


In  (27)  the  notation  means  that  the  upper  inequality  symbol  on  the  left  is  to 
be  paired  with  the  upper  inequality  symbol  on  the  right  and  lower  left  with 
lower  right. 

The  comparison  in  (27)  can  in  special  cases  be  viewed  from  a  slightly 
different  perspective.  First,  view  the  left  side  of  the  second  inequality  in 
(27)  as  the  Fourier  transform  of  a  kernel.  Define 


e  (cj)  = 


m  N?(u) 


f-l  N7(u) 


i  =  l  i 


1/2 


(28) 


We  refer  to  c  as  the  envelope  transfer  function  corresponding  to  the 

envelope  operator  f  for  a  given  strongly  coprime  Lg  in  comparison  with  a 

given  Lq.  If  v'mS  acts  on  Lg,  then  the  pair  of  functions  associated  with 

/mf L  is 

s 


Vm  e ,  v'mNg. 


(29) 


Recall  that  the  pair  for  Lg  is  given  by  (15)  (rewritten  for  convenience) 

mp0,  VmNg.  (15) 

That  is,  the  composition  of  /m£  with  Lg  has  a  noise  power  spectral  density 
equal  to  that  of  Lg.  If,  for  example,  Pg  is  real  and  positive,  then  it 
makes  sense  to  compare  (29)  with  (15).  It  is  easy  to  check  that  the  condition 
uv'me  2  umpg  (on  the  support  of  a)  coincides  with  our  definition  ULg  2  ULg, 

A  A  A  A 

and  umpg  2  uv'me  coincides  with  what  we  mean  by  ULg  2  ULg.  These  two 

inequalities  are  precisely  the  content  of  the  comparison  of  the  right  side  of 

(27).  One  could  say  that  r/m?  is  the  normalization  of  L  to  the  noise 

s 
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power  spectral  density  of  LQ. 


For  either  point  of  view  we  consider  W>  ■  {w  €  Fn  :  e(o)  i.  ^ |Mq(w) | } 

_ 

and  W^,  =  <w  €  R  :  e(w)  i  Vm | mq(«)  | >  -  For  all  U  such  that  a  has  support 
in  W  it  follows  from  (27)  and  the  definitions  that  UL  i.  UL-. 

>  5  v 

Consequently, 

y/rximD  |iin| 

_ — —  =  - - — _ -  >  l,  (30 

y/muiL n)  ,  .  Vm  p 

I  U  — <p  N 


where  QQ  is  used  to  denote  the  linear  operator  corresponding  to  the  transfer 
function  mpg  of  LQ. 

Assume  Pg(0)  *  0  and  define 

nQ  =  (w  €  IRn  :  Vt  €  (0,1)  |M0(tw)|  >  0>. 


Note  that  for  R1  the  usual  definition  of  limiting  resolution  is  sup  fiQ.  If 

~  ~  -1 
supp(u)  is  compact  and  supp(u)  c  fig,  then  UDg  makes  sense,  consequently 

ULq |  llL^ .  Hence,  if  supp(a)  is  compact  and  supp(u)  c  W<  n  nQ,  then  ULg  > 

ULs.  Consequently, 


yA'K(ULs) 

yATKUQ^L 


^  Mo 


<  1. 


(31) 


In  general  the  inequality  cannot  be  extended  to  all  of  W<  a  (lg  because  of 
the  behavior  of  1/jig  on  the  boundary. 

y^RUL  ) 

There  is  no  Information  regarding  ^  )"  1  mplied  by  either  ULg  >  UL( 

or  UL-  i  UL  .  Additional  information  is  needed.  For  example,  it  may  be 
0  3 

sufficient  to  know  the  effect  of  the  so-called  "boost"  ULgi — ►  U0g1Lg.  In 
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(32a) 


(32b) 


For  supp(u)  c  Fn-ng,  or  even  for  supp(u)  n  (Fn-ng)  *  0,  it  is  often 

the  case  in  applications  that  UL-  is  not  defined.  Since  UL  is  defined 
^  0  s 

for  all  U  it  malces  sense  in  such  cases  to  consider  UL  >  UL„. 

s  0 


5.  EXAMPLES:  OPTIMAL  DECONVOLVERS  FOR  CHARACTERISTIC  FUNCTIONS  OF  SETS  IN  Fn 

Collections  of  sets  in  IRn  such  that  the  characteristic  functions  of  the 

sets  in  the  collection  are  strongly  coprime  have  been  described  in  Section  3. 

Let  denote  the  characteristic  function  of  a  set  S.  Consider  the 

2 

case  in  which  the  noise  power  spectral  density  has  the  form  IIz^IIjN^,  where 

llZgl'l  is  the  L*  norm  of  Xg  (equivalently,  the  Lebesgue  measure  of  the 

set).  For  such  a  case,  let  sets  S.,S_ . S  ,  be  chosen  so  that,  for  p.  = 

lcm  l 

Xo  ,  the  p  ,p  , .  ...p  are  strongly  ccprime.  Then,  from  the  Proposition  in 
o.  lcm 

the  previous  section, 

Hj(w) 

D.  (u)  =  - "  —~5  (33) 

■  |p  (u)p 

E  — - - 

1=1  lMJll 

and  (25)  becomes 
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(34) 


Let  SQ  be  any  set,  let  p^  =  Its  characteristic  function,  and 

consider  this  to  be  the  convolver  In  Lg  defined  by  (10)  (l.e. ,  m  parallel, 

identical  convolvers).  Let  the  noise  power  spectral  density  have  the  same 
2  2 

form  as  above,  =  ||p  ||  N  .  From  (27)  and  (28)  one  obtains  an  envelope 
U  U  1  0 

transfer  function  and  the  associated  comparison  for  these  two:  a 

1/2 

convenient  renormalization  by  the  constant  KMqIIj  is  made  in 


2 

For  an  explicit  example  let  Si  c  R  be  the  region  in  a  focal  plane  of 

an  imaging  device  which  corresponds  to  a  single  light  sensitive  detector.  The 

exposure  time  interval  is  assumed  fixed  and  the  image  is  assumed  constant. 

Then  p^  =  is  the  idealized  response  function  of  the  detector.  (The 

i 

actual  shape  of  the  response  function  can  be  incorporated  into  the  mollifier 

0  of  Fig.  3.  That  is,  if  a.  is  the  actual  detector  response  function, 

and  if  the  deviation  of  from  p.  is  due  to,  say,  a  material  diffusion 

process  that  is  common  to  all  detectors,  then  =  p.»0,  where  0  models 

the  diffusion.  Such  a  common  0  would  not  be  deconvolved;  one  would  use 
m 

£  p.*0»u  =  0.  )  Then  p  is  what  is  referred  to  as  the  "detector  MTF.  "  The 

1-1  1 

form  of  the  noise  power  spectral  density  above  is  that  of  many  noise  processes 
in  electro-optical  detectors,  that  is,  that  the  noise  power  spectral  density 
is  proportional  to  the  area  of  the  detector.  Such  is  the  case  for  infrared 
detectors  and  this  proportionality  is  contained  in  the  definition  of  the  fami- 


2.24 


liar  "D*."  The  above  fora  is  also  valid  for  the  so-called  background  limited 

3 

case.  It  also  has  this  fora  for  R  when  the  tine  interval  is  included  as 
the  third  dimension.  Further,  a  background  United  slit  detector  corresponds 
to  the  above  forms  for  R1  with  the  slit  width  as  the  coordinate.  (The 
background  limited  case  consists  of  a  relatively  small  signal  of  interest 
superimposed  on  a  relatively  large  constant  signal  so  that  the  noise  in  the 
signal  of  interest  is  due  to  the  "shot**  noise  of  the  constant  signal.)  (For 
detector  characteristics  discussed  above  see,  for  example,  (20,  Ch.2].) 

In  Fig.  5,  Fig.  6,  and  Fig.  7  the  transfer  functions  for  such  cases  are 
shown.  In  Fig.  5,  a  comparison  is  shown  for  the  example  in  R*  described  in 
the  Introduction.  The  characteristic  functions  p^  and  p^  f°r  two 
intervals  (-1,1)  and  (-v' 2,V2 ),  respectively,  are  strongly  coprime.  The 
envelope  transfer  function  is  shown  and  is  compared  with  the  transfer 

function  for  the  two  identical,  parallel  convolvers  as  in  (35)  where  p^  =  p^. 
The  choice  p^  =  p^  is  used  rather  than  p^  =  p2  in  this  comparison  because 
Pj  is  "better"  than  p^  in  the  sense  that  the  first  zero  of  p^  (i.e.,  its 
bandwidth)  is  greater  than  the  first  zero  of  p^.  Recall  from  the  scaling 
property  for  Fourier  transforms  on  R*  that  p^(x)  =  ^^{V2x)  for  all  x  e  R* 
if  and  only  if  v^p^v^u)  =  p2(w)  for  a11  u  €  R*-  In  making  the  comparison 
we  use  two  (i.e.,  m  =  2)  identical  parallel  convolvers  rather  than  only  one 
so  as  not  to  bias  the  comparison  in  favor  of  the  strongly  coprime  case.  In 
terms  of  the  applications  to  electro-optics,  the  m  strongly  coprime  detec¬ 
tors  must  make  m  sequential  measurements  of  the  optical  signal,  and  the 
noise  for  each  measurement  is  assumed  Independent.  Consequently,  there  is  an 
inherent  noise  averaging  due  to  these  m  measurements,  and  we  wish  to  have 
the  number  of  Independent  measurements  the  same  for  both  cases.  The  case  of 
sequential  measurement  with  identical  detectors  for  noise  averaging  is  exactly 
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what  Is  done  In  so-called  "tine  delay  integration"  for  scanned  infrared 
detector  arrays  (21].  Fig.  5  illustrates  the  consequence  of  the  strongly 
coprime  condition:  the  envelope  response  is  approximately  an  envelope  for  the 
modulus  of  the  other  two  responses  and,  correspondingly,  is  without  zeroes. 
Also,  it  can  be  observed  that  the  envelope  response  decreases  approximately  as 

1/lwl. 

In  Fig.  6  and  Fig.  7  the  envelope  transfer  function  is  shown  for  an 
2 

example  in  R  ,  the  case  of  three  squares  of  side  length  l,v*>,  V2, 

respectively.  The  characteristic  functions  of  these  three  squares  are  strong¬ 
ly  coprime.  The  comparison  (35)  is  illustrated  by  graphing  the  modulus  of  the 

2 

corresponding  transfer  functions  for  two  subsets  of  R  :  the  u^-axis 

2  2 
{u  =  (Wj,o>2)  e  R  :  u>2  =  0}  (see  Fi^.  6)  and  the  diagonal  {w  =  e  R 

W1  =  (see  Fig'  7^'  8raPhs  use  the  Euclidean  distance  as  abscissa, 

2  2  1/2 

M  =  (Wj+i^)  •  The  comparison  illustrated  in  Fig.  6  and  Fig.  7  is  for 

A  A 

Pq  =  Xq  •  (As  before,  *q  has  the  greatest  bandwidth  and  the  scaling  prop¬ 
erty  for  Rn  has  the  form  Pj(x)  =  p2(Jcx)  for  k  >  0  and  for  all  x  €  Rn 
if  and  only  if  Anpj(A  w)  =  p^fw)  for  u  6  ^  The  comparison  is 

essentially  the  same  as  that  for  the  two  intervals  in  R*.  The  difference 

between  the  w^-axis  and  the  diagonal  illustrates  that  approximately  the 

-1  -2 

envelope  response  decreases  as  |u|  along  the  u^-axis  and  as  fu|  along 
the  diagonal. 

From  (35)  (and  as  illustrated  by  the  figures)  the  following  statements 
can  be  made.  These  are  stated  as  "observations"  because  the  results  can  not 
be  given  in  terms  of  explicit  inequalities.  Some  notation  is  helpful.  Define 


n  =  {w  €  R  :  Vt  e  (0,1]  |p  (tu)|  >0)  and  (1  =  ft  n  . 

i  =  l 


Observations.  Let  Pj , Pg, • • • . Mn  «  L  (Rn)  be  strongly  coprime  characteristic 
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functions  of  sets  In  Fn  as  considered  above.  With  each  p^  let  there  be 

associated  as  in  Fig.  3  an  additive  wide-sense  stationary  noise  with  noise 

power  spectral  density  of  the  form  IP.I,^-  Let  L  be  the  configuration  in 

110  s 

Fig.  3  with  deconvolvers  determined  by  the  Proposition.  Let  Lq  be  the 
trivial  configuration  as  in  (10)  with  pQ  *  Pj,  NQ  =  Nj. 


Observation  1:  For  U  with  supp(u)  c  £2,  ULg  a  UL^. 

*  n  m 

Observation  2:  For  H  with  supp(a)  c  Fn  -  U  £2, ,  11L  i  ULn. 

i=2  1  S 

.  m 

Observation  3:  For  1  with  supp(u)  c  U  £2.  -  £2,  11L  <  ULn. 

i=2  1  S  ° 

*  -1 

Observation  4:  For  U  with  supp(u)  compact,  supp(u)  c  £2^ ,  let  0^  be 

the  boost  on  U,  llL^i — (see  (31)  and  (32)).  If  Observation  3  can 
be  neglected  then 


y^K(iiLrt) 


fMR  ( UL  ) 
s 


As  discussed  at  (31)  it  is  not  possible  to  extend  this  to  all  of  £2^, 
for  p1  =  0  on  the  boundary  of  £2^  nowever,  it  still  is  desirable  to  have  a 
means  to  compare  Lg  with  the  more  well  known,  more  thoroughly  studied 
trivial  configurati  s.  In  the  next  section  this  is  accomplished  by  pushing 
the  troublesome  set  (p^  =  0}  out  toward  infinity. 


6.  •  MORE  COMPARISONS:  STRONGLY  COPRIME  VERSUS  CHANGE  OF  SCALE 

Let  L  be  the  same  as  above.  In  the  above  L  was  compared  with  L  , 
s  s  o 

2  2 

where  LQ  was  chosen  to  be  Pj  and  Nj  =  llp^N^  •  In  these  cases  p1  was 
the  "best"  in  the  sense  £2l  c  £2j  1  *  l,2,...,m.  Here  Lg  will  be  compared 
with  a  one  parameter  family  of  such  L.  Define  L^  by  the  trivial  configura- 
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Corollary  to  Observation.  For  the  conditions  in  the  Observation  above,  it 
is  advantageous  to  construct  L  using  sets  that  are  as  large  as  possible. 

Application  of  the  Corollary.  In  parallel  scanned  imaging  systems  with  square 
detectors  wherein  the  systems  are  ranked  using  some  U  meeting  the  conditions 
of  the  Observation  (e.g.,  horizontal  or  vertical  bars),  the  detector  size 
should  be  sufficiently  large  so  that  the  array  of  detectors  fills  the  image, 
and  the  detector  sizes  in  the  array  should  constitute  a  strongly  coprime 
collection.  (This  application  depends  on  sufficiently  high  sampling  rates. 

See  Section  1 . ) 


The  Observation  is  illustrated  in  Fig.  8  for  n  =  2.  For  Fig.  8  Lg  is 
as  in  Fig.  6:  in  the  notation  just  above  Lg  is  configured  from  the  parallel 
convolvers  p<1;>,  ^<y/2>’  ^<V3>’  3X101  ls  ^he  characteristic  function  of  the 


unit  square.  For  this  Lg  the  envelope  transfer  function  is  compared 


with  V3  -  as  In  (35),  for  c  =  1,  0.5,  0.2,  and  0.1.  For  a  =  1 

■»*«»■  1 


see  Fig.  6;  for  a  *  0.5,  0.2,  and  0.1,  see  Fig.  8.  The  observation  In 


(40)  is  clearly  evident.  (Here  we  neglect  Observation  3  of  the  last  section  by 


2.28 


■eans  of  a  broad  interpretation  of  *  in  Observation  1. ) 


The  Observation  (40)  depends  on  the  following  properties.  The  first, 
which  is  again  an  approximation,  is  that  for  Aj  *  {«  €  Rn  :  ■  0,  i  *  J> , 

the  Uj-axis, 

edU(w)  “  CM'1-  (' 

J 


The  second  is  that 


l£*>(w)l  n^K(<kW,l 


IM^JI 


a  £  1,  and  for  u  €  A^,  V^n 


1/2  = 

.i  ,,1/2' 

1 

"*t'i 

|mj(^) 1 

i-l 

„  1/2 

<  C|«| 

'“i'i 

Hence,  for  n  i.  2,  for 


™  - 1/2  -  CM  • 

>»*<*>■! 


Fig.  9  and  Fig.  10  show  two  counterexamples  for  cases  not  addressed  in 

2 

the  Observation.  Fig.  9  is  for  the  case  of  the  diagonal  in  R  ,  and  Fig.  10 

is  for  n  =  1.  The  Observation  fails  on  the  diagonal  D  = 

2 

{u  -  e  R  :  u>^  =  u> because 


ed^2)(w)  *  Clwl  • 


It  fails  for  R  because  (41)  holds. 

If  in  place  of  characteristic  functions  of  cubes  one  uses  characteristic 

,/\  , 

2  -  JM<o>l 

functions  of  disks  on  R  ,  then  the  relationship  between  c  and  vm - — — 

®  a> 11 1 

is  intermediate  between  that  of  the  w.-axis  and  that  of  the  diagonal  for 

J 


e  (u>)  £  C  |  w  | 


The  significance  of  Observation  (40)  is  that  it  provides  a  qualitative 
lower  bound  for  the  performance  of  the  strongly  coprime  conf igurat ion.  To  the 
extent  performance  is  characterized  for  the  1IL^,  the  "envelope"  consisting 

of  the  collection  over  all  a  is  a  lower  bound  for  the  performance  of  BL  . 

s 

All  of  the  above  has  focused  on  performance  away  from  the  origin.  If  the 


figures  are  rescaled  so  that  the  m<a>  appear  fixed  with  a  sequence  of  Lg 
constructed  from  convolvers  of  increasing  support,  the  Observation  indicates 
that  nothing  is  sacrificed  away  from  zero  while  the  envelope  transfer  function 
near  zero  is  substantially  increased.  That  is,  lLL-s  £  represents  a  sub¬ 

stantial  enhancement  near  u  -  0,  not  merely  approximately  identical  perfor¬ 
mance.  On  the  other  hand,  this  uniform  improvement  is  for  the  case  of  II 
supported  by  the  axes.  For  the  cases  off  the  axes  for  cubes  and  for  the  case 
of  disks  there  is  a  trade-off  between  some  loss  away  from  zero  and  the  gain 
near  zero. 


APPENDIX:  REVIEW  OF  DISTRIBUTIONS  OF  COMPACT  SUPPORT 


We  present  here  a  short  review  of  some  properties  of  distributions  of 
compact  support  in  Fn  (sometimes  called  generalized  Functions) .  The  impli¬ 
cit  reference  throughout  this  Appendix  is  the  very  clear  monograph  [22]. 

The  simplest  example  of  a  distribution  of  compact  support  is  obtained 
from  a  continuous  function  F  of  compact  support.  (That  is,  F  is  zero  out¬ 
side  a  bounded  set  in  ff?n. )  One  associates  to  F  a  scalar  valued  linear 
operator  Tp  acting  on  infinitely  differentiable  functions  f  by  the  rule 


<Tp,f> 


C 


F(x)f (x)dx. 


The  fact  that  supp  F  is  bounded  is  what  makes  the  above  integral  finite  and 
well  defined  for  an  arbitrary  f.  This  linear  operator  Tp  is  called  a  dis¬ 
tribution  of  compact  support. 

If  the  function  F  were  not  only  continuous  but  also  continuously  dif¬ 
ferentiable,  we  can  similarly  associate  a  distribution  to  each  of  its  partial 
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derivatives.  Their  relationship  to  Tp  is  given  by  integration  by  parts: 


<T 


a  f 
ax 


•f>  - 1  *  -I 

JRn  J  JRn  J  J 


g 

This  motivates  defining  the  derivative  Tp  of  the  distribution  Tp  by 

J 

the  formula 

<a^W>  :■  -<TF.%>. 

which  does  not  require  that  F  be  differentiable.  One  defines  in  the  same 

al«i 

way  more  distributions  D  T_  =  - T_,  a  =  (a . a  ),  |a|  = 

F  axf1. . .ax“n  F  1 


n 


a  ♦  .  .  .  ♦  a  ,  by 


<D®Tp.f>  :=  (-l)|a,<TF.Daf>. 


Further,  for  any  pair  of  continuous  functions  F.G  of  compact  support, 
and  any  multi  indices  a,  a,  the  operator  D0^-  +  D^T 

r  (j 

<d“tf  +  D^TG,f>  :=  <D“TF,f>  +  <D^Tc,f> 

is  a  linear  operator  on  the  space  of  infinitely  differentiable  functions. 

The  collection  of  all  finite  sums  of  this  kind  is  precisely  the  space  of 
distributions  of  compact  support  in  Rn.  That  is,  for  a  distribution  of  com¬ 
pact  support  p  we  mean  a  scalar  valued  linear  operator  which  can  be  repre¬ 
sented  in  the  form 

m  -  £  D V 

lalSN  “ 

for  some  integer  N  and  continuous  functions  Fft  of  compact  support  ((22), 
p.  741). 
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Let  us  Indicate  how  the  unit  impulse  at  a  point  a  c  F  can  be 

written  in  the  above  for*.  Let 


L(x) 


|x~E 


if  x  <  a 
if  x  i.  a 


and  let  $  be  any  twice  continuously  differentiable  function  of  compact 
support  such  that  for  some  c  >  0  ^(x)  *  1  for  a-c  £  x  £  a  +  c.  Then 

(0L)"  =  0"L  +  20'L'  +  ^L",  and  the  function  $'L'  is  continuous  in  spite  of 

the  discontinuity  of  L'  at  x  =  a  since  s  0  near  a.  Letting  F  =  <pL 

and  G  =  -(^''L+20'L')  we  get 


6 

a 


dx 


— T  +  T 
27  C‘ 


The  space  of  infinitely  differentiable  functions  in  Rn  is  usually 

denoted  by  &  (or  £(Rn )  or  C®(Fn))  and  the  space  of  distributions  of 

compact  support  by  &'  (or  £'(Rn)). 

We  can  essentially  work  with  distributions  of  compact  support  as  if  they 

were  ordinary  continuous  functions  of  compact  support.  In  particular,  for. a 

distribution  T  €  6'  it  makes  sense  to  compute  its  Fourier-Laplace  transform 

T(<;)  for  <  €  Cn.  Namely,  let  £*x  =  £,x,  ♦...+?  x  (x  €  IRn,  <  €  Cn) , 

11  n  n 

then 


T(<)  :=  <T,e 


where  i  =  V-\.  For  instance,  for  the  unit  impulse  at  a  point  a  e  Rn 

we  have  <5  (C)  =  e  ^  a. 
a 

From  the  above  representation  of  the  distributions  of  compact  support  one 
can  see  that  the  function  *T(<[)  has  two  properties.  First,  it  is  an 
analytic  function  in  the  whole  of  Cn  (i.e.,  am  entire  function).  Second, 
there  are  positive  constants  A, B, N  such  that  T  satisfies  everywhere  the 
estimate 
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|T(OI  S  A(1^ICI)M*B,I"?I. 

Im  <  «  (I®  . I®  Cn)-  It  is  usual  to  call  the  space  of  functions  satisfy¬ 

ing  these  two  properties  the  Paley-Wlener  class,  PW(Cn). 

The  classical  Paley-Wiener  theorem  for  functions  of  compact  support  in 

p 

L  cam  be  extended  to  the  case  of  distributions  of  compact  support  as 
follows. 

Theorem  (Paley-Wiener-Schwartz) .  The  Fourier-Laplace  transform  is  a 
one-to-one  correspondence  between  the  spaces  £'(Rn)  and  PW(Cn). 

Recall  now  that  if  F  and  G  are  continuous  functions  of  compact 
support  then  their  convolution  F*G  is  given  by 

F.G(x)  =  |  F(x-y)G(y)dy  =  |c(x-y)F(y)dy 

R 

and  it  is  again  a  continuous  function  of  compact  support.  Moreover,  if  one  of 

Q 

them  is  continuously  differentiable,  say  F,  then  ^ — (F»G)  = 

J 

observation  allows  for  the  definition  of  convolution  of  distributions  of  com¬ 
pact  support  in  such  a  way  that  if  T,S  €  S'  then  T*S  =  S»T  e  S'.  The  unit 
impulse  at  the  origin  6  acts  as  the  unit  of  this  product  T»5  =  T.  Further¬ 
more,  (T*S)  (O  2  T(OS(C).  so  that  the  convolution  becomes  ordinary  pro¬ 
duct  of  the  analytic  functions  T  and  S. 

In  this  paper  we  were  interested  in  solving  the  equation 

Vi'-'V".  *  *• 


where  p, . p  are  given  in  S'  and  we  need  p,,...,v  eg'.  Using  the 

Id  l  m 

Fourier-Laplace  transform  this  equation  is  equivalent  to  the  analytic  Bezout 
equation 
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MjtOfjtO  ♦...  ♦iB(Ofli(C)  ■  i. 

with  vinknowns  f f  €  PW(Cn).  In  particular,  to  solve  this  equation  is 

1  to 

a  e  n 

necessary  that  the  functions  p. . p  have  no  common  zeros  in  C  . 

X  HI 

In  the  text  we  raised  the  question  of  whether  one  could  find  a  deconvol¬ 
ver  v  €  S'  for  a  single  convolutor  p  €  S' ,  i.e., 

A  A 

p«i>  =  6  or  equivalently  p(<;)i/(£)  =  1. 

We  need  that  p(£)  *  0  for  every  £  €  Cn.  It  follows  that  there  is  an  analy¬ 
tic  function  h  in  Cn  such  that 


p(C)  =  e 


h(0 


Since  p  is  in  the  Paley-Wiener  class  one  can  show  without  difficulty  that 


|h(Ol  *  A  +  Bllm  O 


for  some  positive  constants  A,B.  Therefore,  by  the  Liouville  theorem,  h 
must  be  of  the  form 


h(?)  =  -i( >  . . .  ♦  anCn)  +  c 


for  some  a  =  (a,,..., a  )  e  P  ,  c  e  C.  This  proves,  by  the  the  Paley-Wiener- 
l  n 

Schwartz  theorem,  that  p  =  C5  ,  for  some  constant  C  €  C. 

& 

Having  no  common  zeros  is  not  enough  to  solve  the  analytic  Bezout  equa¬ 
tion,  the  necessary  and  sufficient  condition  is  the  following. 

Theorem  [1],  [16].  Given  p,,...,p  €  6'(Fn),  the  necessary  and  sufficient 

1  m 

condition  for  the  existence  of  deconvolvers  y,,...,u  €  S' (Rn)  is  that  there 

1  HI 

are  constants  e>0,  M>0,  C>0  such  that 


1/2 


-C|Im  O 


i  c 


(1*10) 


M 


for  all  <  €  C 
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This  Is  a  purely  exlstenclal  theorem  and  the  work  in  (5}  -  (10)  consists 
In  obtaining  explicit  choices  of  deconvolvers  .  This  is  not 

trivial.  A  portion  of  the  work  (7)  and  111)  also  consists  in  finding  simple 

and  practical  examples  of  p, . p  that  satisfy  the  above  necessary  and 

1  m 

sufficient  condition  for  the  existence  of  deconvolvers. 
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Figure  1.  (a)  Multiple  parallel  linear  operators  with  distinct  distributions 

Single  operator  (b)  and  the  multiple  parallel  operators  with  identical 
distributions  (c). 
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CHAPTER  3 


A  local  version  of  the  two-circles  theorem 


by 


C.A.  Berenstein  and  R.  Gay 
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Introduction 


One  of  the  oldest  questions  in  integral  geometry  has  been 

that  of  recovering  a  function  f  in  Rn  from  the  knowledge  of 

its  average  over  balls.  It  is  easy  to  see  that  unless  f  decays 

sufficiently  fast  at  infinity  the  average  over  all  balls  of  a 

fixed  radius  could  vanish  without  f  being  identically  zero.  It 

is  not  always  possible  to  assume  such  decay  but  a  very  elegant 

result  of  Zalcman  [20]  and,  independently,  Br own-Schreiber-Taylor 

[10],  describes  explicitly  a  countable  set  such  that 

averages  over  all  balls  of  radii  r,,  r2  suffice  as  long  as 

r,/r„  t  E  .  Th is-" two  circles"  theorem  can  be  described  as 
1  '  2  r  n 

saying  that  the  map 


C<Rn)  ♦  C<Rn)  O  C(Rn) 

f  *  {  !  f ( y ) dy ,  /  f(y  )dy) 

B  (  x  ,  r  )  B  (  x  ,  r  2  ) 

is  injective  if  and  only  if  r  /r  ^  . 

(B(x,r)={y  Rn  :  |  x-y  |  < r } )  .  Under  slightly  stronger  conditions  on 
the  quotient  r 1 /r  this  map  has  also  a  continuous  and  explicit 
inverse  [8].  This  result  and  other  variants  ot  the  so-called 

iompeiu  problem  have  bean  generalized  to  symmetric  spaces  'see 
the  surveys  [ 2 l ] ,  [1]  for  positive  results  and  their  limitations) 

In  practical  situations  of  a  tomographic  nature  one  is 
limited  to  balls  that  fit  into  a  fixed  region  ft.  One  could  take 
smaller  and  smaller  balls  when  approaching  the  boundary  3ft  of 
ft,  this  is  roughly  the  situation  when  we  consider  the  case 


Q  ■  unit  ball  of  Hn  ii  the  hyperbolic  apace,  but  it  ia  clear 
that  it  Bight  be  hard  to  accoapliah  if  we  are  dealing  with 
physical  devices  whose  size  cannot  be  made  infinitesimally  small 
or  cannot  even  be  changed  at  will.  It  is  this  kind  of  problem 
that  we  call  a  local  version  of  the  two— circles  theorem.  The 
main  difference  with  the  above  mentioned  results  is  that  we  do 
not  have  any  longer  the  whole  group  of  Euclidean  motions  at  our 
disposal  which  was  the  crucial  ingredient  lying  behind  the  two 
circles  theorem  and  its  generalizations.  The  Inversion  formula 
of  18]  would  allow  us  to  reconstruct  f  away  from  but  gives 

no  indication  of  whether  we  could  change  the  values  of  f  in  a 
collar-like  region  near  without  affecting  its  average. 

There  is  some  recent  work  on  systems  of  convolution  equations  in 
convex  domains  which  deals  with  this  type  of  question  [4]  but  the 
hypotheses  required  are  far  too  restrictive  to  be  satisfied  by 
our  simple  looking  problems.  Nevertheless,  using  a  combination 
of  ideas  from  classical  harmonic  analysis  and  results  of  Cormack- 
Quinto  on  the  Radon  transform  on  spheres  [12]  we  are  able  to 
prove  the  following. 

Theorem .  -  Let  r 1 , r  be  positive  numbers,  r1/r2  ^  En ' 

Q  an  open  subset  of  Rn  such  that  every  point  lies  in  an  open 
ball  contained  in  (2  of  radius  strictly  larger  than  r  ^  +  r 

If  f  €  C(fl)  satisfies 

/  f(y)dy  «  0  for  every  B  (  x  ,  r  .  )  £  fl ,  j  *  1,2 

B  (  x  ,  r  )  3 

j 
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then  £  =  0.  Furthermore,  this  etetment  does  not  hold  if  D 
fails  the  above  geometrical  restriction. 

The  method  of  proof  allows  us  to  generalize  this  theorem 
greatly,  providing  in  particular  new  local  mean-value  theorems 
for  harmonic  functions. 

We  will  like  to  express  our  appreciation  to  Professor  L. 
Zalcman  who  called  our  attention  to  these  problems. 

The  second  author  wishes  to  thank  the  Mathematics  Department 
of  the  University  of  Maryland  for  its  hospitality  while  this  work 
was  carried  out. 

2 .  Preliminaries.  - 

We  will  follow  the  standard  notation  for  distributions  found 
in  [14).  We  denote  B(x,r)  =  { y €r" : | x-y | < r } ( r > 0 ) ,  B(x,r)  its 
closure  and  Xr  the  characteristic  function  of  B(0,r).  Let  ft 


be  an  open  set  in  Rn , 

fl  = 
r 

{  xC  n  :  d  (  x  , 

nC)<r}  . 

For  a  locally 

- 

integrable  function  f  in  an 

open  set 

fi  the 

average 

( 1  )  Af (f , x) 

_L_  / 

f  (  y  )  dy 

u 

n 

n  _  ,  , 

r  B  (  x  ,  r  ) 

is  defined  for  x  £ 

He  re 

(i)  is  the  volume 
n 

of  B  (  0  ,  1  )  . 

If 

we  let  \i  =  X  /<*>  t  ,  we 

r  r  n 

can 

interpret 

this  average  as  a 

convolution  and  hence  it 

makes 

sense  to 

define 

it  for  f  £  V' 

(n ) 

giving  a  distribution  Ar 

if) 

in  V  (fi 

r 

) ,  name  1 y 

A  (f)  *  f  *  v  .  Therefore,  for  uniqueness  questions,  if  the 
r  r 

averages  of  f  are  zero,  by  restriction  ourselves  to  0^  ,  c  >  0 

m 

small,  we  can  assume  f  £  C  .  Henceforth,  all  distributions  with 

vanishing  averages  will  be  assumed  to  be  C“  functions  in  fl. 


the  distributions  defining 


For  r  >  0,  we  denote  by  o 

r 

the  spherical  average 

(2)  A  (f ,x)  ■  /  ,  f(x+ry)dc(y)  -  (o  *f)(x) 

r  Sn_1  r 

do  is  the  normalized  Lebesgue  measure  on  Sn_1. 

* 

For  T  €  F  the  Fourier  transform 

T(t)  «=  <  T  ,e  (x  |  C  )  *  l  x  .  C  •  » 

x  j  3  3 

is  an  entire  function  in  Cn  which  satisfies,  for  some 
A,  N  >  0,  the  estimates 

(3)  |TU)|  <  A(l+|t|)N  exp  (  H  ( Im  C )  )  . 

where  £  =  £  +  in ,  £,  n  €  Rn  ,  Im?  =  n  and  H  is  the 

supporting  function  of  the  support  of  T,  i.e.: 

H(n)  =  Ma x{  ( x  |  n )  : Xf suppT) 

Note  that  H  is  also  the  supporting  function  of  cv(suppT),  the 
convex  hull  of  suppT.  The  Fourier  transform  is  an  isomorphism 

between  the  convolution  algebra  F  (R  )  and  F  (Rn),  the 

algebra  of  entire  functions  of  exponential  type  and  polynomial 
growth  on  the  real  axis. 

A 

A  distribution  T  will  be  called  invertible  (or  T  is 

•  n  “ 

slowing  decreasing)  if  whenever  S  (  F  (R  )  and  S/T  is  an 
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9  ah 

entire  function,  then  there  is  •  distribution  U  9  f  (» )  such 

*  A  A 

that  0  -  S/T,  that  is 

(4)  S  -  T  *  U 

and  we  have  the  identity 

C  5 )  Hs  =  HT  +  Hy 

or,  what  amounts  to  the  same  thing 

(6)  cv( suppU*T)  =  cv(suppU)  +  cv(suppT) , 

where,  for  two  sets  A,B  £_  Rn  we  have  A  ±  B  =  {  x±y ;  x  6  A,  y  €  b}  . 
We  will  need  to  use  that  Wr  is  an  invertible  distribution. 

A 

This  will  follow  from  the  explicit  formula  for  given  below 

and  the  characterization  of  invertible  distributions:  T  is 

invertible  if  and  only  if  there  is  a  positive  constant  a  such 
that  for  all  £  £  Rn 

(7)  Max  {  |t  (  £+n  )|:n  €  Rn  ,  |  h  |  <a  .  log  (  2+ |  £  |  )  }  >  (  a  +  |  £  |  )  3 

The  Fourier  transform  of  a  radial  distribution  T  is  a 
radial  function,  i.e.:  if: 

<T,foA  ^  >  “  <T,f> 
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for  «vary  A  (  0(n)  then 


T(C)  ■  T(A.;) 

for  every  A  €  0(n),  £  f  ,  and  depends,  for  £  €  R° ,  only 

on  | £ | .  Hence  we  consider  the  associated  even  entire  function 

#*•» 

T  of  one  variable  by 

(8)  X  (  |  £  |  )  *  T  {  £  )  and  T  (  C )  -  ?(  (  £  2  + .  .  .  +  £2 

Let  us  call  E  ( Rn ) ,  the  space  of  radial  distributions.  This 

o 

correspondence  establishes  an  isomorphism  between  the  algebras 

E ' ( R° )  and  E  (R).  Using  this  notation  we  have 
o  o 


n-2 


(9)  Hr(t)  = 

n  2  2 

n|) 

Jn  (rt)/(rt)n/2 

2 

n-2 

n-2 

(io)  o ^ ( t ) 

2  2 

r<|> 

n-2 

( r  t )  /  ( r  t )  2 

2 

and,  more  generally,  if 

f  (x)  = 

♦  (lx 

1) 

is  a  radial  function  of 

compact  support 

(11)  f  (  £)  =  f (t)  = 

<2,)n/2 

n-2 

r 
;  0 

♦  <  P  ) 

Jn-2  (  Pt)pn/2dp  (  |  £  |-t) 

t  2 

2 

To  show  that  is  invertible  it  is  now  sufficient  to 

recall  the  asymptotic  development  of  the  Bessel  functions  [19]  on 
the  positive  real  axis 
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(12) 


Jv(t) 


2  tJ/2 


co.lt ♦  0(t'3/2) 


It  follows,  for  |  C  ]  >  1  and  some  C  >  0 


n+  1 


Max{|wr(^+n)|:ri€Rn,|n|<Tf}  >  c|t| 


which  is  the  condition  of  i n ve r t ib i 1 i ty  . 

From  (12)  we  also  obtain  MacMahon ' s  asymptotic  development 

of  the  positive  zeros  a,  of  J 

k  ,  v  v 


0  <  a 


,v 


<  a 


2 ,  v 


<  .  .  . 


(13)  ak  v  =  (  2k  + 1  )•—  +  (2v  +  i)-J  +  0(1  /k  ) 


which  will  be  used  further  on. 


3 .  Series  development  of  mean-periodic  function s  . 
Let  SI  be  an  open  convex  set  in  Rn  and 
K  =  cv  ( suppjj  )  ,  \i  €  £  (Rn).  We  say  that  a  function 

oo 

f  (.  C  (SI- K)  is  mean-periodic  with  respect  to  u  if 


(14)  p*f(x)  =  <y  ,f(x-y)>  *  0  for  all  x  €  Si 

If  an  exponential- polynomial,  that  is  a  finite  linear 
combination  of  terms  of  the  form  x^ 


a  j,  in 

(xJ»x4  ...x  »  j.  €  K,  1  «  k  <  n),  i«  Bean-periodic  with  reapect 

i  n  x 

to  V  then  the  frequencies  C  Bust  satisfy 

* 

p(<;)  •  0  since 


(15) 


( V*e 


(•jO 


)  (x) 


V  (  Oe‘ 


(x  |  O 


when  the  zeros  of  p  are  simple  no  non  constant  monomials 
can  appear.  More  generally  if  a  monomial  x3  appears  with  non 
zero  coefficient  then 


,g  I  3  I  - 

(^-)u(0  =  0 

a  c3 

for  the  corresponding  frequency  C. 

For  n  =  1  there  is  a  well-known  series  development  for 
such  functions  in  terms  of  the  exponential  polynomial  solution  of 
the  same  convolution  equation  (14)  due  to  L.  Schwartz  [18],  [15], 

[13].  The  case  of  interest  for  us  is  n  >  2,  p  invertible.  In 
this  case,  a  development  in  terms  of  integrals  over  the  zero  set 

A  ^ 

of  p  has  been  proved  when  fl  =  R  [6].  For  ft  arbitrary 
convex  set,  a  similar  development  has  been  proved  in  [4]  but  only 
for  a  very  restrictive  class  of  invertible  distributions.  In  all 
these  cases  one  obtains  also  some  knowledge  of  the  behavior  of 
the  terms  involved  in  this  development.  Unfortunately,  the 
distributions  p^,  though  invertible,  do  not  satisfy  the 
conditions  required  in  [4],  as  was  shown  (for  a  different 
reason)  in  [3] ;  moreover  we  are  interested  in  Q  =  B(0,R) . 
Therefore  we  cannot  depend  on  any  of  the  previously  known 
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results.  We  obtain  here  a  aeries  development  without  additional 
information  on  the  coefficients  that  appear  in  it;  nevertheless 
the  existence  of  this  development  is  all  we  need  later. 

Propos i ti on  1  -  Let  fl  be  an  open  convex  subset  in 
Rn  (n>2),  y  £  E  (Rn)  an  invertible  distribution, 

m 

K  =  cv(suppy).  Any  function  f  (  C  (fi-K),  mean  periodic  with 
respect  to  y  can  be  written  as 


(16) 


f ( x )  =  I  P.(x) 

j>1  3 


(xcn-K ) 


with  Pj  exponential-polynomials  also  mean-periodic  with  respect 

OD 

to  y,  and  the  series  is  convergent  in  the  C  -topology  of 
fi-K.  Furthermore,  given  a  sequence  of  Positive 

numbers,  letting  PQ  *  0,  we  can  chose  the  P..  so  that  the  absolute 

value  of  all  frequencies  in  Pj+i  exceeds  the  largest  absolute 
value  of  the  frequencies  in  Pj  by  at  least  s^+^. 

Proof.  Let  us  show  first  that,  for  any  s  >  0,  the  exponential 
polynomials  which  are  mean-periodic  with  respect  to  y  and  whose 
frequencies  lie  outside  the  ball  of  center  0  and  radius  s 
in  tn  are  dense  in  the  space  N  =  {f€C  (fi-K):y*f-0  in  fi}.  N 
is  a  closed  subspace  of  a  Frechet  space  and  we  only  need  to  show 

f 

that  if  v  £  £  (fi-K)  is  orthogonal  to  the  above  exponential- 

A  N* 

polynomials  then  v  is  orthogonal  to  N.  Hence  (v  )  is 

A  ^  ___ 

divisible  by  y  at  every  point  of  C  \  B(0,s).  Since  n  >  2, 

A  v  A 

by  Hartogs '  theorem,  (v)  /y  is  an  entire  function.  Since  y 
is  invertible  there  is  a  distribution  T  (  E  ( Rn )  such  that 
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V  -  U  *  T 


We  need  to  know  where  is  the  support  of  T.  By  (6) 


cv(suppv)  e  cv(suppp)  +  CV(SUppT) 


or 


cv(suppT)  -  K  =  cv(suppv)  £  ft  -  K 

By  the  Hahn  -  Banach  theorem  one  concludes  that 

cv(suppT)  £  ft 

v 

Hence  <v,f>  =  (v*f)(0)  =  (T*y*f)(0)  *  <T,y*f>  =  0  for  f  €  N. 

To  end  Lhe  proof  uf  the  proposition,  v®  pick  an  exhaustion 
of  ft  -  K  by  convex  compacts  sets  K j ,  hence  we  can  find  P,  , 
exponential-polynomial  with  frequencies  lying  in 

A 

{  c€<rn:u  (  C)  =  0,  |  5  1  >s1  }  such  that 

sup | f -P  |  <  1 . 

K1 

Let  =  maximum  of  the  absolute  values  of  frequencies  in 

P1  .  We  can  find  Pj  with  frequencies  in 
|cf«n:yU)  -  0,  |c|  >  *2+°,}  * uch  that 

■ax  eup|oa(f-P, -P_ ) |  <  1/2 

|.|<i  k2  ’  2 
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Continuing  in  this  fashion  we  obtain  the  desired  expansion.  □ 


Rema  rk .  One  can  eliminate  the  requirement  of  p  being 
invertible  by  using  [14,16.4.1}. 

From  (9)  we  know  that  the  zero  variety  of  is  the  union 

of  the  hyper surf aces 

1.7)  t2  =  »  *  >.2.... 


where  X^  =  2/r.  We  disre<3ard  temporarily  the  dependence 

on  r  though  it  will  play  a  role  later  on.  Furthermore  the 
function  pf(t)  vanishes  at  t  =  X^  with  multiplicity  one,  in 
fact 

n-2 

(18)  ft",*4’  *  -"2  2  r(f’r  J(„/2>.,'rt»''‘">n/2 

and  well  known  properties  of  Bessel  functions  show  that  this 
expression  does  not  vanish  for  t  =  X^.  Using  the  asymptotic 
expressions  (12)  and  (13)  we  obtain 


(19) 

du 

0  +  IT  <V 


n  + 1 


n  +  3 


r  n 


2 ( n -  1  ) / 2  r(n/2) (_i)k  +  1/(A  r)2 


o(fc  2  ) 


We  introduce  some  auxiliary  radial  distributions 
formula 


T 


r  ,k 


by  the 


(20) 


Wr(t) 


k 
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Hence  they 


They  are  even  and  entire  since  u^fiX^)  *  0. 
correspond  to  radial  distributions  (in  fact  C1 ' 1  functions)  whose 
supports  are  contained  in  the  support  of  i.e.  B(0,r). 

Furthermore  they  satisfy 


(  21  ) 


(4  *  Xk)Tr,: 


=  - p  and 

r 


(  22  ) 


T  .  (  X  ) 
r  ,  k  k 


"r(V 

2X. 


k  + 1 


n  +  3 


=  const. (-1) 


n  +  5 
2 


0(k  ). 


We  remark  that  these  distributions  have  conspicuously 
appeared  in  previous  work  on  the  Pompeiu  problem  [2],  (7). 

Proposition  2  -  Let  r  >  0  be  fixed.  For  any  p  ,  0  <  p  <  °» ,  we 

can  decompose  in  the  following  form 


(  23  ) 


0 

P 


v  +  y 
P  r 


S 

P 


where  S  is  a  radial  distribution,  whose  support  satisfies 
P 


(  24  ) 


supp 


S  c 
P  - 


b( 0 , Ma x  (  r 


P  )-r) 


and  v  is  given  explicitly  by: 


(25) 


v 

P 


l 

k>  1 


°a(Xk) 

_JL_* -  AT 

X 2  T  ( X  )  r,k 
k  r ,  k ' Ak  ' 


hence  supp  v  c  B(o,r). 
P  ~ 
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Proof 


At  consider  the  series 


(26) 


9(t>  -  I  -  t2  T  k(t) 

k>’  \Tr.*'V 


The  coefficients  X  a  (X)/T  .  (X  )  are  uniformly  bounded 

x  p  x  r ,  x  X 

by  a  constant  depending  only  on  p  as  it  can  be  seen  from  (10), 

(135  and  (22),  since  X^  ~  const,  k.  Therefore,  if 

| t |  <  R  n  >  2R  we  have  |t  T  .(t)|  <  const,  k  which 

x  r  i  k 

guarantees  the  convergence  of  the  series,  and  shows  g  is  an 
even  entire  function.  We  can  obtain  more  precise  estimates  by 
picking  a  sequence  of  circles  of  center  0  and  radii 


R_.  =  (4j+n+5)  /4r,  j  =  1,2,.. 


Decomposing  the  sum  into  those  terms  where 


X,  <  2R ,  and  X,  >  2R .  one  can  estimate  the  second 
*  J  *  3 


sum  over 


I  t  J  =  Rj  by 


Max  I t^  y  ( t ) | . 
|t|=Rj 


c  (p) 

O 


The  first  (finite)  sum  can  be  estimated  by 


C  (p)(Max|  t2y  (t)  |)  (Max  J  - -  ) 

®  I  i.  I  .  *  r  r\  -  \  M  I  \  ^  *  I 


f.jit0<Xk<2Rj|i2-t2 


where  Q  is  the  region  obtained  from  |t|  <  R  by  removing 

j  » ®  3 

disks  of  radius  c,  0  <  c  very  small,  about  *X^.  One  can 
then  see,  without  difficulty,  that  the  last  sum  is  estimated  by 


const.  c 


In  any  case  we  obtain  as  a  final  estimate 


Max  |g(t)|  <  C ( p )  Max  |t2u  (t)|. 

|t|<R. 

Thus  g  defines  a  radial  distribution  of  order  2, 

Vp,  by  v  ^  *  g,  one  can  see  is  given  explicitly  by  (25). 

We  also  have 


with  h  even  entire  function  since  g(±X  )  =  v  (±X  )  *  o  (±X  ) 

K  p  K  p  K 

by  (26).  Since  p  is  an  invertible  distribtion  it  follows 
J  r 

h  =  S  for  some  S  €  E  (R  ).  The  identity  (6)  gives 
P  p  o 


(27)  cv(supp(0p-vp ) )  =  cv(suppSp)  +  cvtsuppy^) . 

There  are  two  cases  to  consider.  If  p  <  r,  then  the  support  on 

the  left  hand  side  of  (27)  is  contained  in  B(o,r)  and 

cvfsupp  S  )  =  | 0 } ,  which  says  S  is  a  polynomial  in  the 

P  P 

Laplace  operator;  if  p  >  r  then  the  left  hand  side  of  (27)  is 
contained  in  B(o,p),  which  says  cv(suppSp)  £.  B(o,p-r).  □ 


Rema rk  The  decomposition  we  have  just  given  in  proposition  2 
works  also  if  we  replace  by  any  radial  distribution.  We 

need  onl\  j  change  ( t/xJc )  by  (t/X^)  q  with  q  convenient 
non  uege  ,  a  integer.  In  particular  there  is  such  a 
decomposit.  with  oQ  -  6,  the  Dirac  ease  at  the  origin  (take 

%  "±1% 

i  >  — > 
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Corollary  3  -  Let  f  be  a  y^  -mean-periodic  function  in 
C  (b(0,R))  ( R>r ) .  Let  |xQ|  <  R  -  r.  Then,  for  any  p, 

0  <  p  <  R  -  IxqI  we  have 


o  (X  ) 

(28)  X  (f,xn)  «=  (v  *f)(0)  =  -  l  -~-z -  A(T  *f)(x  ) 

”  0  p  k>1  X*  ?  ( X  )  r'k  0 

k  r ,  k  k 

Proof.  It  suffices  to  use  (2)  and  (23).  D 


4 .  Local  two-circles  theorem. 

Let  r ^ ,  rj  be  two  positive  numbers  and  consider  the 

distributions  y  ,  y  They  will  have  no  common,  mean- 

r  1  r  2  » 

periodic,  exponential-polynomials  if  and  only  if  y  and  y 

ri  r2 

have  no  common  zeros.  By  (17)  this  occurs  if  and  only  if 
r^/r2  ^  quotient  of  two  zeros  of  J  2 

The  set 


E 

n 


‘k  ,  n/2^a  j  ,  n/2 


:1< j ,k<«} 


is  the  exceptional  set  described  in  the  two-circles  theorem. 
Proposition  4  Let  R  >  r ^  +  r2»  r1/r2  $  En •  The  only 

m 

function  in  C  (B(0,R))  which  is  mean-periodic  with  respect  to 

both  y  and  y  is  the  zero  function. 

r- 

Proof .  We  assume  <  r2»  Let  f  ^  C  (B(0,R))  be  y^ 

mean  periodic.  By  proposition  1  we  have 
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lie 


f (x)  -  l  P .  (x)  ( |x |  <R) 

j>!  3 

where  the  frequencies  appearing  in  the  exponential  sums  P^ 
i  n 


{c*Cn:U  (0  =  0}  =  U  {  Ocn  :  ?2=(a  n/2/r )  2)  “  U  V  . 

r.  ^  s i  k>1^ 


k>  1 


We  fix  now  k  >  1,  and  consider 
C  (B(0#R-r  )),  furthermore 


T 

r 


which  is  in 


(29)  TK  k*f  =  I  Tr  v*p< 

rl'k  j>1  ri'  3 

If  P.(x)  =  l  C.  el(xlCj,£)  then 

J  £  3  '  * 


*P  . 
3 


I 

£ 


'j  »  * 


rr,,k(S. 


i(x  I  C. 


but  T  (  C .  4  0  only  if  €  V  in  which  case  we  obtain 

r,  i  i  *  * 

the  value  T  (X  )  *  0  (where  X  is  computed  with  respect 
r  ^  f  K  K  K 

to  r.j).  Therefore 


(30) 


T  k*f  *  k(Xk}  I  pn  k 

ri'*  ri 'K  *  j>1  3 '* 


where  P,  . 

j 

lie  in  . 
convolve  now 


is  the  sum  of  the  terms  in  P^  whose  frequencies 
This  series  is  convergent  in  C  (B(0,R-r  ))■  We 
with  y  *  We  obtain 


(31  ) 


■r  «V 

2 


j>1 
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since  y^  is  also  a  radial  distribution.  The  expansion  (31)  is 
valid  in  C  (  B (  0  , R-r  1  -r  2  ) )  .  Since  f  is  also  y^  -  mean  - 
periodic  we  have 


0  »  (t  *y  *f)(x)  =  y  (X  )(t  *f)(x) 

r  1  >K  r  2  T  2  K  i'K 


for  | x |  <  R  -  r 1  -  r2«  The  hypothesis  r^/r2  4  En  n°w  implies 

that  y  (A  )  ^  0.  Hence 
r  2  * 


(32)  (t^  k*f)(x)  =  0  for  1 x  |  <  R  -  r^  -  r2 


On  the  other  hand  we  have  (by  (22)) 


(A+A^)(Tr  ,k*f)  =  -(f*Ur  )  =  0  in  | x |  <  R  -  ^ 


hence  T  *f  is  a  real  analytic  function  in  |x|  <  R  -  r  .  We 

E  y  §  K  » 

conclude  that 


(33)  (T  *f ) (x )  *  0  for  | x [  <  R  -  r 

r  ^  #  k  l 

Applying  now  corollary  3,  formula  (28),  we  have 

(34)  A  (f,x)  -  0  whenever  |x|  <  R  -  r^  ,  0  <  p  <  R  -  |x| 

(We  are  allowed  to  take  p  «  0  by  continuity).  In  particular 
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£ ( x )  ■  0  for  | x |  <  R  -  r 

To  finish  the  proof  of  the  proposition  we  need  to  show  f 
is  zero  in  the  remaining  annulus,  we  do  that  using  (34).  It  is 
at  this  point  that  we  use  Cormack-Quinto  (12).  For  any  y 
B ( 0 , R ) ,  consider  R(f)(y)  «  X j y j ^ ( f , y/2 ) .  This  is  the  Radon 

transform  on  spheres  through  the  origin  discussed  in  [12).  We 
want  to  show  Rf(y)  =  0.  We  only  need  to  verify  that  the 
conditions  stated  in  (34)  are  valid.  Here 
p  =  | y 1/2,  x  =  y/2,  hence 

R  -  |  x  |  ■  R  -  *  R  -  p  >  R/2  >  p 


The  only  condition  left  to  see  is  that  |x|  <  R  -  r  .  We 

have  2r^  <  r1  +  <  R  hence  r^  <  R/2  and  R  -  r1  >  R/2, 

therefore  |x|  <  R  —  r ^  holds. 

By  [12,  corollary  2]  f(y)  «=  0.  (We  note  that  in  [12],  they 

®  n  ® 

require  that  f  C  (R  )  while  we  only  have  f  (  C  (B(0,R))  but 

the  proof  of  corollary  2  depends  on  an  explicit  inversion  formula 

for  the  Radon  transform  on  spheres  which  uses,  for  each  y,  values 

of  f  in  a  neighborhood  of  B(0,|y|).)  □ 

Remark .  The  crucial  point  of  the  proof  above  is  (32).  One 

does  not  really  need  the  whole  strength  of  Propositon  1  to  obtain 

it.  One  can  get  by  using  the  density  of  the  exponential 

polynomial  solutions  in  the  sub-space  N  introduced  in 

Propositon  1.  Nevertheless,  we  feel  that  the  proof  is  clearer  using 

the  expansion  (16)  as  we  have  done. 
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tie  want  to  show  that  the  condition  R  >  r1  ♦  r2  is  sharp. 
It  is  easier  to  show  this  under  the  slight  restriction  that 
r2/r1  is  not  too  well  approximated  by  elements  in  En . 

Definition.  -  For  N  >  0,  we  say  that  a  positive  number  is 
N-well  approximated  by  points  in  En  if,  for  every  l  >  1, 
there  are  indices  j,k  such  that 

(35)  I r  -  «k/aj|  < 


where  a,  =  a, 

k  k,n/2 

Proposition  5  For  any  N  >  2, 
approximates  by  En  has  zero  measure 
Proof .  Given  p,q,  0  <  p  <  r 
we  have 


the  set  of  numbers  N-well 
in  (0,») 

<  q  and  v  >  0,  from  (13) 


»  -  (2k+l)v/2  +  (2v+l)*/4  +  0(1/k) 


Therefore,  if  r  satisfies  (35),  for  l  >  1,  we  have 


(36) 


|r.j  -k+Ar+B|  <  C 


for  some  constants  A,B,C.  Hence 

pj  -  C1  <  k  <  q  j  ♦  C2 

for  some  constant  C1 ,  C2  >  0.  Hence  the  cardinal  of  the  set 
of  k  satisfying  (36)  is  bounded  by  (q-p)j  ♦  L,  L  constant  >  0. 
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Now ,  the  set  of  N-well  approximated  numbers  in  (p,q)  is 


(37)  n  U  {r:  P<r<q,  Jr-a^a  |  <  1/*jN} 

1>1  j  ,  k>  1  3 


For  1  fixed,  the  Lebesgue  measure 

|  U  {  r  :  p<  r  <  q  ,  |r-a  /a  |  <  1/S.jN  }  |  <  j  l  ( q  -P^3-~L  <  ~j 
j  * k  J  j>l  j  ' 

(C3  >  0)  since  N  >  2.  Therefore  the  set  (37)  has  zero  measure 

and  by  letting  q=p+1,p€N  we  obtain  the  proposition. 

It  is  interesting  to  compare  proposition  5  with  [8,  Lemma 

2.1]  where  examples  of  numbers  which  are  not  2-well  approximated 

by  En(n=2)  are  discussed.  It  might  be  that  these  include  all 

rationals  ^  1  or  all  quadratic  irrationals  1,  but  no  such  theorem 

seems  to  be  known.  Also,  it  is  easy  to  see  that,  for  N  <  1, 

every  positive  number  is  N-well  approximated  by  ER. 

Proposition  6.  -  Let  r.,,  r2  be  two  positive  numbers  such 

that  r  2/r  •]  is  not  N-well  approximated  by  En .  Denote  by 

X  the  positive  zeros  of  y  .  There  is  a  positive  constant 
k  r  1 

C  such  that 


(  38) 


Proof .  Let  us  denote 
that 


°k , n/2  * 


Recall  that  X. 

k 


aj£/r1  and 


const. 


J 

n/2  2 


(r  t) 
2 


n/2 
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From  the  asymptotic  development  (13)  we  have 

a.  .-a.  *=  x  +  0  ( 1  /k  ) 

k  +  1  k 


Hence,  if  k  is  fixed  and  jk  is  chosen  such  that 
is  minimal  we  have 


r  X  -a  . 
2  k  j 


(39)  ek  =  jr^-a.  I  <  */2  +  0(1/k). 

Let  us  distinguish  two  cases:  <  */4  or  not.  In  the  second 

case  we  have 

|cos(r2lk— <n^— **)  1  «  fcos  (±Ek+(  2jk  +  l  )  j  +  o(£))| 

-  |sin(Ek  +  o(  1  /k) }  |  >  ^  +  o(^)  >  cQ  >  0 


for  large  k.  In  this  case  the  asymptotic  development  (12)  gives 
the  estimate 


n  + 1 

|Sr2(Vi  *  Ci  k  2 


for  some  C1  >  0  and  all  large  k. 

By  hypothesis  we  have  that  for  all  j,k 


*  7  <c’  > 


0) 
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Therefore  >  C3/kN-1  (C3  >  0) . 
the  mean-value  theorem  there  is  a 
such  that 


Suppose  also  e <  */4.  By 

i  between  a.  and  r„X. 

3k  J  » 


(r2V 


n/2 


-  r , 


:  .n/2 


(r  X  -a  . 
2  k  j . 


(Recall  J  (a.  )  =  0.)  Note  that  6  =  |£-a.  |  <  c  <  x/4. 

n/t  3k  ^  J  ^  ^ 

Again  by  (12)  we  have  to  estimate 


cos ( £  - 


ir 

7 


(nt2)w 

4 


cos  (±V(2V 


lZ.  i n t 1  ll 

)  2+  4 


(n+2 )ff 

4 


+  o(  i/k) ) 


=  cos  (  ±  +  j  +0  { 1 /k  )  ) 


*  i  cos  6^  +  0(1 /k ) . 


The  n 

C 

1 1  (Xw)|  >  - ‘ - 7  1C,  >  0). 

r2  "  M-1.2UI  4 

k  2 

Since  N  >  1  the  estimate  (38)  holds  in  both  cases.  □ 

Proposition  7  let  f  be  a  function  in 
L*  (B(0,R)),  g  €  I*]  (B(0,R)),  suppg  c  B(0,r),  g  radial. 

For  | x 0 1  <  R  -  r  and  p  <  R  -  r  -  |xQ|  we  have 

(40)  Xp(f*g,xQ)  »  ( X | # j (f , xQ )*g ( • ) ) (y )  (|y|-p) 
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(The  notation  indicates  that  we  are  convolving  in  the  variable 
denoted  by  a  dot). 

Proof .  Recall  that  the  average  can  a*8°  be  computed 

by 

X  (f,xQ)  =  ^  f ( x  Q+Ay ) d  A 

0  ( n  ) 

where  y  is  any  point  with  |y|  =  p,  0(n)  is  the  orthogonal 

group  and  dA  is  the  normalized  Haar  measure.  Let 

♦  (y )  =  (X | %  j (f ,xQ) *g( •  )) (y ) 

we  have 

♦ ( y )  =  /  (/  f(x  +A( y-x ) dA) g (x )  dx 

R°  0 { n  ) 

=  /  (/  f(x  +A(y-x))g(x)dx)dA 

0  (  n  )  R 

Set  u  «  Ax  then  g(x)  *  g(u)  and  dx  =  du  Hence 

♦(y)  *  /  (/  f (x  +Ay-u )g(u )du)dA 

0  (  n  )  R 

*  V£*9,X0)*  □ 

Corollary  8.  Let  g  be  radial  integrable  function  of 
compact  support  and  a  a  positive  number.  Then 
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-  2 ' tt  ^  ^ 


(41  ) 


g  ( x  )* 


n-2 

2 


(y)  «  g(o ) 


2 _ 

n- 2 


( a | x |  )  '  (o|y|) 

Proof.  Let  £  €  Rn  be  any  vector  with  |£|  *  a,  then 


(42) 


(go 


... 


)  (  y )  *  g  (  i  )  e 


x  (  C  |  y ) 


=  g (a  )  e 


iU|y). 


On  the  other  hand 


X  (ei(5l* },0}  =o  (a) 

p'-  '  p 


n-2 

=  -2  2  T ( n/ 2 ) 


J  .(op) 
n-  2 

2 

n-2  * 
(op )  2 


Applying  now  to  (42)  proposition  7  we  obtain  the  desired 

formula(41).  D 

Proposition  9.  Let  r.,,  r2  be  two  positive  numbers  such 

that  r  2 /r  j  is  not  N-well  approximated  by  En>  Let  R  be  any 

number,  ma x ( r ^ , r 2 ) < R<r 1 +r 2 •  Then  there  is  a  non  zero  radial 

function  f  €  c" ( B ( 0 , R ) )  which  is  mean  periodic  with  respect  to 

U  and  u 
rl  r  2 

Proof .  Let  ♦  €  P()0,r  (),♦  $  0  such  that  supp$  £  [R-r^rjf 
It  follows  from  [16,  theorem  2.1  page  247)  that  +  admits  a 
series  development  of  the  form 
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(43) 


♦  (t) 


I 


k>  1 


n  -  2 

<Xkt)  2 


where  A^  *  n/2^ri*  Tf*^s  *s  Sturm-Liouville  expansion  for 
a  boundary  value  problem  singular  at  t  *  0  and  derivative  equal 
to  zero  at  t  *  r^  It  can  be  seen  by  successive  integrations  of 
parts  that 


(44) 


a^  |  =  0(k  ^ )  for  every  p  >  0. 


Since  ^/r^  is  not  N-well  approximated  by  ER  we  see  that 


bk  -  W/V 
2 


satisfies  the  same  estimates  as  a^  (Proposition  6)  Hence  the 
function 


Jn-2(Vx|) 


(45) 


f  ( x)  =  l  b 


k>  1 


n  -  2 


(Xfc  I *1  ) 


«■>  n  . 

is  a  C  radial  function  in  R  ,  f  f  0,  .  And,  from  corollary 

8,  it  follows  that  f  is  y  me  a n-pe r i od i c .  Furthermore 

r  1 


(y  * f ) ( x )  -  I  b  U  (A  ) 

2  k>  1  2 


♦  (  |  x  |  ) 


Jn12(>fc|xll 

2 _ 

n-2 
2 


<Ajx!> 
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which  is  tcro  in  B(0,R-r2)  and  therefore  the  function  f 

restricted  to  B(0,R)  is  y  and  y  itan-per iodic .  □ 

r  1  r  2 

The  above  propositions  can  be  summarized  by  the  following: 

Theorem  10.  Let  r1  >  0  and  r2  >  0  be  such  that 

r  /r  €  E  .  The  necessary  and  sufficient  condition  on  a  open 
2  i  n 

set  Q  of  Rn  so  that  the  only  distribution  T  €  P ' ( ft )  which 

can  be  mean-periodic  with  respect  to  both  y  and  y  is 

ri  r  2 

T  =  0,  is  that  ft  is  the  reunion  of  balls  of  radii  strictly 
larger  than  r1  +  r2* 

An  amusing  corollary  of  theorem  10  is  the  following: 

Corollary  11  If  z2^r\  *  E2 '  r1  +  r2  <  R  and  f  ^  C(B(0,R)) 

then  the  conditions 

/  f(C)dc  =  0  for  every  z,  |z|  <  R  -  r.(j=1,2), 

3B(z,r  .  )  -1 

imply  that  f  is  holomorphic  in  the  disk  B(0,R). 

5 .  Generalizations 

After  the  paper  was  written,  we  became  aware  of  the  work  of 
J.D.  Smith  [  17]  ,  in  which  local  versions  of  certain  two-circle 
theorems  are  also  proved.  Smith's  results,  which  require 
R  >  2r^  +  r2  are  less  sharp  than  Proposition  4;  nor  does  the 
method  cT  proof  seem  to  generalize  to  the  other  problems  discussed 
in  [20]  ,  e.g.  the  converse  of  the  mean  value  property  for  harmonic 
functions.  The  aim  of  this  section  is  to  show  that  the  methods 
used  above  do  generalize. 
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Definition.  We  say  that  a  radial  distribution  y  of 


compact  support  is  hyperbolic  if: 

(i)  y  is  invertible,  and 

(ii)  there  is  a  constant  C  such  that  every  zero 
X  of  y  satisfies 


| ImX  j  <  C  log ( 2+ | X  |  )  . 

Theorem  1  2  Let  li  1  ,  y  .  .  .  be  a  (possibly  infinite)  family  of 

radial  distributions  of  compact  support,  cv{suppw J  =  B ( 0 , r  ) . 

Suppose  { z£Cn  :y_.  (z  )=0  Vj  }  =  ft,  y  ^  is  hyperbolic,  and 

R  -  r1  >  sup  r^ .  Then  {  f €  V  ( B ( 0 , R ) :  y  ^  *f  =  0  Vj }  =  {  o} . 

Proof :  Due  to  the  condition  on  R  we  can  assume 

f  C  c"°  (B(0,R)  )  as  done  before.  The  proof  that  leads  to  (32) 

can  be  repeated  almost  verbatim  just  using  for  each  Xfc,  zero  of 

y  ,  a  convenient  y.(j>2)  with  y.(X,  )  f  0.  We  obtain 
1  3  3  K 

(47)  Tu  *  f ( x )  -  0  for  |x|  R  -  r  -  sup  r  , 

x ,  s  1  3  3 

where  T.  (t)  =  y  ,  ( t )  (t2-X2)  ,  1  <  s  <  m  ,  nw  = 

K  ,  S  1  X  x 

A* 

multiplicity  of  X^  as  a  root  of  y^ 

On  the  other  hand 

(48)  (-1  ) 8 ( A  +  X 2  )  8  (T^  s*f)  -  y^f  -  0  in  B(0,R-r1), 

therefore,  T.  *f  is  real  analytic,  and  hence 
k,s 
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(49) 


0  in  B( 0 , R-r  ) 


T.  *  f  * 
k  ,s 


as  before.  It  is  at  this  point  we  have  to  be  more  careful  to 
prove  the  correct  version  of  Proposition  2.  It  will  be  replaced 


by  the  following 


Lemma  13  Let  A  =  {A  }  =  set  of  distinct  zeros  of  y  , 

-  V  '  i 


then 


A  =  U  A.,  where  the  A.  are  finite  and  mutually  disjoint 
j=0  3  3 

sets.  There  is  also  a  positive  integer  q  such  that  for  any 
p,  0  <  p  <  •  we  can  write 


(50) 


o  =  v  +  y  *  S  , 
P  pip 


where  v  ,  s  are  radial  distributions  satisfying 
P  P 


(51  ) 


supp  vp  £  B(0,r 1  ) 


(52) 


supp  S  c  B ( 0 ,  Max ( r 1  ,  p )-ri  ) 


Furthermore  , 


(53) 


I  A* 


a  convergent  series  in  E  ( R° ) ,  each  v  a  finite  linear 

o  Pi  j 

combination  of  the  distributions  T  ,  A  €  A  ,  1  <  s  <  m  (if 

JC  ,  S  K  j  X 

n  >  1  then  one  denotes  by  Aqv  not  only  a  finite  linear 

O  PiO 

combination  of  A**  T.  ,  A  i  A  but  also  of 


T  ,  AT 
o .  ■  o  ,  ■ 

o  o 


k ,  s  k  o 

,Aq-T  .  .) 
o  •  ■ 

o 
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Once  this  lemma  has  been  proved,  the  proof  of  Theorem  12  is 
achieved  the  same  way  as  it  was  done  in  Proposition  4  and  we  note 
that  the  hypotheses  imply  2r1  <  R. 

Proof  of  Lemma  1 3  The  proof  of  this  lemma  proceeds  as  in 
Proposition  2  by  interpolating  the  values  of  op  on  the  variety 
of  zeros  of  y  (counted  with  multiplicities).  We  have  to 
repeat  with  due  care  the  procedure  used  in  (13),  (15),  (18)  since 

we  need  the  precise  statement  (51),  (52),  and  (53). 

First  we  note  that  as  in  (5,  lemma  4)  (cf.  also  (11,  p. 

50)),  the  condition  of  hyperbolici ty  and  the  minimum  modulus 
theorem  allow  us  to  construct  a  family  of  a  Jordan  quadrilaterals 
r ,  k  €  Z  symmetric  with  respect  to  the  real  axis  and  enjoying 
the  following  properties: 

(54)  for  some  d  >  0  the  horizontal  sides  lie  on  the  curves 

Im  z  =  ±  log (d+ f Rez  j  )  , 

and  the  vertical  sides  are  arcs  of  circles. 

(55)  0  €  int  which  is  symmetric  with  respect  to  the  origin 

(i.e.  if  z  €  T  then  -z  €  T  also). 

o  o 

(56)  for  k  j  0,  r  is  the  symmetric  of  with  respect  to 

the  origin. 
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(57)  for  j  f  K  >  int  r  fl  int  I*  ^  furtheriore ,  for  some 

positive  number  a  we  have  that  if  z  f  r., 

dist(z,r  )  >  (a+|z|)  *  for  any  fc  yf  j . 

(58)  for  some  positive  constant  b  we  have: 

diam  f_.  <  b(1  +  |z|)b 


and 


length  <  b  (1+|z])^, 


for  any  z£  int  T,  ,  any 


j  • 


(59)  there  is  a  constant  c  >  0  such  that  for  any  j,  and  any 

z  €  T  .  we  h  a ve 
3 

| W 1  ( z )  |  >  (c+  |  z  |  ) “C  , 


and  this  inequality  is  valid  even  for  those  z  such  that 
dist(z,Tj)  <  1/2  (a+|z|)  3  (the  same  a  as  in  (57)) 

m 

(60)  A  c  u  (int  r  )  . 

j—  j 

(61)  for  some  d  >  0:  if  j  >  1,  z  €  ,  then  |z|  >  j/d. 

(62)  A  -  A  int  r,  A.  -  A  (int  I*  n  int  r  . ) .  j  >  1. 

o  o  3  j  -J 
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For  the  sake  of  definiteness  we  will  index  the  points  in  A  so 

that  A  *  0,  and,  for  k  >  1,  either  Re  A  >  0  or  ReA.  ■  0 
o  k  k 

and  Im  A^  >  0,  and,  finally,  A  =  -A^. 

Now  consider  the  even  entire  function 


(63) 


2  q  ** 

f (t)  =  t  4  w1 (t) , 


for 


a  positive  integer  to  be  chosen  conveniently  later  on. 


We  note  that  if  t  t  int  T.  U  int  I*,  then 

T  3  3 


(64) 


♦  j  (t) 


2*i  [ 

3 


/  -E— 
Jr  f(s) 


(s)  ^ 
_  ds 


— - —  +  -  f  — E - 

s- 1  2*1  j,  f  (  s  ) 


(s)  , 

) _  ds 


s- 1 


(where  we  disregard  the  second  term  if  j  »  0)  is  an  even 
function  which  is  a  linear  combination  of  terms  of  the  form 


2  2  -  s 

(t  -A  ) 
k 

1  <  s  < 


for  A,  €  A.  and  1  <  s  <  m,  if  k  >  1  , 

k  j  k 

>  +  2q  if  .  k  =  0.  Hence  6.  can  be  defined  as  a 

o  3 


rational  function  throughout  t  and  the  function  f(t)  $..(t)  is 
an  even  entire  function.  We  want  to  show  now  that  q  can  be 
chosen  so  that 


(65) 


g ( t )  -  l  f(t)  ♦  (t) 

3-0  J 


is  in  E  ( R  )  and  the  series  converges  in  the  topology  of 
o 

*w  *  n 

E  ( R  ) . 

o 

In  fact,  we  have  that  for  |lm  t|  <  log (d+ | Re t | )  there  is 
some  N  >  0  such  that 
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(66/ 


|op(t)  |  <  C ( p )  ( 1+  |t|  )N 


and  also 

(67)  | y 1 (t) |  <  Co(1+|t|)N. 

Therefore,  for  some  N1  >  0  sufficiently  large,  if 
dist(t,  int  r  U  int  T  )  >  1  we  have  by  (66),  (59)  and  (58), 
that  with  respect  to  an  arbitrary  point  z  €  int  r  ,  which  we 
can  take  it  to  be  the  point  in  the  positive  real  axis  closest  to 
the  origin. 


|  4»  j  ( t )  I  <  (Nj  +  |z|)  ^zl  2q  <  const,  j 


by  (61)  (just  take  2q  >  N1  +  2).  Therefore,  under  the  same 
condition  on  t  we  have 


(68) 


|f  (  t)  ♦  (t)  | 


<  C. 


.  -2 
3 


d+|t|)M 


r 1 I Int | 
e 


Using  the  condition  (58)  on  the  diameter  of  r  and  (67),  this 

estimate  remains  valid  throughout  £,  after  possibly  increasing 

C, ,  M.  This  shows  that  the  v  €  E  ( Rn )  defined  by 
i  p  o 


V  (t)  -  g ( t ) 

P 

satisfies  (51).  It  is  also  clear  that  the  distributions  v 

P.  j 

such  that  t2<*  v  «  f(t)  ♦.(t)  have  the  properties  required 

P  #  j  3 
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To  end  the  proof 


by  (S3)  (with  special  care  taken  if  m  j*  0). 


of  the  lemma  we  only  have  to  show  that  g(t)  -  o  (t)  is 

P 

divisibly  by  w^ft),  the  rest  is  the  same  as  in  Proposition  2. 

Note  that  if  t  €  int  T.  U  int  T  then  we  have 

3  -3 


♦  j  ( t ) 


1  r  a ( s )  ds 

2*i  p  +r  f(s)  s-t 

j  -  j 


is  holomorphic  and,  we  pick  a  new  residue  at  s  =  t  yielding 


f (t)  (t)  =  f (t)  $ . (t)  +  o(t) . 

3  3 

This  concludes  the  proof  of  Lemma  13  and  Theorem  12.  □ 

We  give  here  the  local  version  of  Delsarte's  theorem  for 
harmonic  functions. 

Corollary  14  Let  H_  =  {  C/n  :  £ /n  €(  0  ,•  )  ,  o.(£)  =  o.(n)=l|.  If 
-  n  1  1 

R  >  r j  +  r2»  ri/r2  i  Hn>  an<^  u  is  a  continuous  function  in 
B ( 0 , R )  satisfying 

u(x)  =  X  (u,x),  for  all  x,  |x|  <  R  -  r  , 

r  1  1 

and 

u(x)  -  X  (u,x),  for  all  x,  |x|  <  R  -  r  , 

T  2  2 

then  u  is  harmonic  in  B(0,R). 

Proof  t  Prom  the  asymptotic  development  of  the  Bessel  functions 


end  the  formula  (10),  it  follows  that  the  radial  distributions  y 


defined  by 


y  .  (t) 
J 


j 


<or  ( t ) - 1 ) /t 2 
3 


are  hyperbolic.  The  hypothesis  on  r,,/r2  guarantees  these  two 
entire  functions  have  no  common  zeros.  Theorem  13  shows  now  that 
the  distribution  Au  is  zero  in  B(0,R). 

Remark  As  mentioned  in  [20],  Delsarte  proved  this  theorem  in  Rn . 
He  also  showed  that  Hr  is  finite  and  =  {l}.  Hence,  at  least 

for  dimension  3,  any  pair  of  distinct  positive  value  r^,  r2 
would  work  in  the  above  corollary. 

The  several  other  results  in  [20]  can  now  be  carried  over  to 
the  local  case  without  difficulty.  It  remains  as  an  open 
question  for  the  moment  the  elimination  of  the  in ve r tibi li ty 
condition  on  y1  ,  which  could  probably  be  done  following  the 
Euclidean  summation  method  of  [6].  More  interes ting, in  our  view, 
is  to  try  to  extend  this  theorem  to  non-compact  symmetric  spaces 
of  rank  1  or  even  to  the  Euclidean  group  thus  obtaining  a  local 
version  of  the  Pompeiu  problem  considered  in  [9]. 

As  an  example  of  this  let  us  mention  the  following 
corollary  of  Theorem  13. 

Corollary  15  Let  R  >  /n  a,  if  f  €  l|  (B(0,R))  has  zero 
integral  over  any  n-cube  of  side  a  contained  in  B(0,R),  then 

f  ■  0  a  .  e . 

Proof  Following  the  ideas  from  [9]  we  see  we  can  consider  all 
radial  distribution*  u  whose  Fourier  tranforms  are  of  the  form 
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(69)  /  XqUOTUO  dk  -  1(0 

O  (  n  )  W 

where  Q  is  the  cube  (~a/2,a/2]n  and  T  is  a  distribution  of 
compact  support  in  the  ball  B(0,e),  t  +  /n  a  <  R.  Then,  for  any 
such  \t  ,  cv(Suppu)  B(0,r),  and  f  will  satisfy  the 
equations : 

W  *  f  =  0  in  B(0,R-r). 


Since  this  set  of  distributions  generates  the  same  closed  ideal 
in  E  ( R° )  as  those  are  considered  in  [9,  p.  602],  then  their 
Fourier  tranforms  have  no  common  zeros  [9,  section  9].  It  only 
remains  to  find  a  distribution  that  plays  the  role  of  y  1  in 
Theorem  13.  The  easiest  one  is  obtained  when 


An  easy  computation  shows  that  in  this  case,  for 


IJ  =  average  over  0(y)  of  - - - —  Xq»  we  have 

3  X1  ' *  *  3  Xn 

(70)  U^tt)  =  const.  t^n^2*  +  1j3n  2^n  a*/2) 

2 

which  is  clearly  hyperbolic.  (For  n  ■  2,  this  can  be  obtained 
from  Soni  le  second  finite  integral  (19,  p.  376].) 
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CHAPTER  4 


INVERSION  OF  THE  LOCAL  POMPEIU  TRANSFORM 


by 


*  **  ** 
C.A.  Berenstein  ,  R.  Gay  and  A.  Yger 


To  appear  in  J.  Analyse  Mathematics 


4.0 


I 


1 .  Introduction. 

The  prototype  of  the  Pompeiu  transform  and  of  the  problems  we  will 
consider  here  Is  the  Two-Disks  problem  we  discuss  presently.  Let  E  be  the 
collection  of  positive  quotients  of  zeros  of  the  Bessel  function  J^.  This 
countable  set  is  precisely  the  exceptional  set  for  the  following  theorem 
(proved  independently  in  (12]  and  [21]). 

Theorem.  Let  r  r  >  0,  r  /r  *  E.  Then  a  continuous  function  In  the  plane 
is  identically  zero  if  and  only  If  all  Its  averages  over  any  disk  of  radius 
either  r  ^  or  r vanish. 

If  we  let  * , *  be  the  characterist 1c  functions  of  the  disks  B ( 0 . r ^ ) , 
BfO.r^)  we  are  saying  that  the  map 

(1.1)  P  :  f  i — > (x1*f.X2*f) 

is  injective.  In  (7],  (8J  it  was  shown  that  a  stronger  arithmetical  condition 
was  necessary  and  sufficient  for  the  existence  (and  the  explicit  construction) 
of  two  distributions  vyv2  comPact  support  such  that 

(1.2)  ul*(*l*f  *  +  P2*(*2*f  *  =  f 

Because  be  distributions  (the  map  P  Is  smoothing),  it  is  more 

/  2  00 

natural  to  consider  P  in  the  Frechet  space  6  =  6(R  )  of  C  functions. 

Then  we  have  that  r\/'r2  *  ^  and  only  if 

P  : 

is  injective.  The  existence  of  vi>v2  solving  (1.2)  is  equivalent  to  Im  P 

being  closed,  hence  P  is  right  invertible.  One  can  in  this  case  also  show 

that 

(1.3)  Im  P  *  { (g, h)  €  (£(R2))2  :  p2«g  =  |ij*h  > 
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» 


and 

P_1  :  In  Pt—>S 
(g.h)  i — *u1*g  +  i>2*h. 

The  exact  arithmetical  condition  is  Just  the  existence  of  a  constant 
A  >  0  such  that 

(14)  *rl/r2  ”  *  X71  A* 

for  any  pair  £,v  >  0,  Jj(£)  =  (i>)  =  0. 

Since  (1.2)  is  an  explicit  deconvolution  formula,  L.  Zalcman  asked  in 
[20]  whether  some  sort  of  explicit  reconstruction  formula  would  also  exist 
under  the  sole  condition  r^/r^  *  E.  The  theorem  below  answers  this  question 
in  the  affirmative.  Naturally  coupled  to  this  question  is  the  following  local 
problem.  Suppose  we  know  the  averages  of  f  on  disks  of  radii  r ^  and  r^. 
only  when  those  disks  lie  in  some  fixed  disk  B(0,R),  could  these  data  deter¬ 
mine  f  in  B(0,R)?  A  priori,  formula  (1.2),  if  valid,  would  only  determine 
f  in  B(0,R')  for  some  convenient  R'  <  R.  This  is  just  a  consequence  of 
the  fact  that  (1.2)  is  given  by  convolutions,  and  do  not  have 

support  at  the  origin.  This  is  the  nature  of  what  was  called  elsewhere  the 
local  Pompeiu  problem.  It  is  shown  in  [21  that  if  r^  +  <  R,  rj/r2  * 

then  f  is  determined  in  the  whole  disk  B(0,R)  by  those  local  averages.  We 
will  show  here  that  one  has  even  a  local  reconstruction  formula. 

These  local  theorems  are  also  valid  for  other  sets  instead  of  disks,  when 
one  takes  all  rigid  motions.  The  injectivity  was  shown  in  [3J.  We  give  here 
a  method  to  construct  a  local  inverse  in  this  more  general  situation. 

Finally,  we  would  like  to  mention  that  even  for  the  injectivity  question 
2 

in  F  ,  it  is  very  hard  to  see  when  one  can  work  just  with  translations  of  a 
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finite  family  of  sets.  The  antecedent  of  several  such  results  Is  the  Three- 
Squares  Theorem  (SI  where  It  Is  shown  that  If  the  averages  of  a  continuous 
function  f  In  the  plane  vanish  on  every  square  of  sides  either  r^r^  or 
r^  (with  sides  parallel  to  the  axes),  then  f  is  Identically  zero  if  and 
only  if  ri>r2,r3  are  O-linearly  Independent.  In  [8],  [9],  explicit  decon¬ 
volution  formulas  were  found  for  this  case  under  extra  arithmetical  assump¬ 
tions  on  the  triple  ri,r2,r3'  ^he  *oca*  version  of  the  Three-Squares  Theorem 
requires  that  we  place  ourselves  In  a  square  of  side  R.  In  14]  it  is  shown 
that  if  R  >  rj + r2 + r3’  then  the  uniqueness  still  holds.  Its  proof  is 
akin  to  the  methods  used  In  this  paper. 

The  motivation  for  this  work  lies  in  trying  to  find  an  algorithmic  decon¬ 
volution  approach  (with  due  care  for  error  bounds  and  noise  behavior)  in  a 
situation  where  part  of  the  scene  is  obstructed  from  our  view  but  the  object 
we  are  looking  for  lies  very  close  to  this  obstruction.  In  the  case  of  the 
Radon  transform  this  is  sometimes  called  the  Hole  problem  or  the  Bagel  problem 
[15],  It  is  interesting  to  note  that  these  two  problems  are  related,  since 
one  way  to  find  an  inversion  formula  is  to  use  Cormack-Quinto’ s  Spherical 
Radon  Transform  and  its  explicit  inverse. 


2.  Preliminaries. 

We  recall  some  notation  and  basic  properties  of  the  Pompeiu  transform 
Let  be  a  collection  of  compact  sets  in  Fn  of  positive  measure 

and  let  M(n)  be  the  group  of  Euclidean  motions  in  Rn,  then  the  (global) 
Pompeiu  transform  (associated  to  Ej . 

(2.  1)  P  :  C(Fn)-»C(M(n))N 


E^)  is  the  map  (10) 
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given  by 


(2.2)  (Pf)(g)  - 

Civen  am  open  set  U  £  F°  we  can  define  open  sets  £  M(n)  by 

(2.3)  Cj  =  (g  €  M(n)  :  gEj  £  U) 

(these  sets  could  be  empty).  We  define  the  local  Pompeiu  transform 


* 

f 

fdx . 

]  fdx 

L  gEl 

*en 

N 

(2.4)  P  :  C(U) — »  ©  C(G  ) 

J=1  J 


defined  exactly  by  the  same  formula  (2.2).  The  family  . E^  is  said  to 

have  the  Pompeiu  property  (respectively,  the  local  Pompeiu  property  with 
respect  to  U)  whenever  the  map  (2.1)  (respectively  (2.4))  is  injective.  The 
most  interesting  case  is  the  case  of  a  single  set  with  relatively  nice 
boundary.  One  cam  prove  the  following  theorem. 

n  —  c 

Theorem  1  [1],  [6].  Let  Q  be  an  open  bounded  set  in  R  ,  E  =  n,  E 
connected  and  dE  Lipschitz.  Then  E  (or  Q)  has  the  Pompeiu  property  if 
and  only  if  there  is  no  positive  eigenvalue  for  the  overdetermined  Neumann 
problem 


(N) 


fin  ♦  ou  =  0 


in  Q 
in  an. 


As  a  corollary,  one  obtains  two  classes  of  examples: 

(1)  A  single  ball  never  has  the  Pompeiu  property; 

(2)  A  set  n  satisfying  the  hypotheses  of  Theorem  1  whose  boundary  is 
not  a  real  analytic  hypersurface  has  the  Pompeiu  property. 

Examples  of  sets  having  the  Pompeiu  property  and  also  having  real  analy¬ 
tic  boundary  exist,  for  instance,  any  ellipsoid  which  is  not  a  ball.  Since  a 
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single  ball  does  not  have  the  Pompeiu  property,  the  following  theorem  gives 
the  best  possible  answer. 

Theorem  2  (12],  (211.  Let  Z  =  {£/tj  :  £  >  0,  tj  >  0,  J  ._(£)  =  J  frj)  =  0). 

n  n/d  n/2 

A  pair  of  balls  B(0,r^),  BfO.r^)  has  the  Pompeiu  property  if  and  only  if 
rl/r2  *  Z„- 

To  consider  the  local  Pompeiu  property  it  is  clear  we  have  to  start  with 
a  family  that  satisfies  the  global  Pompeiu  property.  We  restrict  ourselves  to 
open  sets  U  that  can  be  covered  by  balls  of  a  fixed  radius  R.  In  this  case 
we  might  as  well  suppose  that  U  is  the  ball  B(0,R).  In  this  setting  it  Is 
easy  to  state  the  local  version  of  Theorem  2. 

Theorem  3  (2].  Let  r  ,r  >  0,  r,/r„  «  Z  and  R  >  r,  ♦  r_.  Then  the  pair 

id  id  n  12 

B(0,r1),B(0,r2)  has  the  local  Pompeiu  property  with  respect  to  the  ball 
B(0, R) .  Under  the  additional  restriction  that  r^/r^  is  badly  approximated  by 
elements  of  Z^,  the  condition  R  >  is  also  necessary. 

It  is  quite  possible  that  R  >  r^  + r^  is  always  a  necessary  condition 
in  Theorem  3.  When  considering  a  single  set  E  having  the  global  Pompeiu 
property,  in  order  to  decide  whether  E  has  the  local  Pompeiu  property  with 
respect  to  B(0,R),  it  is  natural  to  try  to  measure  R  against  a  value  r 
such  that  E  S  B(x0>r)  for  some  xQ.  The  condition  that  replaces  rj  + r2  < 

R  is  2r  <  R.  We  are  compelled  to  introduce  an  extra  technical  condition  due 
to  the  fact  that  we  know  nothing  about  3E.  This  condition,  which  we  will 
explain  below,  is  called  hyperbol ici ty.  There  a re  a  number  of  simple  condi¬ 
tions  that  imply  it  (see  (3]),  for  instance, 

(1)  Near  an  extreme  point  xJ(  x}  c  E  n  3B(xQ,r),  E  coincides  with  a 
polyhedral  angle  with  vertex  Xj.  The  fact  that  the  walls  or  edges  are 
straight  plays  no  role.  For  instance,  E  is  a  cube  in  Rn. 

(2)  There  is  an  extreme  point  xJt  Xj  €  E  n  3B(x0,r),  such  that  Xj  is 
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a  point  of  strict  convexity  for  flE  and  OE  Is  of  class  C  near  Xj. 

2 

(3)  In  R  ,  near  a  point  xJt  Xj  *  E  n  3B(x0>r),  8E  Is  a  Jordan 
curve,  sufficiently  smooth. 

Theorem  4  [ 3 ] .  Assume  E  £  Bix^.r)  has  a  hyperbolic  point 

Xj  €  E  r\  3B(x0>r).  Assume  further  that  E  has  the  global  Pompeiu  property 

and  2r  <  R.  Then  E  has  the  local  Pompeiu  property  with  respect  to  B(O.R). 

We  remark  that  the  conditions  rj  * r2  <  ^  in  theorem  3  and  the  corres¬ 
ponding  2r  <  R  in  Theorem  4  mean  that  the  open  sets  C.  (respectively  g) 

J 

of  M(n)  are  not  to  small.  On  the  other  hand,  they  impose  no  restriction, 
when  U  is  the  complement  of  a  closed  convex  set  K.  This  is  akin  to  say 

that  the  Pompeiu  transform  possesses  the  Hole  property,  as  is  the  case  for  the 

Radon  transform  (in  that  case  K  must  be  compact). 

The  proofs  of  Theorems  1  and  4  depend  on  the  introduction  of  a  countable 
family  of  radial  distributions  associated  to  the  set  E.  For  that  purpose  we 
need  to  introduce  the  concept  of  the  radialized  version  of  a  distribution  (or 
circular  symmetrizat ion)  with  respect  to  a  point  xQ.  Namely,  let  T  e  D'(Rn) 
then  R  T  =  XT  is  defined  by 

(2.5)  <R  T,^?>  =  <T(x),  0(k(x-x_  )+X-)dk>, 

Xq  u  u 

JS0(n) 

where  0  e  7)(IRn),  dk  is  the  Haar  measure  on  S0(n).  In  case  T  =  x £.  E  a 
compact  set  in  IRn,  then  R  T  is  a  function  of  compact  support  with  values 

Xq 

in  [0,1]  which  is  a  radial  function  with  respect  to  x  ,  R  T(x)  = 

0  Xo 

tpi  lx-x0l ), 

(2.6)  f>( t )  =  <r(E  n  8B(x  ,t)), 

where  o -  Is  the  normalized  surface  area.  We  consider  tp  as  a  radial  Borel 
function  in  Rn,  of  compact  support,  and  therefore  its  Fourier  transform 
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Is  an  entire  function  In  Cn  in  the  Paley-Wlener  class,  which  is  radial, 
i  .  e. ,  if  we  let 

(2.7)  ♦(£)  =  ?<P(C.0 . 0) 

then  *  is  an  entire  function  in  C,  even,  satisfying  the  inequality 


(2.8)  l*(<)l  £  A  exp(r|Im  <|). 

whenever  E  £  B(xQ,r).  In  fact,  #  is  the  Fourier-Besse  1  transform  of  the 
radial  function  <f> 

00 

*(C)  =  (2n)n/2  <pi  t)  J(n_2)/2{<t)tnldt’ 

•*0 

where  j.(z)  =  J^(z)/z^,  being  the  Bessel  function  of  order  l.  We  also 
denote  ♦  =  <p.  Conversely,  ♦  determines  by 

(2.10)  MCj . = 

The  smallest  r  which  appears  in  (2.8)  determines  the  smallest  disk  of  center 

_  mm 

xQ  such  that  E  £  B(xQ,r).  We  call  exterior  radius  of  E,  r  =  r  (E),  the 
minimum  value  of  r,  with  respect  to  all  xQ,  such  that  this  inclusion 
holds. 

For  fixed  x^  and  any  a  €  INn,  we  can  construct  a  corresponding  entire 

function  of  one  complex  variable  9^  by 


(2. 11 ) 


oc  d 
where  D  =  - 


l«| 


3x 


•  CO  ==  (*  (D %))  (O. 

a  x0  E 


It  is  convenient  to  introduce  the  distributions  p  = 

a 


!R  D V,  which  are  radial  with  respect  to  xn.  These  definitions  are 

Xo  E  u 

justified  by  the  following. 


Lemma  4  1121.  E  has  the  global  Pompeiu  property  If  and  only  If  the  functions 
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•  ,  a  (  l/1,  have  no  common  zeros. 

a 

An  even  entire  function  I  of  a  complex  variable  is  said  to  be  hyper¬ 
bolic  if  the  following  three  conditions  are  satisfied. 

(i)  For  some  constant  A  >  0.  it  satisfies  the  estimate 

!*(<)!  <  A( 1*ICI )Aexp(A|Im  <|). 

(ii)  For  some  constant  a  >  0  and  any  real  value  £ 

max{  |#(£+tj)  |  :  1 7j  |  <  a  log(2-*j£l)>  £  (a+|£|)  a. 

(iii)  For  some  constant  c  >  0  we  have  that 

0(0  =  0  =>  |  1m  Cl  *  c  log(2+|<|  ). 

The  first  condition  is  simply  that  *  €  ?(£'(R)).  The  second  is  usually 
called  invertibility.  it  means  that  the  principal  ideal  *?(£'(R))  Is  closed. 
It  coincides  with  the  property  that  4>H(C)  n  ?(g'(R))  =  4>?(6'  (R)).  The  third 
one  says  that  all  the  zeros  of  ♦  are  almost  real. 

In  [3],  it  is  shown  that  a  number  of  natural  conditions  on  E  imply  that 
the  radial ized  function  of  x^>  is  hyperbolic. 

The  idea  of  Lemma  4,  and  a  fortiori  Theorems  1  and  4,  is  that  everything 
is  reduced  to  the  study  of  the  radial  case,  as  follows.  Let  us  fix  x^, 
which  we  take  to  be  the  origin  for  simplicity.  Then  the  condition 

fdx  =0  V  g  e  M(n)  (or  V  g  €  G) 

JgE 

can  be  rewritten  as  an  infinite  system  of  convolution  equations 
(2.12)  V  p  €  S0(n) :  x  F*f(x)  =0  V  x  €  Rn  (or  |x|  <  R-r). 

Pt. 

The  proof  of  the  injectivity  of  the  (local)  Pompeiu  transform  is  bused  on 
the  principle  that  if  there  is  a  non-zero  solution  of  (2.12),  we  can  assume  it 


to  be  C*°  and  show  there  Bust  be  a  non-zero  radial  C*  function  g  In 
(respectively  In  B(0,R)1  solution  of  the  denumerable  system  of  convolution 
equations 

(2.13)  V  a  €  Nn  p  *g(x)  =0  V  x  e  Rn  (or  |x|  <  R-r). 

a 


where  the  u  are  the  distributions  defined  In  (2.11).  In  particular,  they 
a 

involve  derivatives  (see  112),  13]  for  the  details). 

In  order  to  find  an  inversion  formula  for  the  local  Pompclu  transform  we 
must  make  this  principle  a  bit  more  precise.  This  depends  on  the  following 
lemma. 


Lemma  5.  Let  E  £  B(0,r)  and  U  =  B(0,R),  R  >  2r.  Then  for  every  x  such 
that  | x I  <  R-r,  a  e  Nn,  and  any  f  €  £(B(0,R))  we  have 


(2. 14) 


p  •f(x)  = 
a 


S0(n) 


<6(a)(y),P(f )  q  Xlky j)>dk 


where  M(n)  is  considered  as  the  group  of  (n+l)x(n+l)  matrices  of  the  form 
*||,  k  e  SO  ( n ) ,  x  e  Rn,  and  Rn  is  identified  to  the  affine  subspace 


(x  =  1}  of 
n+1 


r,n+l 


Proof.  The  identity  (2.14)  is  an  easy  computation  and  the  only  thing  to 

k  ^  x— ky 

q  1  S  is  con“ 
tained  in  B(0,R)  (recall  that  supp  6(°^(y)  =  (0)).  This  is  evident  since 
(-k  J)E  £  B(0,r).  Hence,  translation  by  x  keeps  it  inside  B(0,R). 


Given  a  compact  set  E,  let  E^  =  (x  +  t  :  x  e  E)  the  translate  of  E 
by  x  e  IRn.  Consider  the  functions 


0  (t)  :=  (K  E  )  (t). 
T  xo  T 


It  is  easy  to  see  they  are  real  analytic  with  respect  to  the  parameter  t.  In 
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analogy  with  Lemma  4  we  have 

Lemma  4 * .  E  has  the  global  Pompelu  property  If  and  only  If  there  Is  an  open 
set  J  £  Rn  such  that  the  family  of  functions  has  no  common  zero. 

Another  minor  technical  lemma  will  be  useful  to  us. 

Lemma  6.  Let  P  be  a  homogeneous  harmonic  polynomial  of  degree  k  and 
the  normalized  measure  on  the  sphere  8B(0,p)  (0  <  p),  considered  as  a 

radial  distribution  in  pn.  Then,  with  |x|  =  r,  we  have 


,k  2-n 


(2.15)  P(*)<Tp(r>  -  ^4^77  TT-  r>(fc]np2-lx|2)k-1lB(0  P,(x)l. 


where  fi  =  2nn/^/r(n/2)  is  the  surface  of  the  unit  sphere  In  Rn. 
n 


Proof.  If  we  consider  the  Fourier  transform  of  P(x)cr  (r).  we  have 

P 

*<P<*>yr»«>  -  p[ig^ . 

.k  An-2)/2rln.ti[d  1  .  Pi  "772. 

■  1  2  r<2)P[^JJ<„-2)/2(p'V '•  ■‘V 


2)/2<pvV'- •*<>• 


where  j^(z)  =  Jy(z)z  We  can  now  apply  an  identity  from  [16,  p.  126}  to 
obtain 


,fej  . 

[sgJ(n- 


2i/2(p/<6  =  p(c>[s  In]  ■»<„-; 


(-1|kp2kP(OJ(n/2))k.1(p/^...K7 


[2  2 

On  the  other  hand,  the  function  j  (  n/2 ) +  k- 1  ^ pv^  1 + '  '  ' +l^n  ^ 
transform  of  a  rather  simple  function.  Namely,  if  t  =  |^|, 


is  the  Fourier 


,  2  .  ,2.k-l  ,  .  -ix-C  , 

(P  ~|X|  ]  ZB(0,p)(x)e  dx 
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r,  2  2,k-l  , 

(p  -r  ’  J(n-i 
J  A 


,  .  x  n/2 . 
,(rt)r  dr 


(n/2)*2k-l  f  _  .  , 

(2">  P,Tn^7F-  <>-*  >  J.„-2>/2«>ts)s  ** 

1  Jo 


.n/2„k-l,.  n*2k-2  .  , 

(2n)  2  ( k- 1 ) ! p  J(n/2)+k-l(pt)' 


as  can  be  seen  from  [14,  p.  688,  formula  6.567.1).  Hence  If  we  let 


TOO  =  Tpik(x)  =  <P2-lxl2)k-1lB(0  p)(x). 


we  have 


9(P(x)«r  (r) )  (<)  =  (-1  )kik  - p2~nP(C)T(0 

p  nn/22k(k-D- 

=  (-nk  rfn/2)P -  ,(pf|-]T(r),(a. 

nn/22k(k-l)! 


Therefore 


P(x)tr  (r) 
P 


n  _k 

n  2 


I-”11  pfS-.lT 


which  is  the  formula  (2.15). 


In  what  follows  we  will  need  to  determine  the  radial  distribution  in  R  , 

u  =  p  ,  ,  that  has  as  Fourier-Bessel  transform 

n,  k,  a,  p 


(2. 16) 


>*(l)  =  J(nX2)*k-l(pU/(t2‘“2)’ 


where  a  is  a  zero  of  the  numerator.  It  is  easy  to  find  the  even  function  <f> 
in  IR  whose  Fourier  transform  coincides  with  p.  Therefore,  what  we  need  is 
to  have  a  systematic  way  of  obtaining  p  from  ip.  In  fact,  this  type  of 
relation  underlies  all  the  proofs  about  the  local  properties  of  the  Pompeiu 
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transform. 


We  have  the  following  setup.  Let  £q(F  1  be  the  space  of  radial  dlstrl 
n  rj  * 

butlons  In  F  and  {S' (R) n)  the  space  of  entire  functions  In  the  Paley- 

Wiener  class  which  are  radial  when  restricted  to  Rn.  Then,  if  9  denotes 

n 

the  Fourier  transform  In  Rn,  9£  the  Fourlei — Bessel  transform  in  Pn,  we 

n 


(2. 17) 


S'  (Rn) 
o 


(S' IF  )) 


S'  (IR) 
o 


-i - >  (G'(F)) 

M  ° 


where  SF(t)  =  F(t,0,...,0)  (=  F(<)  for  any  C  e  C  such  that 

< j  .  .  .  ♦  =  t ) ,  S  V )  =  f  (V^+.  .  .  +<^).  The  isomorphism  B  =  B^  is 

defined  by  the  commutative  diagram  and  it  is  called  a  transmutation  operator. 
We  remark  that  B  is  an  isomorphism  of  convolution  algebras,  as  such  it  has 
been  originally  studied  by  Bochner,  Leray,  Levitan  111],  117]. 

It  is  easier  to  find  first  B  *  explicitly.  Namely,  let  p  be  a  smooth 


radial  function  of  compact 

support  in  R  , 

then  if  r 

00 

r 

CO 

S9  (t)  = 
n 

-itxj  , 
e  dx  ^ 

pfr)dx2 

J 

1 

-oo 

J 

-oo 

-oo 

Comparing  with  the  above  diagram  we  see  immediately  that 


00  03 


V>(Xj )  =  B 


1M(x1 )  = 


p(r)dx„. . . dx 
2  n 


“00  -00 


,("-0/2  r  r2 
2  Jn 


22\  n-2  , 

,+y  )y  dy 
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.  2,(n~1)/2  f 

“  r(^)  J, 


p(s)s(s2-x2)(n'3)/2ds. 


We  Introduce  a  change  of  variable,  so  that 

. { n- 1 ) /2 

n- 
~2 

o  r> 

Let 


(n-l)/2  f" 

*><V  =11PT— 

r(V-)  J. 


,  /-w  2  2 « ( n-3)/2  , 

p(vs)(s  -Xj)  ds. 


‘1 


2  2 
$(Xj ) ,  v  =  Xj ,  then 


(n-1 )/2  f 
n _ 

T(~)  . 


4>(v)  = 


M(v's)(s-v)(n  3)/2ds. 


To  invert  this  expression,  we  distinguish  two  cases.  First,  the  case  n  is 
odd,  n  =  2£+  1,  £  >  1.  Then 


n  Li.  J  JV=S 


Hence 


(2. 18) 


p(s)  =  (B 


2£+ 


x=s 


2 


03 

(t-v) 

L  J  1 

*  V 

The  case  of  n  =  2£,  £  >  1,  is  a  bit  harder.  First  we  find,  using  Laplace 
transforms,  that 

p(\/s)  = 

.  UVJ  I  I  1 

L'v  -*v=s 

Returning  to  the  original  function  <p  we  get 


(2.19)  p(s)  =  (B^^pUs) 


x=s 


f  1  d  T* 

r  00 

,  2  2, -1/2  ,  ,  . 

(t  -x  )  ^(T)xdT 

2x  dxj 

■  X 

We  can  now  apply  these  formulas  to  find  the  distributions  p  =  p 


n,  k, a, p 


A  .  1  3 


defined  by  (2.16). 

Lei  v  «  1/c^  :*  2l’v^r(y+^),  then 


Vpt)  «  c  vp 


-2t> 


-its,  2  2,u-( 1/2)  . 

e  (p  - s  )  <is. 


-P 


Hence,  if  J  (pa)  =  0,  we  rewrite 

V 


Vt]  _  i  rJ/pt)  J/pt)i 

,  2  2  2a  t-a  "  t*a  ' 

t  -a  >-  J 


ju(pt) 

t-a 


c  p 
v 


-2v 


-its  -ias  _  _  , . 

e  -e  ,2  2.i»-(l/2)  . 

- i - (p  -s  )  ds, 

t-a 


~P 


and 


J  (pt) 

v  -2v 

— -  =  c  p 

t+a  v 


e  ltS-eiaS,  2  2,i>-(l/2)J 

- - (p  ~s  }  ds- 


-P 


Therefore 


Vpt) 

t2-a2 


c  p 
v 


■2y  -p 


2a 


e  itS0  (s)ds, 
a,  v 


where  0^  ^(s)  *s  ^he  even  function  which,  for  s  >  0,  is  given  by 


(2.20) 


0  (s )  -2 
a,  y 


,  w  2  2,v-(l/2). 

sin  a(s-u)(p  -u  )  du. 


To  find  p  we  apply  the  transmutation  operator  to  tl.e  function 

~2v 

c  p 

V 

If  =  - - -  0 

a,  v  2a  a,v 


For  instance,  for  n  =  2, 
(2.21) 


r1  'i 

1,2.k.a,p(r)  '  ",k  2  3tIO.pl(r) 

C.  71  , 

,  2  2, -1/2 
(t  -r  ) 

f1  c2,k-( 1/2)  ,  _,jr 

(l-f[  )  cos  a(x-p£)d£ 

r 

T/p 

dx. 
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Finally,  we  recall  a  staple  property  of  the  Bessel  functions. 


^emma  7.  Let  a^.a^.a^  >  0  and  m  €  2  .  Let 


0(z)  =  Jn/2(alz)Jn/2(a2z)J(n/2)*m-l(a3z) 

Then  there  are  constants  L  ,  A  >0.  such  that  for  any  integer  t  *  L 

mm  m 

there  exists  p ^  f 1 1  so  chosen  that  if  either  |z|  =  or  |Im  z|  >  1 


(2.22) 


,  .  ,  .  .  .-( 1/2) (3n+l+2m)  (ai+a2+a3)  I  lm  z| 

|  e  ( z  )  |  >  A  1 2  I  e 

m 


Here  n  >  2  is  considered  fixed  throughout.  If  we  are  given  6  >  0  in  such 
a  way  that  a.  €  (S,S_J]  then  the  constants  A  ,  can  be  made  to  depend 

v) 

only  on  6  and  m.  They  actually  become  explicit  in  the  proof. 

Proof.  The  function  e  is  an  even  entire  function,  therefore  to  prove  (2.22) 
we  can  assume  we  have  Re  z  >  0.  In  this  case  we  can  use  the  following  asymp¬ 
totic  development  of  the  Bessel  function  J  (u  >  1)  (see  [13]) 


r — —  / —  1 1  m  z  I 

(2.23)  |J  (z)  -  /  —  cos(z  -  ^(2i'+l ) )  I  <  3/  ^ 


f|4r2-l|  fn  |  4i'2-  1  I  ’J ") 

l  8  '  eAPl  8'z'  JJ 


Clearly,  when  |zl  >  gl4t'2-l|  we  obtain 


r~  r~  „  I  Im  z  I 

|j  (Z)  -  /?_  cosiz  -  ;(2u*l))|  <  p  /||4u  -  II— T 
v  nz  4  82  jz|3/2 


On  the  other  hand,  if  V  =  {(21+l)n/2  :  f  €  Z),  d(z,V)  =  mi n{  1 , d i st ( z ,  V) } 
then  the  cosine  satisfies  the  Lojasiewicz  inequality 


,  .  1  ,,  1  Im  z  I 

I  cos  z |  >  — d(z, V)e 
ire 


Considering  the  expression  (2.23)  we  see  that  we  need  to  know  how  many  zeros 
does  the  product  cos(ajZ  -  ^( n+1 ) )cos( a2z  -  ^( n+1 )  )cos(a3z  -  ^( n*2m- 1 ) )  has 
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In  an  interval  U. fMl;  as  we  have  at  aost  such  zeros,  we  can 

find  a  value  e  ]/./*l{  such  that  the  distance  from  p ^  to  any  of  these 

zeros  Is  at  least  5  - - — .  Hence,  on  the  circle  |zl  =  p„  we  have 

c.  n+aj+aj^aa  * 

Icosla.z- J(n*l))|  *  1 —  eailIm  z< 

1  4  2e  ( n+aj  ■*a2+a3) 

On  the  other  hand,  if  |Im  zl  >  1,  then  If  o  =  ninll.a  } 

J  J 

1  1  71 1  , 1  \  \  1  -v  a,  I  Im  z  | 

cos(a,z  -  7(n+l ) )  >  —  e 

14  ne 

that  is,  the  previous  inequality  is  valid  both  for  lz|  =  p^  and  for  |lm  z| 

i  1.  Therefore,  on  every  such  circle  as  well  as  on  |lm  z|  >  1  we  have  that 

aj | Im  z I 

|J  (a  z)|  >  - *£■ - — 

2ev/2n(7r+a1  +a2+a3)  | a^z  |  1 

2  2 

if  we  assume  |z|  >  ^  e  Tt(n'fa,i*£^+a3_I.rL,  The  same  kind  of  estimates  hold  for 

4  ajat 

the  other  two  factors  in  6.  Let 


(2.24) 


3  2  2  2 

Lm  =  4  e  n^n+ai  +  a2+a3)(n+2m)  /min  {a^,a^,a  ) , 


then  for  l  >  L  ,  |z|  =  p  for  llm  z|  >  1,  |z|  >  L  we  have 

m  t  m 


(2.25)  |0(z) |  > 


- ( n+ 1 ) /2  - ( n+ 1 ) /2  ~((n+l)/2)-m  (aj +a2+a3) | 1 m  z| 

jl _ J. _ ^3 _ ttltt2a3e _ 

t  •->_ ,/tT t  _  , „  .  1 ,3,  .  1  /2 ( 3n+2m+  1  ) 


[2ev/2n(?r-*a1  +  a2+a'i ) )'  Iz 


Remarks .  1)  The  sequences  p^  depends  both  on  the  a  and  on  m. 

2)  We  do  not  really  use  the  full  strength  of  the  estimate  (2.23).  It 
would  be  enough  to  use  Weber’s  estimates  [19,  §  7.33). 

We  conclude  this  preliminary  section  with  upper  bounds  on  the  Bessel 


functions,  which  depend  on  the  fact  that  the  order  v  is  either  an  integer  or 
half  an  integer.  From  [19,  §2.2  and  §3.32)  we  have  for  Be  z  >  0  and  q  an 
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Integer  fc  0 


where 

that 


P 

q 


J  (z)  *=  -  cos(q0-z  sin  0)do 
q  ii 

J0 

Vll/2),Z)  -  ('1)qNqf  elZC°S  6pq(C°S  6)  Sin  6  dG’ 

J0 

is  the  Legendre  polynomial  of  degree  q.  From  [14,  8.917)  we  know 


IP  (cos  6)1  <  1. 

q 


Therefore 


(2.26) 


IJ  (z)|  <  e 

q 


1 1  m  2 1 


1 Jq+(1/2)(Z) 1 


<  /I 


.1/2  I  I m  z I 
z  |  e 


1  /2 

which  has  the  advantage  of  being  a  global  inequality,  but  loses  |z|  and 
|z|  respectively,  from  the  asymptotic  expansions.  This  is  a  minor  loss  for 
us,  and  we  will  use  this  estimate  below. 


3 .  Reconstruction  of  a  function  in  the  case  of  two  disks. 

Let  2  =  :  £  >  0,  7)  >  0,  J  (£)  =  J  .-.(tj)  =0).  It  is  known  that 

n  n/2  n/2 

if  r.  >  0,  r„  >  0  and  r,/r_  t  Z  .  then  the  Pom.ru  iu  transform 
12  1  2  n 


P  : 


(* 


B(0,r  ) 


■r.* 


B(0 


P2) 


is  injective.  Moreover,  this  condition  on  r^/r^  ’s  also  necessary  for  the 

injectivity  112).  In  case  i, /r„  is  badly  approximated  by  elements  of  Z 

l  z  n 

that  is,  in  case  there  are  constants  c  >  0,  N  >  0  such  that 


Therefore,  under  the  arithmetic  assumption  13.1)  any  function  (or  distribu¬ 
tion)  f  in  Rn  can  be  reconstructed  from  the  knowledge  of  its  Pompeiu 

transform  Pf  =  (*D,n  ^.Zorn  x»f).  This  process  (3.2)  is  known  as 

d i u , r  j  j  d(u, r2 ) 

deconvolution.  In  case  we  only  have  a  function  f  defined  in  a  ball  B(0,R), 
with  R  >  r  +  r2>  even  if  we  use  the  deconvolvers  vyv2’  we  wi  ^  only  be 
able  to  reconstruct  f  in  B(0, R-fr^+r^)  )  •  For  this  reason,  we  will  give  a 
different  kind  of  explicit  inversion  process  for  tnis  local  problem.  It  turns 
out  that  the  arithmetic  condition  (3.1)  will  not  be  necessary,  only  rj/r2  * 

Z  .  In  particular,  this  answers  on  the  affirmative  a  question  posed  by 
Zalcman  121]  for  the  case  R  =  oo.  In  fact,  our  procedure  is  basically  a 
method  to  approximate  f  with  an  estimate  of  the  error.  In  this  sense,  we 
are  dealing  with  several  possible  procedures  that  depend  on  the  way  we  want  to 
represent  a  given  unknown  signal  f.  For  the  sake  of  concreteness  we  will 
describe  this  in  detail  when  we  represent  f  in  terms  of  the  Fourier  expan¬ 
sion  in  B(0,R)  by  spherical  harmonics.  In  a  forthcoming  paper  we  will 
explain  the  modifications  necessary  if  one  wants  to  use  wavelets  or  the  Radon 
transform  on  spheres. 

OP 

To  begin,  let  us  point  out  that  if  f  is  not  C  in  BlO.R),  then  we 


can  first  consider  f*^,  supp  £  B(0,c),  0  <  c  «  1,  and  since  Pf  is 

given  by  convolutions  we  will  have  P(f»y  )  ~  Pf«y  .  The  first  component  of 


Pf,  *g(0  r  j»f.  is  known  in  B(O.R-rj). 


hence  that  of  P ( f • y ^ ) , 
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v  .•f*«  ,  will  be  known  In  B(O.R-c-r, ).  Similarly  for  the  second 

*B(0.rj)  rc’  1 

component.  Hence,  replacing  R  by  R-c,  and  then  taking  c  so  small  that 

CD 

R-c  >  rj  *  r2>  ^  can  &ssume  f  is  C  . 

Given  a  function  f  €  Cc°(B(0,R)),  wc  have  a  Fourier  expansion  of  the 


(3.3) 


f(x)  = 


W  - 

a  .  (r)S  (wl 
/  .  m, k  m, k 

m=0  k 


)  ,  x  =  rw.  w  €  S 


where  {S  )  is  an  orthonormal  basis  (with  respect  to  the  normalized 
m,  k  k 

measure  of  the  unit  sphere)  of  the  spherical  harmonics  of  degree  m.  That  is, 


r  S  (w)  =  H  .(x)  is  a  harmonic  homogeneous  polynomial  of  degree  m.  The 
m ,  k  m ,  k 

dimension  of  the  vector  space  SpanHS^  ^1^)  is  0( m  2)  [18).  The 


coefficients 


a  (r) 
m,  k 


.  cn_1 


f(ru>)S  .(wider  (w) 
m,  k  1 


satisfy  an  inequality  of  the  form 


a  ( r )  <  C  — 

m,  k  N  N 

m 


for  any  N  >  0.  Another  way  to  write  down  this  formula  is  as  the  action  of 

the  distribution  r  m<r  (x)  on  a  smooth  function 

r 

(3.4)  a  ,  (r)  =  <r  "V  (x),  H  ,(x)f(x)>  =  <r  mH  ,(x)<r  (x),f(x)> 
m,k  r  m,k  m,  k  r 


.m  2-n-m 

- <h  ,  J-  T  .(Ixi).nxi 

-'(.-Din  "■l‘exJ  r'” 


where,  as  given  by  Lemma  6, 


T  (lx|)  =  (r2-  lx|2)n’“1*Rfn  Jx). 
r,  m  8(  0,  r  J 
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the  idea  now  will  be  to  write  T  as 

r,  a 

Tr.»  *  gl*zB(0,r1) +  82*2B(0.r2)* 


with  control  on  the  supports  of  gj  and  g2-  In  this  way  we  reduce  the 
problem  to  a  decomposl t Ion  of  radial  functions.  We  fix  some  notation  first. 

Choose  a  strictly  Increasing  sequence  of  positive  numbers  »:  with  limit 

D 

— - — -  l.  and  a  corresponding  increasing  sequence  of  radii 
rl  +  r2 

R*  =  (r1+r2)ll+ck)’  k  ~  1(  R0  =  °' 

We  want  to  decompose  the  radial  distribution  T  ,  and  th  .*  decomposi- 

r ,  m 

tion  will  depend  on  the  value  k  such  that  r  e  We  eca11  that 

the  Fourier-Bessel  transform  of  T  is  given  by 

r ,  m 


T  (t)  =  (2n)n/22n'_1 
r,  m 


(m-1 ) !  r 


n+2m-2 


j(n/2)+m-l 


r  ( t ) . 


Proposition  8. 
For  every  r,m 
tri but  ions 


Let  r1>r2  >  0,  «  Z^,  +  r2  <  R>  fixed  as  above, 

with  -  r  <  there  are  two  sequences  of  .'adial  dis- 

v whose  Fourier-Bessel  transforms  satisfy 


<3  51 

c'  -(n/2)-m+(3/2)„  .. -( n/2  ) -m+ ( 13/2 )  Rk|Im  z| 
^  j-r  I!  zli  e 


f  or 


l  > 


2 

cm 


(3.6) 

(3.7) 


lp£(z)l 

l^(z)| 


^  „  -(n/2)-m+(3/2)..  „-m+2  (r;>+3k)|lm  z| 

<  c"r  II z |!  e  “ 

<  c„r-(n/2)-m+(3/2)||^jj-m+2e(r,+Bk)  |  Im  z| 


In  this  statement,  we  use  the  notation  ||z||  =  max{  1.  |z| } ,  c,c' 
strictly  positive  constants  depending  only  on  r  fr2,R,c  and  n. 


being 

while  c " 
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depends  also  on  l  and  a.  The  quantity  0fc  denotes 

Proof .  We  are  going  to  apply  Lemma  7  with  aj  «  rJt  a2  *  r2>  a3  *  0fc.  The 

2  2 

corresponding  value  L  i  ^  e  S  cm*Vcjv  We  let  q  be  an  even 

vw 

integer  to  be  fixed  later  on.  To  simplify  the  notation  w;  denote  f ( z )  = 
J(n/2)+m-l(rz)' 

The  Idea  of  the  proof  Is  to  profit  from  the  simple  fact 


(3.8) 


1 

2ni 


ICI  =p, 


ftp 


<K  = 


ff(z)  if  |z|  <  p 
0  if  izl  >  p 


l. 


The  usual  method  of  interpolation  theory  consists  in  rewriting  (3.8)  as 
fol lows 


(3.9) 


f  ( z )  = 


1 

2ni 


f(p  Kqe(U-zqe(z)) 


IP=p, 


<qe(0 


C-2 


<K 


<ar2q9(z) 


lci=p, 


f.,r  Izl  <  p, 


<qe(p 


and 


(3. 10) 


2ni 


f(p  Kqe(O-zq0(z) 


Kl=p, 


<qe«) 


C-z 


<K 


zq0(z) 


2n  i 


KMp, 


P-  for  |2|  >  p 
<%«)  ?"2  ' 


We  define  an  even  entire  function  by 


(3. 11) 


Hl(z)  =  2tri 


f(p  (cqe(<)-zqe(z)) 


I<I=p 


_  Oe(0 


<-z 


dC. 


that  is,  the  interpolation  of  f  on  the  zeros  of  zq0(z).  Note  that  is 
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even  Independently  of  the  parity  of  the  Integer  q;  we  took  q  even  to 
simplify  the  computations  below,  but  this  Is  not  necessary. 

We  are  going  to  show  the  family  is  a  bounded  family  in  S'  (R). 

For  that  purpose  we  recall  the  upper  bounds  for  f  and  0  that  we  obtain 
from  (2.26), 

(3.12)  |f(z)|  <  err_((n-3)/2)_n'||zir{ln“3)/2)'n'er|lm  zl 

(3.13)  |0(z)|  <  eRk(r  r  )"(n” 1 }  /2g“ ( ( n“3  5  /2  )_m||z|(~ <3(n" 1  1/2  ]  ~m*  ]  eRk  1 1 m  2  I 

I  2  k 


We  note  that  for  n  even  these  two  estimates  are  too  big  by  a  power  ||z||1/2 

Rk 

car.  be  replaced 


3/2 

and  tl z ||  respectively.  Also  the  constants  er  and 


R  2  2 

by  e  .  Choose  /  >  cm  /«-£  >  >  2,  then  for  Kl  =  Pf  we  have 

|e(<)|  >  c  )/2)~m  -( 1/2)  (3n+2m+l  J^R,,  1  Im 

lk  ( 


where  =  Cjfr^^.R)  is  given  by 

(r^)  (n+1)/2»in{l,r1,r2,Pk>3 

c  -  - - . 

(2eV5n(n+R)) 

For  Izl  >  p£  we  can  use  (3.10)  to  estimate  H£.  We  let  q  =  n  3  if  n 
is  odd,  q  =  n  +  4  if  n  is  e\~n,  then 


|H£(z)1  < 


r  ( (n-1 )/2)+m 
e  Pk 


c 


1 


-( (n-3)/2)-m 
r  P 


-q+n+3 

e 


I  z  I q  1 0  ( ;• )  | 
Izl-P/ 


s  r  r 


~( (n-3)/2)-m  |z|-(n/2)-»*(n/2) 

UhF,  e 


where 


OD  _  ^ 

e  rjr2Pk(2ev^n(n+R) ) 

c2  -  - - 3 -  5  c2(rl,r2,R) 

min(  1 ,  Tj , r2*  Pj  > J  2 
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(we  have  taken  into  account  here  the  fact  that  the  upper  estimate  for  8  can 

3/2 

be  improved  by  | zfl  when  n  *  even).  For  |z|  <  p(  we  use  (3.9),  so  that 


(3. 14)  |Hj(z)|  <  I f (z) I  ♦  c^r 


-( <n-3)/2)-m  >z||~ (n/2  ]'m* ( 1 1/2  }eRk  1  1  m  z| 

P^~  I  z  I 


These  estimates  are  fine  when  |z|  <  £-1  or  |zl  >  £  *  1,  since  then 
||z|  ~Pf\  *  1.  In  this  case  we  have 


(3.  )5)  |H{(z)l  <  zl 


In  the  remaining  annulus,  £-  1  <  |z|  <  £  +  1,  we  can  use  the  maximum  prin¬ 
ciple  to  obtain 


(3. 16) 


|Hj(z)|  < 


C2r 


-( (n-3)/2)-m. 


, - ( n/2 )-m+ (11/2)  Rk|z| 
e 


The  new  constant  c"  is  not  substantially  different  from  c^.  and  it  is 


.  4  ,  .  ..  nil  fllzll)  (n/2  )+m-(  11/2)  . 

independent  of  m  since  if  m  +  ^  -  j? bey  - 

<  c^t  Independent  of  m  since  £  >  cm2/c2.  The  other 
values  of  m  are  handled  by  the  same  reasoning.  We  conclude  that  (H^  is 
bounded  in  Exp  ( C )  =  (entire  functions  of  exponential  type).  To  conclude 
that  (H^  is  bounded  in  £'(1R)  we  need  estimates  on  the  real  axis.  From 
(3.14)  we  know  that  if  x  €  F\{£-1,£+2J  then 


(3. 17) 


-( (n-3)/2)-m  „-(n/2 ) -m+ ( 1 1/2) 

IH^(x)  I  <  2c^r  l.xl! 


As  it  was  done  in  [2]  we  write 


H^(x)  =  H^(x)  -  H^+3(x)  +  H^+3(x), 

hence,  by  (3.17),  the  last  term  satisfies  the  correct  inequality  when  x  € 
[£-l,£+2],  and  we  only  need  to  estimate  the  difference  (H^-H^^Mx).  Let 
F.  be  the  Jordan  curve  which  is  the  boundary  of  the  region  (Re  z  Z  0, 

V 
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s  | z  |  S  1 1 »  z I  S  1}  and  Tg  its  symmetric  with  respect  to  the 

Imaginary  axis.  They  are  chosen  with  a  convenient  choice  of  orientation  so 
that 


(3.  18) 


Hj(x)  -  H^+3(x) 


4  [  1  f  (p(gqe(Q-xqe(x) ) 

2ni[J  J  J  <qe(OK-x) 

1 1 


For  £  €  and  x  €  I  £-1,1+2],  we  have 


Cqe(Q-xqe(x)  max  jd  (wqe(w))| 

^  •*  dw 


C-x 


|w-x|<5 


„  ,  6R,  . - ( n- 1 ) /2  - ( ( n-3 ) /2 ) -m .  . - (n/2 )-m+ ( 1 1/2 ) 

<  c3e  {r,r2)  Ixl  ; 


by  the  same  reasoning  that  led  to  (3.  16)  the  constant  is  independent  of 

m.  For  C,  €  and  x  €  [l-\,l+2],  we  can  use 

l<qe(0-xqe(x)  I  <  Kqe(C)l  ♦  lxqe(x)l 

^  ,  6R,  n-1 )/2  - ( ( n-3 )/2)-m .  , -(n/2)-m+ ( 11/2) 

<  2c' e  (r,r_)  p.  lx! 


We  conclude  that  for  x  €  l£-l,£+2]  we  have 


I  H^t  x ) |  <  c^r 


-(n-3)/2-m„  ,, - ( n/2 ) -m+  (1 1/2 ) 


which  is  therefore  valid  everywhere.  By  applying  the  proof  of  the  Phragmen- 

4 ,  „  ..  ,  i  .  (  n/2)  +  m-(  11/2)  iRkz  .  s  .  . 

Lindelof  theorem  to  the  function  (z+-)  e  H.(z)  on  the 

Ul  t 

quadrant  Re  z  >  0,  Im  z  >  0  we  see  that 
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IHf(z) 


<  c.r 

4 


-( (n-3)/2)-m..  - ( n/2 )-m+  ( 1 1/2 )  Rk!lm 

llzll  e 


where  c^  is  independent  of  m.  This  argument  works  in  the  other  three  quad¬ 
rants  This  shows  that  this  sequence  {H^  is  bounded  in  £'(IR).  Moreover, 
we  can  show  its  limit  is  f  and  one  can  estimate  the  error  by  the  previous 
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procedure . 


If  |z|  <  p formula  (3.10)  yields  (see  (3.14)): 


(3. 20)  IH^(z)  -  f (z) I  <  c^r 


..  ..-(n/2)-m*(ll/2)  Rk|Im  z| 
-((n-3)/2)-m  |[z|| _ e_^ _ 

Pt- Izl 


<  2_2r~(  (n-3)/2)-mj(z||-(n/2)-mM13/2)eRlll  1m  z| 


if  |z|  <  //2.  For  |  z  |  >  £/2  we  have  from  (3.19)  and  the  estimates  for  f 
that 

(3.21)  | H£( z >  -  f  ( z )  I  <  |H^(z)  I  ♦  If  (z)  I 

C4  -( (n-3)/2)-m.,  .-(n/2)-m+ ( 13/2  )  Rk|lm  z| 

<  r  z  e 


c'  is  independent  of  m  and  bigger  than  2c  ,  hence  the  last  estimate  is 
4  <- 

valid  everywhere  as  stated  in  (3.5). 

We  can  describe  and  v ^  explicitly  using  the  Residue  Theorem.  Let 

us  do  it  for  the  case  the  sequence  (c^)  has  been  chosen  with  the  additional 

requirement  that  the  pairs  jn/2^rjZ^  anc*  J(n/2)+m-l(Pkz)  (J  =  !’2)  have 
no  common  zeros.  We  note  that  if  there  were  a  common  zero  the  multiplicity  of 
this  zero  for  the  product  of  the  two  functions  would  be  exactly  two  and  the 
formula  below  can  be  modified  easily.  From  the  definition  of  we  can 

write  for  example  when  n  is  odd 


*Yz)  ~2  ^n/2(r!z)')n/2(r2z>e3.*(zl  +  2  ^n/2(  r1z)J(  n/2  )  +  m- 1 ( Pkz  )e2.  f(  z  1 


where 


V  _ fU^ _ Jn/2(riz) 

gl,f  z  l_  n+3 ,  ,  .  .  to  \  ri J'  ,~(r  a) (z-a) ’ 

M<P,  *  Jn/2(r2°,J(n/21.m-l(Sk<‘)  ’  n/2  1 


',n/2(rl°l  0 
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and  analogous  formulas  for  g^  ^  and  S3  P  is  the  polynomial  of  degree 
S  n  ♦  2  (In  fact,  P  must  be  an  even  function)  given  by 


no  - — — ] 

Cn  3e«)  i 


We  define 


^(z)  2  ^n/2(,‘2z)g3.<(z)  *  2  J(n/2)  +  m-l(,3kZ)e2./(z) 

vt(z)  =  z  j(n/2)4m-l(Pkz)gl^(z)  +  P(z),j(n/2)(rl2)J(n/2)  +  m-l(Pkz) 


Hence 


(3.22) 


Hf(z)  =  ^(z)Jn/2(rlz)+D<(2)Jn/2(r2z)' 


and  the  estimates  (3.6)  and  (3.7)  follow  from  the  asymptotic  estimates  of  the 
Bessel  functions.  Even  though  we  need  to  replace  n  +  3  by  n  +  4  in  case 
n  is  even,  the  estimates  remain  the  same  by  the  remark  made  after  (3.13).  □ 

The  previous  proposition  leads  to  the  following 

Theorem  9.  Let  r^,^  >  0,  r  /r  «  Z^,  +  r2  *  R>  <ck*k>1  a  strictJy 

R 

increasing  sequence  of  positive  numbers  with  limit  - -  1,  R  = 

rj+r2  k 

(i~l+r2)  ( 1+e^)  (k  >  1).  Rq  =  0.  For  any  r  >  0,  re  [ P^_5 .  R^[  ,  and  S  a 

normalized  spherical  harmonic  of  degree  m  there  are  two  sequences  of 

distributions  t)^,  8^  of  order  <  n  +  3  with  compact  support  in  B(0,R-r^) 

2 

(respectively  B(0,R  r^))  such  that  for  (  >  cm  and  any  function  f  in 
C°°(B(0,  R)  ) 


(3.23)  | 


f(r,)S(.W»|-<Ir,B(  0.r,),r>'<*f*B(  0.r,)'f>' 
,n-l  1  2 


7  _,-N_-(n-3)/2  ,5 


<  'j  (R-D  r 


max  | - f  ( x )  | , 

la  UN  5xa 


|x|<R^ 
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i  2  1 

where  N  ■  *  ^43*Sc  and  »  >  0  Is  a  constant  depending 

on  Tj ,r2>R, n.Cj. 

Proof ,  Let  mil  and  let  H  be  the  harmonic  polynomial  homogeneous  of 
degree  m  such  that  H(pu)  =  p^Sfw),  w  €  Sn  p  >  0.  From  (3.4)  we  have 


C 


2-n-m 


f  (rw)S(w)d<r(w)  = 


.n-1 


2m  1 ( m- 1 ) ! fi 


<T( I x | ) . H  ^  f(x)>. 


2  2  in _ 1 

where  T(|x|)  =  (r  -|x|  )  *D,n  We  have  seen  that  the  Four) er-Besscl 

d  ( u  t  r ) 


2-n-m 


transform  of 


2m-1(m-l)!fi 


T( |x| )  is  given  by 


.  _(n-2)/2r,n.  m .  .  .  . 

th~ *2  F  2  r  J(n/2)-nn-l  rt  ' 


We  can  apply  Proposition  8  and  find  explicit  radial  distributions  v ^  with 

),  supp  v £  <  BCO.R^-r^)  such  that  the  radial  distribution 
whose  Fourier  Bessel  transform  is  given  by 

Vz)  =  2^n_2^/,2n^)rm|j(n^2 )+(n_j f rz)  "  ^n/2^rlz ~  2(n/2 )  * r2z z ’) 

satisfies 


supp  p£  £  BCO.f^-rj 


'Vz)l  -  7 


c  - ( n-3  )/2..  ,,-(n/2)-m+(  13/2)  Rkllm  z| 


(f  >  const.  m  ). 


If  we  define  the  distributions 


V»>  -  l^HfMMf«l*l) 

(  n  laxj  ( 

1  n 


.  ,  ,  c-nVt.fa  )  , ,  , , 

V*1  "  -r— -"IrK'11'" 

rznn  1  1 


then 
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f  f(ru)S(U)d<r(u)-<«t.*B( o  *B(0 ,r ,,’f>  *  ‘VH(li]f> 

Jsn-1  1  2 

We  will  now  estimate  <x<>H^jf>  using  Plancherel’s  formula.  Let 

J  1  2 
be  a  smooth  radial  function  which  Is  Identically  one  on  the  ball  B(O.^R-*=R  ), 

-  2  1 

with  support  in  B(0,  .  0ne  can  find  them  ln  such  a  that 

max  <pR(x)|  <  1  q(R'V~q  -  yq(R'r)’q- 

x,  I  a  I =q  5xa 

The  constants  y  are  Independent  of  R  and  k  Then 

q 

|<x<iHgLjf>|  =  l<VH[^](V)>l  5  ^UCI)lll(5)ll(vkf)‘(C)lcl? 


c  -(n-3)/2  f  „  .  ( n/2 ) 
jr  |  I! pi! 


*(11/2)  |S(«)|  I  («>kf )  (pulldo-(w)dp 


„  c'  -(n-3 )/2 

-  rr 


On  the  other  hand 


|S(w)|  dcr(u) 


_  GO 

1/2  ■ 

J0 


llpll  ^n/2^  +  (  1 1/2^  max  l(vkf)  (pu)ldp. 
weSn_1 


(£)l  i  y'll^H  N(R-r)  N  max  l^—f(x)|. 
k  "  | x I  SR'  ox 


where  N  >  and  R^  =  §^*5^' 

This  concludes  the  proof  of  Theorem  9  for  m  >  0.  For  m  =  0  everything 
works  the  same  way,  except  that  T  is  now  the  normalized  measure  of  the 
sphere  of  radius  r,  with  Fourier-Besse  1  transform  equal  to  j(n_2)/2^rl^  up 
to  a  constant.  o 


The  estimate  (3.23)  leaves  us  with  a  little  problem  for  small  values  of 
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r.  In  the  case  0  &  r  <  R-  (rj+rg).  **  can  usc  ***  «*pHcU-  formula  to 
approximate  6  In  the  form 


60  *  )'al  +  ZB(0,r2)*a2’ 


a  radial  distributions  of  fixed  ordfr.  supp  £  B(0,r2).  supp  a £ 
J 

B(0,r  ).  Hence  we  compute  directly 


f(x)  E  I°l,ZB(0.r1)*fl  4  Ia2**B(0.r2)*fl(x) 


for  1x1  <  R-  (r^+r^,  the  Precise  estimate  of  the  error  being  of  the  form 
as  in  (3.21)  (without  the  dependence  on  r  and  R-r);  we  use  as  a  function 

0  just  the  product  ^n/2^r \z^n/2^r2Z  * ' 

If  we  choose  the  sequence  in  Theorem  9  so  that  every  zero  of  0(z)  = 

Ji1/2(rlz)Jn/2(r2z)'}(n/2)+m-l(PkZ)  iS  slmPle>  the  distributions  u(,  and  a 
fortiori  11^  and  5^  can  be  made  completely  explicit.  We  only  need  to  know 
the  explicit  form  of  a  radial  distribution  whose  Fouri er-Bessel  transform  has 
the  form 

J(n/2)4m-l(_^>  (m  >  n 
z  -a 


when  a  is  a  zero  of  J(n/2)+m_1 <P* >  • 


These  distributions  have  been  computed  in  section  2;  recall  for  example 

J(n/2)+m-l (pt  ] 


that  when  n  =  2,  p  (t)  = 

m,  a,  p 


.2  2 

t  -a 


with 


*V«,p(r)  =  '^T5  *lo.pl(r) 


,  2  2, -1/2 
(t  -r  ) 


, i  *2.m-( 1/2)  ,  . 

( 1-C  )  cosa(-t-p£)df; 


iJr/p 


dx . 


The  sequence  of  approximations  for  the 
f(reie)e±im0d0,  m  >  1,  re  (R^.R^ 


Fourier  coefficients 

can  be  written  as 
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where 


(<'‘r,»,«'XB(0.r1)*f>  4  <*r.fc.rZB(0.r2)*f>)tew* 


II 


(-1) 


m  m 


r,  m,  l 


2nr, 


r  fa  a  V" 
[ax  "  3yj 


(-1) 


m  m 


[3x  3yj  l 


‘  I 

r,  m,  l 


2nr. 


and 


lB(0,r2) 


«5(U* 


J  fra) 
m 


"r2  0<£"r(  «\j;(r2«)J1(r2«)j;(Bk<.1 


W>-° 


(B2-M2)m"3 


21"  1  (m-1 ) !  n 


"B(0.pk) 


I 


j  (ra) 
m 


0<«<p,  a  r2Jl(r!*)Jm(Sk“)Ji<r2“> 


l.a.r. 


Jl<r2c)=0 


V,  =  -A' 


02-l-l2)m_1 


2™  1 ( m- 1 ) ! n 


P?"  *B(0‘V* 


I 


j  (ra) 
m 


0<^(  “SrlJl(r2“,W)Ji(r.a)  ‘'“'H’ 


jJ(r1a)=0 


♦(wv2) 


,_2  .  .2. m-1 

*B(0.r,)  (Pk_l*  1  } 

•  - = - —  Xz 


„  2  m-1  .  ,.102m  B(0,fi 

2nr,  2  n(m-l)!B,  '  k 

1  K 


2  4 

The  polynomial  Sk  ♦  Tkz  +  Ukz  is  defined  as 
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sk 4  V 4  V 


-  Res  [j(rO 

<*0  L 


^oV.-c5 


t°J1(rJOJ1(r2Oj,iO|<C 


i 


4.  Local  Inversion  of  the  Pompelu  transform  for  the  case  of  a  square. 


In  this  section  we  will  show  how  to  reconstruct  a  function  f  defined  In 

2 

the  disk  B  =  B(0, R)  of  R  in  terms  of  its  averages  over  squares  of  side 
2a,  with  2V2a  <  R.  We  will  also  see  that  this  inversion  procedure  general¬ 
izes  to  practically  any  set  E  with  the  Pompeiu  property. 

The  idea  is  to  reduce  ourselves  to  the  situation  of  the  previous  section. 
Let  Q  be  the  square  [ -a, a] < l -a, a ] ,  Xq  its  characteristic  function.  Then 


(4.  1  ) 


T  = 


a2 

5  *q  _ 

Oxdy  ~  5(a,a)  5( -a, a)  + 5( -a, -a )  *(a,-a)' 


It  is  immediate  that  5?^T  =  0.  On  the  other  hand,  this  is  not  true  for  some 
derivatives  of  T.  We  compute  directly  the  Fourier-Bessel  transform  of 


* 


0 


'd^_ 

•Sz1 


-t|  .  Let  t  be  a  positive  number,  choose  <;  e  C^,  |<;|  =  t,  and  iden¬ 


tify  z  =  x  +  iy  with  (x,y)  €  IF  ,  then 

,k  ~  r^k 


h>d?T)] <u  -  <*oC-?T)‘z;-e“7' 


=  <^— rT(z) , K  (e~12*^)> 
dzk  ° 


=  < — j-T(z) ,  <r.  .  ( t  )> 
azk  |z| 


=  <--j-T(z),J  (t  |z  J  )> 
Sz 


=  (-l)k<T(z),~ J  (t|z|  )> 
3zK 
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(as  shown  In  Lemma  6) 


t^<T(2),  |  Jk(t|2|)>. 


That  is. 


,4'2)  feT)  tU  ‘  J?)\2V,(""llkj 


if  k  *  2  (mod  4) 
(v^ta)  if  k  s  2  (mod  4) 


In  particular. 


(4.3) 


(4.4) 


K(JT))ll)  * 
KS‘t))  ,u  =  ^ 


(t)  =  -ia2t4J2(v^at) 


t12jg( V2at ) 


The  two  functions  above  have  only  the  origin  as  a  common  zero  of  order 
4,  by  a  well  known  property  of  the  Bessel  functions  (see  (19,  p.  485)).  On 
the  other  hand,  the  function  (ftQ;tQ)  If  we  use  these  three  functions 

with  a  reasoning  similar  to  that  in  the  last  section  we  would  need  R  >  3v^2a 
to  obtain  an  inversion  formula  which  will  be  even  a  bit  more  complicated.  (We 
will  return  to  this  point  below.  )  We  sidestep  this  question  by  working  only 

fg2  -»  -j 

with  the  radial  distributions  and  — gT  to  first  reconstruct 

S dz  '  S3z  ' 


hi.  Namely,  consider  the  two  radial  distributions  Pj.p  such  that 


p  ( t )  =  -ia2j  (v^at ) 


hil)  - 


whose  Fourier  Bessel  transforms  have  no  common  zeros  and  which  are  supported 
in  B(0,V^a).  One  can  modify  Lemma  7  by  considering  the  function 


6(t)  =  p, (t)p_(t)J  (c2v£t). 

1  e.  ffl 


with  0  <  c  <  R  -  2v^a .  The  inequality  (2.22)  is  replaced  by 
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(4.5) 


le(z) I  *  U,-ll/2)(3*2.)e2vS(l«c)ll.zl 

Whenever  I z )  =  p'  or  lln  zl  2  1,  l  *  L' .  for  some  tl'.L'  >  0. 

t  m  m  m 

The  proof  of  Theorem  9  furnishes  distributions  such  that  for  any 

function  f  €  C*(B(0, R) ) 


'k 


f(re16 )e~lm6d0  -  <U^.p]»f>  -  <8^.M2»f> 1 


/  o  1  1 

<  \  (R-r)  max  | - : 

1  |a| <N  axa,aya2 


f(z)| 


where,  as  before,  the  distributions  depend  also  on  r.m.  have  support 

in  B(R-V^a),  their  order  and  the  constant  N  do  not  depend  on  r  or  m; 

2  i  _ 

R^  has  the  same  meaning  as  in  Theorem  9,  that  is  =  — R  ^  v^2a(  1+j:  ). 

2 

Applying  (4.6)  to  A  f,  we  can  compute  its  Fourier  coefficients  in 

2  6 

B(0,R)  in  terms  of  the  functions  ^  =  [^o(^~2T]]*f  and  *2  =  (^o[^~6T]]*f 


in  B(R-V2a),  since  A  p  = 


and  A  = 


By  Lemma  5,  we 


can  compute  the  functions  VyV2  in  terms  the  P°n>P<?iu  transform  of  f.  In 

p 

other  words,  we  can  (xplicitly  recover  the  Fourier  coefficients  of  A  f  In 
terms  of  P(f). 

We  need  now  to  show  how  to  compute  the  Fourier  corff i ci ents  of  f  in 
2 

terms  of  P ( f )  and  A  f.  This  can  be  accomplished  by  observing  that 

(9?0*q)  (°)  =  area(Q)  =  4a2 

and  that  (z)  has  an  asymptotic  behavior  of  the  same  type  as  that  of 

the  Bessel  functions,  as  It  was  shown  in  [3,  Section  6].  Hence  we  can  apply 


the  same  procedure  as  above  to  the  two  radial  distributions  ??  r  and  Ad 

CTQ  0 
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and  In  this  way,  we  can  approximate  the  Fourier  coefficients  of  f  In  B(0,R) 
as  closely  as  we  want  in  terms  of  the  "data"  P(f). 

Remark  1 .  The  procedure  for  the  local  Inversion  of  the  Pompelu  transform 
associated  to  a  square  Is  valid  for  a  very  general  class  of  domains  having  the 
Pompelu  property.  In  fact,  let  E  be  a  compact  subset  of  Rn  having  the 
global  Pompelu  property.  Let  us  assume  that  for  some  base  point  xQ,  the 
boundary  of  E  contains  a  hyperbolic  point  p  In  the  sense  of  (3, sect  Ion 
2).  One  can  assume  that  SBlx^.r)  n  E  =  (p>,  where  r  =  disKx^.p)  and 
E  £  B(x0>r).  This  is  Just  a  minimal  regularity  condition  on  the  extremal 
point  p  of  E.  By  Proposition  5  in  [3],  we  have  that  near  the  real  axis  we 
have  the  asymptotic  development 

(4.7)  (ft  *r)  (t)  ~  cos(  tr-i>  ) , 

x0  E  tVo  1 


where  are  some  positive  constants  related  to  the  geometry  of  dE  near 

p.  Outside  a  horizontal  strip  there  is  a  corresponding  good  lower  bound. 

This  behavior  is  the  same  as  that  of  the  Bessel  functions. 

Let  p,  be  ft  r_.  If  we  assume  that  R  >  3r,  let  cn  >  0  be  such 
1  x0  E  0 

that  3r  +  3c„  <  R.  For  any  >:  >  0  we  can  consider  the  set  E  which  is  the 
0  e 

translate  of  E  along  the  unit  direction  x^p  by  c,  then 


(4.8) 


(ft  Xr-  )  (t)  ~  cos  ( t  ( r+c ) -r  ), 

VEc  t1'0  ’ 


since  the  only  extremal  point  of  E^  of  interest  is  just  p + cxQp.  It  is 

then  immediate  that  we  can  choose  c  <  cn  so  that  p  and  (ft  x_  )  have 

0  1  x_  E 

0  e 


only  finitely  many  common  zeros.  Let  us  choose  with  such  a  property, 

and 


z . z  the  common  zeros  of  p  and  p  The  argument 

C  X  _.  L>  f  i  l*  1  ^ 


Mo  =  *  X 

0  ~c 


in  ( 3, Proposition  5]  shows  that  one  cannot  only  translate  along  the  ray  x^p 
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but  even  along  a  snail  open  cone  Z  with  vertex  0  in  such  a  way  that  the 


translate  of  p  remains  the  only  extremal  point  in  Et  r\  dB(x0>  |pM-x0l )  for 
t  c  Z.  By  Lemma  4'  the  corresponding  functions  have  no  common  zero  If 


For  each 


we  can  choose  distinct 


such  that 


<*  f(z.)  *  0. 

0  Tj  J 


hence  there  is  a  linear  combination  of  these  finitely  many  functions  which 


will  not  vanish  at  any  of  the  points  This  linear  combination  p^ 

will  automatically  have  a  good  asymptotic  development.  The  radial  distribu¬ 


tions  (with  respect  to  xn)  p  , p0,p„  are  supported  by  the  balls  B(x  ,r  ), 

U  i  6  U  1 

B(x0>r2),  B(xQ,r3)  respectively,  ri  + r2 + r3  <  R  and  the  method  of  the 
previous  section  applies  with  small  modifications.  Lemma  7  allows  us  to  com¬ 


pute  p  ,»f  in  terms  of  P(f), 
J 

any  ball  of  radius  R  >  3r. 


in  this  way  one  obtains  an  inversion  formula  in 


Remark  2.  In  comparing  the  previous  remark  with  the  construction  for  the 
square  one  notes  that  the  construction  for  the  square  could  be  further  simpli¬ 
fied  if  we  knew  that  and  have  no  common  zeros.  Though 

we  believe  this  is  true  we  have  not  been  able  to  prove  it. 
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We  have  demonstrated  the  feasibility  of  building 
super-resolution  systems  using  multiple  detectors.  The 
deconvolution  method  provides  a  real  time  linear 
implementation  of  the  reconstruction  problem  which  is  robust 
with  respect  to  noise  and  perturbations  of  the  overall 
system . 

The  next  task  that  is  starting  to  be  carried  out  by 
doctoral  students  of  Dr.  Berenstein's  in  Maryland,  Drs.  Gay 
and  Yger  in  Bordeaux,  and  Dr.  Taylor  in  Michigan,  is  to 
incorporate  postprocessing  features  adequate  for  particular 
tasks.  For  instance,  in  many  problems  we  need  not  only  to 
reconstruct  the  pixel  by  pixel  values  of  signal  to  a  desired 
degree  of  resolution  but  also  we  need  to  automatically 
identify  it,  e.g.  using  segmentations,  Voronov  diagrams,  etc. 
How  to  incorporate  this  image  (post)  processing  into  the 
deconvolutions  is  the  challenge. 

Another  thing  that  still  remains  to  be  implemented, 
since  it  requires  considerable  manpower,  is  to  write  a  user 
friendly  menu  to  design  and  simulate  multidetector  systems. 

It  is  for  this  task  we  feel  the  newer  completely  numerical 
algorithm  of  finding  deconvolvers  might  be  most  useful. 

The  local  Pompeiu  problem  that  "peeled-of f "  the 
deconvolution  problem  presents  new  challenges.  The  nature  of 
the  inversion  formulas  for  Chapter  4,  indicates  that  the 
deconvolution  and  local  deconvolution  formulas  must  be  well 


adapted  to  implementation  in  terms  of  wavelets.  If  this  were 
the  case,  it  will  also  help  to  produce  a  significant  data 
reduction  and  reduction  of  processing  time  in  ATR  and  similar 
problems. 

Finally,  the  purely  algebraic  problems  of  the  size  and 
degree  of  solutions  to  the  Bezout  equation  have  applications 
to  robotics,  control  of  distributed  parameter  systems, 
computational  geometry,  etc.,  and  they  show  the  unexpected 
payoffs  of  the  use  of  the  powerful  methods  of  several  complex 
variables  in  applied  mathematics. 


